THE CUT-OFF PHENOMENON FOR BROWNIAN 
MOTIONS ON SYMMETRIC SPACES OF COMPACT TYPE 

PIERRE-LOiC MELIOT 

Abstract. In this paper, we prove the cut-off phenomenon in total variation distance for the Brownian 
motions traced on the classical symmetric spaces of compact type, that is to say: 

(1) the classical simple compact Lie groups: special orthogonal groups SO(n), special unitary groups 
SU(n) and compact symplectic groups USp(n); 

(2) the real, complex and quaternionic Grassmannian varieties (including the real spheres and complex 
or quaternionic projective spaces when (? = 1): SO(p+q)/(SO(p) X SO(g)), SU(p+<?)/S(U{p) X U(q)) 
and USp(p + 9)/(USp(p) x USp((j)); 

(3) the spaces of structures: SU(n)/SO(n), SO(2n)/U(n), SU(2n)/USp(n), and USp(n)/U(n). 

In each case, we give explicit lower bounds for d-rvCMt) Haar) if t < tcut-off = ologn, and explicit upper 
bounds if t > tcut-off- This gives in particular an answer to some questions raised in [SC04, CSC08]. 
Our proofs are mainly inspired by those of the papers [CSST08, Ros91, Por96a, Por96b], in which the 
same kind of results is shown for products of random transpositions in 6„, and for products of random 
rotations in SO(?i). 
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1. Introduction 

1.1. The cut-off phenomenon for random permutations. This paper is concerned with the ana- 
logue for Brownian motions on compact Lie groups and symmetric spaces of the famous cut-off phenom- 
enon observed in random shufhes of cards (c/. [AD86, BD92]). Let us recall this very beautiful result in 
the case of "natural" shuffles of cards, also known as riffle shuffles. Consider a deck of n ordered cards 
1,2, ... ,n, originally in this order. At each time fc > 1, one performs the following procedure: 

(1) One cuts the deck in two parts of sizes m and n ^ m, the integer m being chosen randomly 
according to a binomial law of parameter ^ : 

¥\m^M] = -^f Tr V 
^ ^ 2" \M J 

So for instance, if n = 10 and the deck was initially 123456789X, then one obtains the two blocks 
A = 123456 and B = 789X with probability ^{^°) = |i| - 0.21. 

(2) The first card of the new deck comes from A with probability (cardA)/n, and from B with 
probability (card i3)/n. Then, if A' and B' are the remaining blocks after removal of the first 
card, the second card of the new deck will come from A' with probability (card A')/(n — 1), and 
from B' with probability (card i?')/(n — 1); and similarly for the other cards. So for instance, by 
shuflling A = 123456 and B = 789X, one can obtain with probability l/(g°) ~ 0.0048 the deck 
17283459X6. 

An equivalent description of this shuffle of cards is as follows: one chooses a part P of |1, rt] with uniform 
probability and if M = cardP, then one inserts the M first cards of the deck in P, and the n — M 
last cards of the deck in |l,n] \ P. Denote ©„ the symmetric group of order n, and cr^*''^ the random 
permutation in (3„ obtained after k independent shuffles. One can guess that as k goes to infinity, the 
law P^*^) of (T^''") converges to the uniform law on S„, that is to say that for any fixed permutation a one 
has 



There is a natural distance on the set ^((3„) of probability measures on (3„ that allows to measure 
this convergence: the so-called total variation distance d-Ty. Consider more generally a measurable space 
X with cr-field B{X). The total variation distance is the metric on the set of probability measures ^{X) 
defined by 

dTvi^i,^)=snp{\^iiA)~lyiA)\, A e Bix)} e[0,l]. 

Notice that the convergence in total variation distance is a much stronger notion than the weak conver- 
gence of probability measures: for instance, if (a;„)„gN is a sequence of real numbers converging to some 
a; G M, then the Dirac measures dx^ converge weakly to Sx, whereas dTv(^a:„,^x) is equal to 1 for any 
n (z N. On the other hand, if fi and v are absolutely continuous with respect to a third measure dx on 
X, then their total variation distance can be written as a ^^-norm: 



dTvifJ-,!^) = ^ 



da dv 
ax ax 



dx. 



That said, let us return to the convergence of shuffles. It turns out that with respect to total variation 
distance, the convergence occurs at a specific time fccut-offi that is to say that dTv(IP^'^\ U) stays close to 
1 for fe < fccut-off, and that dxv (]?*■''■' jU) is then extremely close to for k > fccut-off- More precisely, in 
[BD92] (see also [CSST08, Chapter 10]), it is shown that: 

Theorem 1 (Bayer-Diaconis) . Suppose k — ^it^ \ogn + 6. Then, 

dTY{V^^\V) = 1-2(1) (^^\ +0(n-^/^\ , with (t){x) = ^ [ ds. 



OO 
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5*0 for 6 negative, the total variation distance is extremely close to 1, whereas it is extremely close to 
for positive. 

When n = 52, the cut-ofT time is '11°^ ~ 8.55, so "8 shuffles^ are sufficient to obtain an almost uniform 

' 2 log 2 ' 

permutation of the cards of a standard deck". The cut-off phenomenon has been proved for other shuffling 
algorithms {e.g. random transpositions of cards), and more generally for large classes of finite Markov 
chains, see for instance [DSC96, Dia96]. It has also been investigated by Chen and Saloff-Coste for Markov 
processes on continuous spaces, e.g. spheres and Lie groups; see in particular [SC94, SC04, CSC08] and 
the discussion of §1.4. However, in this case, cut-offs are easier to prove for the .if^-^^-norm of pt{x) — 1, 
where ptix) is the density of the process at time t and point x with respect to the equilibrium measure. 
The case of the ^^-norm, which is (up to a factor 2) the total variation distance, is somewhat different. 
In particular, a proof of the cut-off phenomenon for the total variation distance between the Haar measure 
and the marginal law fit of the Brownian motion on a classical compact Lie group (say, a unitary group) 
was apparently not known — see the remark just after [CSC08, Theorem 1.2]. The purpose of this paper is 
precisely to give a proof of this ^^-cut-off for all classical compact Lie groups (special orthogonal groups, 
special unitary groups, and compact symplectic groups), and more generally for all classical symmetric 
spaces of compact type. 

In the two next paragraphs, we describe the spaces in which we will be interested (§1.2), and we precise 
what is meant by "Brownian motion" on a space of this type (c/. §1.3). This will then enable us to explain 
the results of Chen and Saloff-Coste in §1.4, and finally to state in §1.5 which improvements we were able 
to prove. 

1.2. Classical compact Lie groups and symmetric spaces. To begin with, let us fix some notations 
regarding the three classical families of simple compact Lie groups, and their quotients corresponding 
to irreducible simply connected compact symmetric spaces. We follow here most of the conventions of 
[Hcl78, Hel84]. For every n > 1, we denote U(n) — U(n,C) the unitary group of order n; this is the set 
of complex matrices of size n x n such that 

Here, W — is the complex conjugate of the transpose of U. Similarly, the orthogonal group of order 
n is the set 0(n) — 0(n, M) of real matrices of size n x n such that 

CO* = 0*0 = /„; 

and the compact symplectic group of order n is the set USp(n) = USp(n,H) of quaternionic matrices of 
size n X n such that 

SS* = S*S = In. 

Here, S* — denotes the quaternionic conjugate of the transpose of S, the conjugate of a quaternion 
w + ix +iy + kz being w — ix —jy — kz. The orthogonal groups are not connected, so we shall rather work 
with the special orthogonal groups 

SO(n) SO(n, M) = {O e 0(n, M) | dct O = 1} . 

On the other hand, the unitary groups are not simple Lie groups (their center is one-dimensional), so it 
is convenient to introduce the special unitary groups 

SU(n) = SU(n, C) = {[/ G U(n, C) | det C/ = 1} . 

Then, for every n > 1, SU(rt, C), SO(ri,M) and USp(7i,H) are connected simple compact real Lie groups, 
of respective dimensions 

dimR SU(n, C) = - 1 ; dimR SO(n, M) = -^^ — - ; dimR USp(n, H) = 2n^ + n 

^Actually the common say is that seven shuffles are sufficient; indeed the 0(n~^/*) is quite big for n = 52, and one can 
compute dTv{P''^\U) ~ 0.61, dTv(P'''\U) ~ 0.33 and cZtvCP'^^U) ~ 0.17. We shall explain later why the total variation 
distance is decreasing with k; see Lemma 5. 
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and with respective Lie algebras 

su{n) = su(n, C) = {M G M(n,C) | + M = and trM = 0} ; 

50{n) =so(n,M) = {M e M(n,E) | M* + M = 0} ; 

usp(n) = sp(n,H) = {M e M(n,H) | M* + M = 0} . 

The special unitary groups and compact symplectic groups are simply connected; on the other hand, for 
71 > 3, the fundamental group of SO(rt,K) is Z/2Z, and its universal cover is the spin group Spin(n). 

Many computations on these simple compact Lie groups can be performed by using their representation 
theory, which is covered by the highest weight theorem; see Theorem 12 in §2.2. We shall recall all this 
briefly in Section 2, and give in each case the list of all irreducible representations, and the corresponding 
dimensions and characters. It is well known that every simply connected compact simple Lie group is: 

• either one group in the inflnite families SU(n), Spin(n), USp(n); 

• or, an exceptional simple compact Lie group of type Eg, Ey, Eg, F4 or G2. 

We shall refer to the first case as the classical simple compact Lie groups, and as mentioned before, our 
goal is to study Brownian motions on these groups. 

We shall more generally be interested in compact symmetric spaces; let us recall briefly what they are. 
A locally symmetric space is a (connected) Riemannian manifold X such that every point x X has a 
symmetric neighborhood N^, on which is deflned an isometry s : with 

for every geodesic 7 such that 7(0) = x. This is equivalent to the vanishing of the torsion and the 
invariance of the curvature tensor by parallel transport (see [Hel78, Chapter 4]). A symmetric space is 
a locally symmetric space whose local geodesic symmetries can be extended to global isometrics of X. 
In that case, the group of isometrics of X endowed with the compact-open topology is in fact a real Lie 
group, and if G is the connected component of idx in Isom(X), then: 

Theorem 2 (Cartan- Ambrose-Hicks) . The space X is isomorphic (as a Riemannian manifold) to G/K, 
where K is the stabilizer of a point x X and a compact subgroup of G; and (G, K) is a symmetric pair, 
which means that K is included in the group of fixed points G^ of an involutive automorphism of G, and 
contains the connected component {G^)^ of the identity in this group. Moreover, X is compact if and 
only if G is compact. 

This result reduces the classiflcation of symmetric spaces to the classification of real Lie groups and their 
involutive automorphisms. More precisely, consider a simply connected symmetric space X — the simple 
connectedness can be forced by taking the universal cover of the manifold. It is always a direct product of 
Xi X X2 X • • • X of irreducible simply connected symmetric spaces. So, consider an irreducible simply 
connected symmetric space, of compact type. Two cases arise: 

(1) The isometry group G = K x K is the product of a compact simple Lie group with itself, and K 
is embedded into G via the diagonal map k ^ {k,k). The symmetric space X is then the group 
K itself, the quotient map from G to X ~ isT being 

G^ K 
g = (^1,^2) ^ kik^^. 

In particular, the isometrics of K are the multiplication on the left and the right by elements of 
K X K, and this action restricted to iiT C G is the action by conjugacy. 

(2) The isometry group G is a compact simple Lie group, and X is a closed subgroup of it. In this 
case, there exists in fact a non-compact simple Lie group L with maximal compact subgroup K, 
such that G is a compact subgroup of the complexified Lie group L^, and maximal among those 
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containing K. The involutive automorphism 9 extends to L'^ , with K ^ = and the two 
orthogonal symmetric Lie algebras {g,deO) and {l,de9) dual of each other. 

The classification of irreducible simply connected compact symmetric spaces is therefore the following: 
in addition to the compact simple Lie groups themselves, there are the seven following infinite families 

Gr{p + q,q,R) — SO(p + g)/(SO(p) x SO{q)) withp,q > 1 (real Grassmannians); 

Gt{p + q, q, C) — SU(p + q)/S{\J{p) x U(g)) with p,q> 1 (complex Grassmannians); 

Gr{p + q, q, H) = USp(p + g)/(USp(p) x USp(q)) with p,q > 1 (quaternionic Grassmannians); 

SU(n)/SO(n) with n>2 (real structures on C"); 

USp(n)/U(n) with n > 1 (complex structures on H"); 

S0(2n)/U(n) with n>2 (complex structures on M^"); 

SU(2n)/USp(n) with n>2 (quaternionic structures on C^"); 

and quotients involving exceptional Lie groups, e.g. P^(0) = F4/Spin(9); see [Hel78, Chapter 10]. For 
the two last families, one sees \J{n) as a subgroup of S0(2rt) by replacing each complex number x + iy 
by the 2x2 real matrix 



and one sees USp(n) as a subgroup of SU(2n) by replacing each quaternion number ?« + ia; + jy + kz by 
the 2x2 complex matrix 



actually USp(ri, H) is then the intersection of SU(27i, C) and of the complex symplectic group Sp(2ri, C). 
We shall refer to the seven aforementioned families as classical simple compact symmetric spaces (of type 
non-group); again, we aim to study in detail the Brownian motions on these spaces. 



1.3. Laplace-Beltrami operators and Brownian motions on symmetric spaces. In the following, 
we denote driK{k) or dk the Haar measure of a (simple) compact Lie group K, and drjx{x) or dx the Haar 
measure of a compact symmetric space X = G / K , which is the image measure of dr^Q by the projection 
map TT : G — > G / K . In both cases, this is the unique probability measure on the space that is G-invariant. 
We refer to [IIcl84, Chapter 1] for precisions on the integration theory over (compact) Lie groups and 
their homogeneous spaces. 

There are several complementary ways to define a Brownian motion on a compact Lie group X or a 
on compact symmetric space G/K, see in particular [Lia04b]. Hence, one can view them: 

(1) as Markov processes with infinitesimal generator the Laplace-Beltrami differential operator of the 
underlying Riemannian manifold; 

(2) as conjugacy-invariant continuous Levy processes on K, or as projections of such a process on 



(3) at least in the group case, as solutions of stochastic differential equations driven by standard 
(multidimensional) Brownian motions on the Lie algebra. 

The first and the third point of view will be specially useful for our computations. For the sake of 
completeness, let us recall briefiy each point of view — the reader already acquainted with these notions 
can thus go directly to §1.4. 




(1.1) 




(1.2) 



G/K; 
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To begin with, let us precise the Riemannian structures chosen in each case. We recall that the adjoint 
representations of a Lie group G and of its Lie algebra q are defined by 

Ad : G ^ GL(fl) 

g^^{Adg:X^{d,Cg}{X))- (L3) 

ad : g — > End(g) 

X ^ {a&X = {dekd}{X)) . (1.4) 

Then the Lie bracket is [X, 1^] — a.dX{Y). Given a simple compact Lie group K. there is up to a scalar 
constant a unique scalar product on the Lie algebra t = T^K that is Ad(i^)-invariant. A way to construct 
it is to take the mean (Ad fc(-) | Ad fc(-)) dk of an arbitrary scalar product. Alternatively, one can take 
the opposite of the Killing form B{X,Y) = tr(adX o adF), which is negative definite for a simple 
compact group. Then, one can transport this scalar product (• | •)( on any tangent space (• | 

{X I Y)^^^ - {{dkLk-.}iX) I {dkL,-.}{Y)), , (1.5) 

where Lfc-i is the multiplication on the left by k^^. Thus one obtains a Riemannian structure on K 
which is invariant by G ^ K x K, and up to a scalar the unique Riemannian structure with this property. 

When K = SU(7i) or SO(n) or USp(n), the Killing form on { is a scalar multiple of the bilinear form 
X (E)Y 1-^ di{tr(XY)) — the real part is only needed for the quaternionic case. Following the conventions 
of [Lcvll], we shall always consider the following invariant scalar products on {: 

{X\Y)=-^^itriXY)), (1.6) 

with /3 = 1 for special orthogonal groups, j3 — 2 for special unitary groups and unitary groups, and j3 — A 
for compact symplectic groups. Notice that the bilinear form in (1.6) is only proportional to minus the 
Killing form, and not equal to it; for instance, the Killing form of SO(n,K) is 

2j? — 4 

{n~2)tv{XY) = {X\Y), 

n 

and not — | Y). However, the normalization of Formula (1.6) becomes asymptotically the same as 
one-half the Killing form, and on the other hand, it enables one to relate the Brownian motions on 
the compact Lie groups to the "standard" Brownian motions on their Lie algebras, and to the classical 
ensembles of random matrix theory (see the stochastic differential equations at the end of this paragraph). 



Similarly, suppose that X = G/K \& a, simple compact symmetric space (of type non-group). We take 
an Ad(G')-invariant scalar product (• | ■)^ on q (unique up to a scalar), and we denote y the orthogonal 
complement of J in g. Notice that y can be identified with the tangent space oi X — G/K aX eK by the 
map diTe, where Tr{g) — gK. One transports as in Equation (1.5) the restriction to y of the scalar product 
(• I •)g to any tangent space T^X: one thus gets a G-invariant Riemannian structure on X, which is in 
fact unique up to a scalar constant. It is called the Riemannian structure induced by the Riemannian 
structure of G. 

From now on, each classical simple compact symmetric space X — G / K will be endowed with the 
Riemannian structure induced by the Riemannian structure of G = SU(n) or SO(n) or USp(n) that was 
given by Equation (1.6). One may argue that this is not necessarily the best choice of normalization for 
these quotients: in particular, when G = SO(n + 1) and K — SO(n) x S0(1) = SO(n), the Riemannian 
structure defined by the previous conventions on the n-dimensional sphere X = §"(M) differs from the 
"usual" metric" by a factor \/n + 1 . But on the other hand, it will be very convenient for computations 
to have the same normalizations for compact Lie groups and for their symmetric spaces (we did not want 
to define specific conventions for each family). 



'By "usual", we mean the restriction of the standard euclidian metric of K' 
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The Laplace- Beltrami operator on a Riemannian manifold M is the differential operator of degree 2 
defined by 

A/(m)= 5^'(Vx.Vx,/(m)-Vv,^x,/(m)), 

l<i,j<d 

where {Xi, . . . , Xd) is a basis of TmM, {g^-')i,j is the inverse of the metric tensor {gij = {Xi \ Xj)j, jv/)«-J ' 
and denotes the covariant derivative of a vector Y along a vector X and with respect to the Levi- 

Civita connection. In the case of a compact Lie group K, this expression can be greatly simplified as 
follows (see for instance [Lia04b, §2.3]). Fix once and for all an orthonormal basis (Xi,X2, ■ ■ ■ , Xd) of t. 
On another tangent space T^.K, one transports each X^ by setting 

xlik) = {dMix,)eTkK, 

where Rk is the multiplication on the right by k. One thus obtains a vector field X^ = which is 
left-invariant by construction and right-invariant because of the Ad(/ir)-invariance of the scalar product 
on i. Then, 

1=1 ' 

Definition 3. A (standard) Brownian motion on a compact Riemannian manifold M is a continuous 
Feller process {mt)t£R^ whose infinitesimal generator restricted to "^"^{M) is ^ A. 

In the following, on a compact Lie group K ov & compact symmetric space G/ K, we shall also assume 
that Too = e or toq = eK almost surely. We shall then denote jit the marginal law of the process at 
time t, and pf{k) = ^j^ik) or p^{x) = ^^{x) the density of /Ltj with respect to the Haar measure. 
General results about hypoelliptic diffusions on manifolds ensure that these densities exist for i > and 
are continuous in time and space; we shall later give explicit formulas for them (c/. Section 2). 

By using the geometry of the spaces considered and the language of Levy processes, one can give 
another equivalent definition of Brownian motions. The right increments of a random process (5t)tGR+ 
with values in a (compact) Lie group G are the random variables rf = gj^gt, so gt = gs for any times 
s < t. Then, a left Levy process on G is a cadlag random process such that: 

(1) For any times — to < ti < ■ ■ ■ < t„, the right increments rl°,rl^, . . . j?'*""^ are independent. 

(2) For any times s < t, the law of r| only depends on the difference t — s: r^_^ r|. 

Denote Pt the operator on the space "^(G) of continuous functions on G defined by {Ptf){g) = ^[f{99t)]\ 
and /Xj the law of gt assuming that g^ ~ cq almost surely. For G G, we also denote by the operator 
on ^(G) defined by Lhf{g) = f{hg). If {gt)tsM+ is a left Levy process on G starting at g^ — eg, then: 

(1) The family of operators {Pt)tm.+ is a Feller semigroup that is left G-invariant, meaning that 
Pt o Lh = Lh o Pt for all ft. G G and for all time t. Conversely, any such Feller semigroup is the 
group of transitions of a left Levy process which is unique in law. 

(2) The family of laws {^J-t)tl£R-^ is a semigroup of probability measures for the convolution product 
of measures 

(M*^)(/)= / f{gh)dfiig)d,.ih). 

Hence, fis * fJ-t = t^s+t for any s and t. Moreover, this semigroup is continuous, i.e., the limit in 
law lim4_j.o /i* exists and is the Dirac measure 5f.- Conversely, given such a semigroup of measures, 
there is always a corresponding left Levy process, and it is unique in law. 

Thus, left Levy processes are the same as left G-invariant Feller semigroups of operators, and they are 
also the same as continuous semigroups of probability measures on G. In particular, on a compact Lie 
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group, they are characterized by their infinitesimal generator 

defined on a suitable subspace of "^(G). Hunt's theorem (c/. [Hun56]) then characterizes the possible 
infinitesimal generators of (left) Levy processes on a Lie group. In the case of a Levy process with 
continuous paths, it can be stated as follows: 

Theorem 4 (Hunt). The domain of the infinitesimal generator of a continuous left Levy process traced 
on a compact Lie group G contains "^^^(G), and it writes uniquely as 

where the diffusion part {a'^^)i<i.j<d is a non-negative symmetric matrix, and the drift part ic^)i<i<d is 
any vector. 



Suppose then that {gt 



tGR+ 



is a continuous Levy process on a simple compact Lie group G, starting 



from e and with the additional property that {hgth 



and {g, 



have the same law in '^(M+, G) 



for every h. These hypotheses imply in Equation (L8) that {a^'-')i.j = iaSij)ij for some constant a > 0, 
and that (c*)^ = 0. Thus, on a simple compact Lie group K, up to a linear change of time t i— ^ at, 
a conjugacy-invariant continuous left Levy process is a Brownian motion in the sense of Definition 3. 
Similarly, on a simple compact symmetric space G/K, up to a linear change of time, the image {gtK)t^TR^ 
of a conjugacy-invariant continuous left Levy process on G is a Brownian motion in the sense of Definition 
3. This second definition of Brownian motions on compact symmetric spaces has the following important 
consequence: 

Lemma 5. Let fit be the law of a Brownian motion on a compact Lie group K or on a compact symmetric 
space G/K. The total variation distance dTv(A'tj H^ar) is a non-increasing function oft. 

Proof. First, let us treat the case of compact Lie groups. We define the convolution product of two 
functions on K by 

(/i*/2)(A;)= / fi{h)f2{h-^k)dh. 



K 



If /i,/2 are in Ji'^{K,dr]K), then /i * /2 is again in .^^{K), with * /2||j^i(if) < ||/2||^i(k) 



Now, since fig+t — f^s*fJ't, the densities of the Brownian motion also satisfy pf^^ 



*Pt^ ■ Consequently, 



< 



II 



\iP. 



K 



l)*pfl 
K 



WP's - lIljSfH^) WVt \\i£^{K) = WPs - l||j^i(G) =2dTvifJ-s,VK)- 

The proof is thus done in the group case. For a compact symmetric space X — G/K, denote pf the 
density of the Brownian motion on G, and the density of the Brownian motion on X. Since the 
Brownian motion on X is the image of the Brownian motion on G by tt : G — > G/if, one has: 



yx = gK, pf{x)= I p'r{gk)dk. 

Ik 



As a consequence. 



Ibf+t - MIj^Hx) = / \pf+ti9K) ~l\dg 



K 



{p%t{gk)-l)dk 



dg 



{pf{h-'gk)-l)pf{h) dkdh 



KxG 



dg 



a 



{pfih-'gK)-l)pf{h)dh 



G 



dg 



< / \Psih-'9K) - 1)1 \pfih)\ dhdg = llpf - \\p?\\j^HG) ^ \\Ps ~ M\j^HX), 

'GxG 
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SO dTv{^s+t,Vx) < rfTv(Msj??x) also in the case of symmetric spaces. Notice that our proof works 
for any projection on a homogeneous space G/K of a Levy process (not necessarily continuous) traced 
on a compact Lie group G. Later, this property will allow us to compute estimates of dTy{iit,r]x) 
only for t around the cut-off time. Indeed, if one has for instance an (exponentially small) estimate of 
1 — dTv{fj,to,r]x) at time = {1 — e) icut-off , then the same estimate will also hold for 1 — dT\j{^t,^x) 
with t <to. □ 

Remark. Actually, the same result holds for the ^P-noim oi pt{x) — l, and in the broader setting of Markov 
processes with a stationary measure; see e.g. [CSC08, Proposition 3.1]. Our proof is more elementary 
since it only uses the triangular inequality | f{x)dx\ < |/(a^)| dx and algebraic manipulations with 
the semigroup of measures. 

A third equivalent definition of Brownian motions on compact Lie groups is by mean of stochastic 
differential equations. More precisely, given a Brownian motion {kt)tetL+ traced on a compact Lie group 
K, there exists a (trajectorially unique) standard d-dimensional Brownian motion {Wt)teR+ on the Lie 
algebra £ = Yect{Xi, X2, ■ . . , Xd) such that, for every function / € ^^(iir), the following stochastic 



dx 



2 



■ ds. 



differential equation holds: 

^ /' ^ „ .... . ;,,.) + E /' ^ ■ 4 S /' 

i=l "'''=0 ^■'^i Js=0 ^■'^i ^ J s. 

where the symbol o corresponds to Stratonovitch's stochastic integral, and the symbol • corresponds to 
Ito's stochastic integral. So for instance, on a unitary group U(n, C), the Brownian motion is the solution 
of the SDE 

UQ = In ; dUt=iUfdHt-^Utdt, 

where (Ht)teR+ is a Brownian hermitian matrix normalized so that at time t = I the diagonal coeffi- 
cients are independent real gaussian variables of variance 1/n, and the upper-diagonal coefficients are 
independent complex gaussian variables with real and imaginary parts independent and of variance l/2n. 
Similarly, on SO(n,M), the Brownian motion is the solution of the SDE 

TL — 1 

Oo^In ; dOt = Ot ■ dAt - Of dt, 

2n 

where {At)ti£R^ is an skew-symmetric Brownian matrix normalized so that at time t = 1 the upper- 
diagonal coefficients are independent real gaussian variables of variance 1/n. In the general case, let us 
introduce the Casimir operator 

d 

C = ^X,(»X,. (1.9) 

1=1 

This tensor should be considered as an element of the universal enveloping algebra U{V}. Then, for every 
representation tt : A' — > GL(V^), the image of C by the infinitesimal representation dir : t/(4) — ^ End(V^) 
commutes with (i7r(g). In particular, for an irreducible representation V , dir^G) is a scalar multiple Kyidy 
of idy. Suppose K is a classical simple Lie group. Then its "geometric" representation is irreducible, 
so X]f=i -^i = cig /„ if one sees the Xi's as matrices in M(n, M) or M(n, C) or M(n, H). The stochastic 
differential equation satisfied by a Brownian motion on K is then 

ko ^ eK ; dkt = kt dBt -I- -^h dt, 

where Bt — X^iLi -^i is a standard Brownian motion on the Lie algebra t. The constant is given 
in the classical cases by 

_ n'^ - I _ 2n + l 

\ J n 2n 

see [Levll, Lemma 1.2]. These Casimir operators will play a prominent role in the computation of the 
densities of these Brownian motions (c/. §2.2), and also at the end of this paper (§4.1), see Lemma 27. 
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1.4. Chen-Saloff-Coste results on .jSf'-cut-offs of Markov processes. Fix p G [1, oo), and consider 
a Markov process X = {xt)t£Kj^ with values in a measurable space {X, B{X)), and admitting an invariant 
probability ?/. One denotes fXt^x the marginal law of Xt assuming xq = x almost surely, and 



(X) = max 



X 



d4M^ 
drj 



with by convention 



(X) = 



ifp=l, 
XD if p > 1, 

when fit J. is not absolutely continuous with respect to 77. This is obviously a generalization of the total 
variation distance to the stationary measure. In virtue of the remark stated just after Lemma 5,ti-^ ^^{X) 
is always non-increasing. A sequence of Markov processes (X^"-')„gN with values in measurable spaces 
S(X*^"'))„gN is said to have a max-^P -cut-off with cut-off times (t("^)„gN if 



lim sup (X(")) = 



lim inf (X(")) = lim sup ^(X^")) M > 

n-J-oo \^t<(i_e)t(") 



for every e E (0, 1) — usually M will be equal to 2 or +00. A generalization of Theorem 1 ensures that 



these ^P-cut-offs occur for instance in the case of rifhe shuffles of cards, with i'") 
p e [l,+oo). 



3 log 71 



for every 



In [CSC08], Chen and Saloff-Coste shown that a general criterion due to Peres ensures a .if ^'^ ^-cut-off 
for a sequence of Markov processes; but then one does not know necessarily the value of the cut-off time 
f^"^ Call spectral gap A(X) of a Markov process X the largest c > such that for all / G ^^{X,ri) and 
all time t, \\{Pt — 77)/||^2(x) < e^*^ where {Pt)tm+ stands for the semigroup associated to the 

Markov process. For simplicity, we shall assume that these operators are normal; this in particular 
the case for Brownian motions on compact Lie groups or symmetric spaces. 



Theorem 6 (Chen-Saloff-Coste). Fix p e (l,cx)). One considers a family of normal Markov processes 
(X'^"))„gN with spectral gaps A*^"-*, and one assumes that limt^oorff(X(")) = for every n. For eo > 
fixed, set 

= inf{i : (X(")) < eo}. 
The family of Markov processes has a max- -cut-off if and only if Peres ' criterion is satisfied: 

lim A(")t(") = +C50. 

n— >C30 

In this case, the sequence (i*-"-')ngN gives the values of the cut-off times. 

Notice that a lower bound on t^"^ also ensures the cut-off phenomenon; but in this case, the cut-off 
time remains unknown. Nevertheless, an important application of this general criterion is the following 
result (see [CSC08, Theorem 1.2], and also [SC04, Theorem 1.1 and 1.2]): 

Corollary 7 (Saloff-Coste). Consider the Brownian motions traced on the special orthogonal groups 
SO(n,M), with the normalization of the metric detailed in the previous paragraph. They exhibit for every 
p £ (1,00) a cut-off with t*^"-* asymptotically between 21ogn and 41ogn — notice that t*^"^ depends on p. 

Indeed, the spectral gap stays bounded and has a non-negative limit (which we shall compute later), 
whereas t^"^ was shown by Saloff-Coste to be a O(logn). Similar results are presented in [SC04] in 
the broader setting of simple compact Lie groups or compact symmetric spaces, but without a proof of 
the cut-off phenomenon (Saloff-Coste gave a window for i^") for every p € [l,-|-cx)]). The main result 
of our paper is that a cut-off indeed occurs for every p £ [l,-fcx)(, for every classical simple compact 
Lie group or classical simple compact symmetric space, and with a cut-off time equal to logn or 21ogn 
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depending on the type of the space considered. In particular, the main improvements in comparison to 
the aforementioned theorems are the foUowings: 

(1) the case p = 1 is now included; 

(2) one knows the precise value of the cut-off time. 

To obtain this result, we more or less had to compute t^'^\ Hence this paper is mostly a case-by-case 
analysis, with precise explicit estimates and a few computational tricks (in particular we depart a lot 
from the generality of Theorem 6) . 

1.5. Statement of the main results and discriminating events. The precise statement of our main 
result is the following: 

Theorem 8. Let iit be the marginal law of the Brownian motion traced on a classical simple compact 
Lie group, or on classical simple compact symmetric space. There exists positive constants a, 7^, 7a, c, 
C and an integer uq such that in each family, for all n > hq, 

Ve e (0,1/4), dTv(Mt,Haar) > 1 — if t ^ a (1 - e) logn; (1.10) 

c 

V£e(0, 00), dTv(Mt.Haar) < if t ^ a {I + e) logn. (1.11) 

n ^ 

The constants a, 7f, and 7^ are determined by the type of the space considered, and then one can make 
the following choices for tiq, c and C : 



K or G/K 




a 


7b 


la 


"0 


c 


C 


SO(n,M) 


1 


2 


2 


2 


10 


36 


6 


SU(n,C) 


2 


2 


2 


4 


2 


8 


10 


USp(n,H) 


4 


2 


2 


2 


3 


5 


3 


Gr(n,g,M) 


1 


1 


1 


1 


10 


32 


2 


Gr(ri, q, C) 


2 


1 


1 


2 


2 


32 


2 


Gr(n,(7,H) 


4 


1 


1 


1 


3 


16 


2 


S0(2n,M)/U(n,C) 


1 


1 


2 


1 


10 


8 


2 


SU(n,C)/SO(n,M) 


2 


1 


2 


2 


2 


24 


8 


SU(2n,C)/USp(n,H) 


2 


1 


2 


2 


2 


22 


8 


USp(n,H)/U(n,C) 


4 


1 


2 


1 


3 


17 


2 



As the function t 1— )■ cZtvIa**! Haar) is non-increasing in t, the aspect of this function in the scale 
t oc logn is then always as on Figure 1. 

dTv(Mt:Haar) 




a log n 



Figure 1. Aspect of the function t 1— > (ixv (/^t 7 Haar) for the Brownian motion on a 
classical simple compact Lie group or on a classical simple compact symmetric space. 
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The constants c and C in Theorem 8 can be sHghtly improved by raising the integer np; the restriction 
n > riQ will only be used to ease certain computations and to get reasonable constants c and C. A 
result similar to Theorem 8 has been proved by Rosenthal and Porod in [Ros94, Por96a, Por96b] for 
random products of (real, or complex, or quaternionic) reflections. Our proofs are really inspired by their 
proofs, though quite different in the details of the computations. For the upper bound (1.11), it has long 
been known that if A(X„) denotes the spectral gap of the heat semigroup associated to the infinitesimal 
generator L — |A, then for n fixed, the total variation distance d'YY{^J.t,'nXn) decreases exponentially 
fast: 

We refer for instance to [Lia04a], where such a bound is proved for any "non-degenerate" Levy process 
on a compact Lie group. Consider now the family of spaces (X„)„gN, and suppose that C(X„) — C , 
and that A(X„) stays almost constant to A — this last condition is ensured by the normalization (1.6). 
Then, one obtains for t = {1 + e) j \ogn the bound 

C 

Thus in theory, the upper bound (1.11) follows from the calculations of the constants C(X„) and A(X„) 
in each classical family. It is very hard to find directly a constant C(X„) that works for every time t. 
But on the other side, by using the representation theory of the classical simple compact Lie groups (c/. 
Section 2), one can determine series of negative exponentials that dominates the total variation distance; 
see Proposition 16. In these series, the "least negative" exponentials give the correct order of decay 
\{Xn). It remains then to prove that the other terms can be uniformly bounded. This is tedious, but 
doable, and these precise estimates are shown in Section 3: we shall adapt and improve the arguments 
of [Ros94, Por96a, Por96b, CSST08]. 

As for the lower bound (1.10), it is obtained by looking at discriminating events, that have a probability 
close to 1 with respect to a marginal law /it with t < tcut-off, and close to with respect to the Haar 
measure. In the case of riffle shuffles, the sizes of the runs of a permutation enable one to discriminate 
a random shuffle of order k < fccut-off from a uniform permutation. Indeed, call run of a permutation 
w = W1W2 ■ ■ -Wn a maximal'^ subsequence Wi^Wi^ . . .Wi^ with ii < 12 < ■ ■ ■ < i,- and Wi.+i = Wi. + 1 
for every j. For example, the permutation 17283459X6 has two runs, namely, 123456 and 789X. Now, 
consider for instance a random permutation cr of size n = 52 obtained by shuffling 3 times the deck. Each 
shuffle splits a run into at most 2 runs, so a has at most 8 = 2"^ runs, and the probability that one of them 
is of length 1 is extremely small. On the contrary, a random uniform permutation has with probability 
close to 1 at least one run of length 1. Indeed, 

(7 has a run of length 1 <^=^ 3i e \2,n - l\ , a^^{i - I) > a^^{i) > a^^{i + I). 

For a given i E |2,n — 1], the probability under U of the event Ai = {a^^{i — 1) > tT^^(i) > cr^^{i + 1)} 
is |. Moreover, if |i — j | > 3, then Ai and Aj are independent under U. As a consequence. 



U [a has a run of length 1] > 1 — U 



n 

»e[2,n-ll 
j=3fc+2 



= 1 - 



5xLtJ 
6 



Thus, the event A — {a has a run of length 1} separates the uniform measure U from the measure 
and by improving this argument, one can show similar inequalities for P^'^^ with k < fccut-off- 

Now, in the case of a Brownian motion on a classical compact Lie group, this is the trace of the 
matrices that allows to discriminate Haar distributed elements and random Brownian elements before 
cut-off time. Indeed, consider for instance a random unitary matrix [/„ of size n, taken under the Haar 



■^By "maximal", we mean that one cannot find a larger subsequence in the sense of the inclusion; however, runs may have 
distinct lengths. 



14 



PIERRE-LOiC MELIOT 



measure or under the marginal law fit of the Brownian motion at a given time t. Then, tr J7„ is a complex 
valued random variable, and we shall see that 

E[|trC/„-m|2] < 1, 

where m is the mean of trUn', and this, for any n > 1 and any time t > if Un ^ fJ-t- However, m ~ 
under the Haar measure, whereas |m| ^ 1 for t < ^cut-off- So, the trace of a Brownian unitary matrix 
before cut-off time will never "look the same" as the trace of an Haar distributed unitary matrix. Up to 




Figure 2. Density of the trace tr[/„ of a random unitary matrix, with C/„ ~ Haar for 
the left peak, and C/„ ~ Mt<tcut-off ^'^^ the right peak. 



a minor modification, the same argument will work for special orthogonal groups and compact special 
orthogonal groups — in this later case, the trace of a quaternionic matrix of size n is defined as the trace 
of the corresponding complex matrix of size 2n, cf. the remark at the end of §1.2. 

Over the classical simple compact symmetric spaces, the trace of matrices will be replaced by a zonal 
spherical function "of minimal non-zero weight"; these minimal zonal spherical functions are also those 
that give the order of decay of the series of negative exponentials that dominate dTvif^t, Haar) after the 
cut-off time. This argument for the lower bound was already known, since it has been used successfully 
in [SC94] to prove the cut-off phenomenon over spheres: we have simply extended it to the case of general 
simple compact symmetric spaces (cf. Section 4). 

An important consequence of Theorem 8 and its proof is that one also has a max-^^-cut-off for every 
p G [1, oo]. Moreover, the value of the cut-off time is known when p £ [1, 2]. More precisely: 

Corollary 9. For every p G [l,+oo], the family of Brownian motions (X„)„£n traced on simple compact 
Lie groups {Kn)nefi one of the three classical families (respectively, on simple compact symmetric spaces 
of type non-group (X„)„gN in one of the seven classical families) has a max-^^ -cut-off . If p e [1,2], it 
is with respect to the sequence i*^"^ = 21ogn (respectively, t'") = \ogn). 

Proof. The upper bound in Theorem 8 will be shown by using Cauchy-Schwarz inequality and estimating 



the ^■^-norm of 



1 



which can be written as a series Sn{t) of negative exponentials. Section 3 is 



dr)>. 

devoted to the proof of the fact that Sn{t) is small after cut-off time, and on the other hand, the same 
series trivially goes to infinity before cut-off time, because some of its terms go to infinity (consider for 
instance the term indexed by the "minimal" label identified in Lemma 17). Thus, our proof of Theorem 
8 implies readily a ^^-cut-off; and since the Brownian motion is invariant by action of the isometry 
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group, it is even a max-^^-cut-ofF. We can then use [CSC08, Theorem 5.3] to obtain the existence of a 
max-^P-cut-ofF for every p g (1, +00], and the comparison theorem of mixing times [CSC08, Proposition 
5.1] to get the value of the cut-off time when p is between 1 and 2. When p = +00, [CSC08, Theorem 
5.3] also gives the value of the cut-off time: it is 41ogn in the group case, and 21ogn in the non-group 
case. However, when p g (2, +00), one still does not know the value of the mixing time: one has only the 
window alogn < f*^"^ < 2Q!logn. □ 



Organization of the paper. In Section 2, we recall the basics of representation theory and harmonic 
analysis on compact symmetric spaces, with a particular emphasis on explicit formulas since we will need 
them in each case. All of it is really classical and of course well-known by the experts. Nevertheless, in 
order to be as self-contained as possible, we have decided to include in our paper some details of the proof 
of the totally explicit formulas appearing in Theorem 15. It also enables us to fix the notations related 
to the harmonic analysis of the classical compact Lie groups and compact symmetric spaces. In Section 
3, we use the explicit expansion of the densities to establish precise upper bounds on \\pt — '^\\^''(x,ri)] 
by Cauchy-Schwarz we obtain similar upper bounds on dTylptTV)- The main idea is to control the 
growth of the dimension of an irreducible spherical representation involved in the expansion of pt when 
the corresponding highest weight grows in the lattice of weights (§3.2). The crucial fact, which was 
apparently unknown, is that precisely at cut-off time, the quantity 

{D^y g-tcut-off-B„(A) jj-^ |-]^g group case, 
£)A g-tcut-oftBn(A) jj-^ |-]^g non-group case, 

stays bounded for every n and every A; Z?^ being the dimension of the irreducible or spherical irreducible 
representation of label A, and — i?„(A) the associated eigenvalue of the Laplace-Beltrami operator. Com- 
bining this argument with a simple analysis of the generating series 

A partition z>l 

this is sufficient to get a correct upper bound after cut-off time. 

Section 4 is then devoted to the proof of the lower bounds. We use in each case a "minimal" zonal 
spherical function (the trace of matrices in the case of groups; see §4.1), and we compute its expectation 
and variance under Haar measure and Brownian measures (§4.2). A simple application of Bienayme- 
Chebyshev's inequality will then show that the chosen zonal spherical function is indeed discriminating. 
An algebraic difficulty occurs in the case of symmetric spaces G/K of type non-group, as one has to 
compute the expansion in zonal functions of the square of the discriminating zonal function, and this 
is far less obvious than in the case of irreducible characters. The problem is solved by writing the 
discriminating zonal function in terms of the coefficients of the matrices in the isometry group G, and 
by computing the joint moments of these coefficients under a Brownian measure. The combinations of 
negative exponentials appearing in these formulas are then in correspondence with the expansions of the 
squares of the discriminating zonal spherical functions. 

Acknowledgements. Many thanks are due to Yacine Barhoumi, Philippe Biane, Florent Benaych- 
Georges, Reda Chhaibi, Djalil Chafai', Kenneth Maples, Ashkan Nikeghbali and Simon Pepin-Lehalleur 
for discussions around the cut-off phenomenon and the representation theory of Lie groups. 
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2. Fourier expansion of the densities 

In this section, we explain how to compute the density pf{k) or p^{x) of the marginal law of the 
Brownian motion traced on a compact Lie group if or a compact symmetric space X = G/K. This 
computation is done in an abstract setting for instance in [Lia04a] or [Appll], and we shall give at the 
end of this section its concrete counterpart in each classical case, see Theorem 15. The main ingredients 
of the computation are: 

(1) Peter- Weyl's theorem and its refinement due to Cartan, that ensures that the matrix coefficients 
of the irreducible representations of K (respectively, of the irreducible spherical representations of 
G) form an orthogonal basis of ^^{K. -q) (respectively, of ^^{G/K^ rj)). and that the irreducible 
characters (respectively, the zonal spherical functions) form an orthogonal basis of ^^{K,ri)^ 
(respectively, of ^^{G / K,t])^); see §2.1. 

(2) the classical highest weight theory, that describes the irreducible representations of a compact 
simple Lie group and give formulas for their dimensions and characters; see §2.2. 

On these subjects, we refer to the two books by Helgason [Hcl78, Hcl84], and also to [BD85, Var89, FH91, 
Far08, GW09] for the representation theory of compact Lie groups. We shall only recall what is needed 
in order to have a good understanding of the formulas of Theorem 15. 



2.1. Peter- Weyl's theorem and Cartan's refinement. Peter-Weyl's theorem deals with the exten- 
sion to arbitrary (possibly non-commutative) compact topological groups of the well-known Fourier theory 
for functions on the circle T. Namely, if / G ^^(T), then of course 

/w = E/w«™' ; ii/iii^^(T) = E|/(")r 

where f{n) — J^^ /(e~'^) e"'^|| is the n-th Fourier coefficient of /. 

Let i^T be a compact (Lie) group, and K be the set of isomorphism classes of irreducible complex linear 
representations of K. Each class A € if is finite-dimensional, and we shall denote V'^ the corresponding 
complex vector space; : K ^ UCV^) the representation morphism"'; D'*' = dime the dimension of 
the representation; x'*'(') = tr /?'''(•) the character; and x'^(-) = x^{')/D^ the normalized character. An 
Hermitian scalar product on End(l/'*') is (M | A^)End(y^) = iT::{M^ N). For every class A and every 
function / G ^"^{K), we set 

7(A) = / f{k)p\k)dk- 
Jk 

this is an element of End(y^). We refer to [BD85, Far08] for a proof of the following results. 
Theorem 10 (Peter- Weyl). The (non- commutative) Fourier transform 

T ■ ^^{K,ri) End(V^^) 
xeR 

XeR 

realizes an isometry and an isomorphism of (non-unital) algebras, where ^'^(K) is endowed with its usual 
scalar product and the convolution product, and End{V^), the Hilbertian completion of the algebraic 



One can always endow with an Hermitian scalar product such that K acts by unitary transformations. 
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direct sum, is endowed with the Hilhertian sum of the scalar products {■ \ •)End(y^)- ''■f f ^ -^"^{K), 
then 



J{k) = tr (/(A) p^k)) = Y tr (^^ f{h-'k) p\h) dh^ 



(2.1) 



XeK XeK 

2 



l/II^^W - E \\fW - E D'tr (/(A)t/(A)) (2.2) 



XeK XeK 



where /(A) = /-(A) = /(fc-i) p^{k) dk. 



X\f)^ l^fik)x\k)dk^tT{fi\)) ; x^(/) 



Suppose now that / is in ^^{K,ri)^ , the subalgebra of conjugacy-invariant functions such that 
f{hkh^^) = f{k) for every h,k E K. Then for every class A, /(A) = x^(f) idyx, where 

xHf) 

So in this case, the Fourier expansion (2.1) and the Parseval identity (2.2) become 

/(fc) = E(^'f x'(r)x'(fc) ; = E ix'mi'' 

XeK XeK 
and in particular, the irreducible characters form an orthonormal basis of ^^{K)^ . 

Cartan gave a statement similar to Theorem 2.1, but for ^^{G / K,rj) where X — G/K is. a, simply 
connected irreducible compact symmetric space. Call spherical an irreducible representation {V^^p^) 
of G such that {V^)^ , the space of vectors invariant by p^(K), is non-zero. Then, it is in fact one- 
dimensional, so one can find a vector of norm ||e'^|p — 1, unique up to multiplication by z G T, such 
that {V^)^ — Ce^. Denote then "^^{G/K) the set of functions from G to C that can be written as 

/(5) = /.(.9) = («|p^5)(e^)) with«eF\ (2.3) 
Such a function is right- X-invariant, so it can be considered as a function from G/K to C. 

Theorem 11 (Cartan). Let G^ be the set of spherical irreducible representations of G. The Hilbert 
space ^^{G / K,r]) is equal to the orthogonal sum ©_;^gg/f '^^{G/K). This decomposition corresponds to 
the Fourier expansion 



XeC- 

for f eJ^HG/K). 



f{gK)= J2 D^trff f{h-^gK), 



p^{h)dh\ (2.4) 



One recovers Equation (2.1) from (2.4) by looking at symmetric spaces of type group [K x K)/K: 
spherical representations in G^ are then the ® = End(V^'^) with XeK, and the corresponding 
spherical vector is up to a scalar idyx . 

In each space "^"^(G/if), the space of left ii'-invariant functions is one-dimensional, and it is generated 
by the zonal spherical function 

cb\gK)^{e'\p\g){e')). 

These spherical functions form an orthogonal basis of ^^{X)^ when A runs over G^ . So, a iiT-invariant 
function writes as 

f{gK)^ E D^<f'\n^^i9K), 
xeG'< 

where (f>^{f) = j^^^^ f{x) 4>^{x) dx = {e^ | f{gK) p^{g){e^) dg) . To conclude with, notice that the 
decomposition given by Theorem 11 is the decomposition of J^^{G/K,vi) in common eigenspaces of the 
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elements of ^{G/K), the commutative algebra of differential operators on X that are G-invariant. Thus, 
there are morphisms of algebras : 2i{G/K) C for every A G such that 

T{f^) = c\T) 

for every T e &{G/K) and every e '^^{G/K). 

2.2. Highest weight theorem and Weyl's character formula. The theory of highest weights of 
representations enables us to identify K or G^ , and to compute the coefficients c^{L) associated to the 
Laplace-Beltrami operator. If G is a connected compact Lie group, its maximal tori are all conjugated, 
and every element of K is contained in a maximal torus T. Denote W = Norm(T)/T the Weyl group 
of G, and call weight of a representation ^ of G a group morphism a; : T — )■ T such that = {u G 
y I Vi € T, t ■ V — uj{t) ■ v} 0. Every representation y of G is the direct sum of its weight subspaces 
V^, and this decomposition is always M^-invariant. Besides, the set of all weights of all representations 
of G is a lattice Zfi whose rank is also the dimension of T. In the following, we fix a VF-invariant scalar 
product'' on the real vector space Rfl = Zfl ®z IR, and a closed fundamental set G for the action of the 
Weyl group on Mfi. We call dominant a weight uj that falls in the Weyl chamber G. 

A root of G is a non-zero weight of the adjoint representation (1.3). The set of roots $ is a root 
system, which means that certain combinatorial relations are satisfied between its elements: for instance, 
if a £ <&, then — q G $, and if a, /3 S $, then the set of fc e Z such that a + € <i> is an interval. There 
is a unique way to split $ in a set of positive roots and a set $_ = — $+ such that 

G = e I Va e $+, {x\a)> 0}. 

Call simple a positive root a that cannot be written as the sum of two positive roots; and simple coroot 
an element a = j^f^ with a simple root. Then, a distinguished basis of the lattice Zi7 is given by the 
fundamental weights voi^W2, ■ ■ ■ the dual basis of the basis of coroots. Hence, the sets of weights 
and of dominant weights have the following equivalent descriptions: 

Zf7 = Zwi jx e Ml] I Va e $, 
i=l ^ 

Dom(Zn) = Nti7j = |a; e ]Rf7 I Va e $, 

Example. The weight lattice of G = SU(3, C) is drawn on Figure 3; T — {{x^y^z) e T"^ | xyz — 1}; 
a-i{t) = xy^^, a2{t) = yz^'^; Wi{t) = x, W2{t) = xy. 



(a I a) 

{x \ a) 
(a I a) 



e z 




roots 

simple roots 
weights 

fundamental weights 
Weyl chamber 



• • • • 



Figure 3. The weight lattice of type A2, corresponding to the compact group SU(3, C). 



E.g., the dual of the scalar product given by Equation (1.6), where Rf! is identified with i* by mean of i-> deui for 
LJ £ ZQ. Actually, one has to choose precisely this scalar product so that Equation (2.7) has the correct normalization. 
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In any case, the Weyl group is a Coxeter group generated by the symmetries along the simple roots 
ai, q;2, • ■ • , ar- In particular, it admits a signature morphism e : W {ill- 
Theorem 12 (Weyl). Suppose that G is a semi-simple simply connected compact Lie group, and consider 
the partial order induced by the convex set C on Mf2. Every irreducible representation V of G has a unique 
highest weight Wq for this order, which is then of multiplicity one and determines the isomorphism class of 
V. Thus, the dominant weights parametrize the irreducible representations ofG. Moreover, the restriction 
to T of the irreducible character associated to a dominant weight A is given by 

where p is the half-sum of all positive roots, or equivalently the sum of the fundamental weights. This 
formula degenerates into the dimension formula 

Z?^=dimy^ = n^f^±iit^. (2.6) 



Notice that Equation (2.5) gives the values of the characters on any element g G G, since they are 
all conjugated to an element t G T, unique up to action of the Weyl group. Thus, Theorem 12 makes 
Equation (2.1) essentially explicit in the case of a conjugacy invariant function on a (semi-)simple compact 
Lie group K; in particular, we shall see in a moment that the highest weights are labelled by partitions 
or similar combinatorial objects in all the classical cases. 

The case of a compact symmetric space X — G/K of type non-group is more involved as it requires 
the knowledge of the involutive automorphism underlying the symmetric pair. Denote 9 an involutive 
automorphism of a semi-simple simply connected compact Lie group G, with K ~ G^ (see Theorem 2). 
Set P ~ {g G G \ 9{g) = g~^}; one has then the Cartan decomposition G ~ KP. In addition to the 
assumptions of Theorem 12, one has now to suppose that the maximal torus T C G is chosen so that 
T^ — T and P H T is a maximal torus in P. This will be specially important for the harmonic analysis 
on Grassmannian varieties, see the next paragraph. 

Theorem 13 (Cartan-Helgason). Under the previous assumptions, the spherical representations in G^ 
are precisely the irreducible representations in G that are trivial on K C\T — T^ . 

This subgroup of T ~ T*" is always the product of a subtorus T''-'" with an elementary abelian 2-group 
(Z/2Z)*; as we shall see in §2.3, this corresponds to additional conditions on the size and the parity of 
the parts of the partitions labeling the highest weights in G^ (in comparison to G) . The corresponding 
zonal spherical functions (j)^ do not have in general an expression as simple as (2.5); see however [HS94, 
Part 1]. For most of our computations, this will not be a big problem, since we shall only use certain 
properties of the spherical functions — e.g., their orthogonality and the formula for the dimension 
— and not their explicit form; see however §4.1. 

The last ingredient in the computation of the densities is the value of the coefficient c'^(i) such that 

for every function f^ either in M^{K) ~ Vcct({fc H' {p^{k))ij, 1 < i,j < D^}) in the group case, or 
in '^^{G/K) in the case of a symmetric space. In the group case, by comparing the definition of the 
Casimir operator (1.9) with the definition of the Laplace-Beltrami operator (1.7), one sees that c^{L) is 
also ^ , where kx is the constant by which the Casimir operator G acts via the infinitesimal representation 
dp^ : U{i) End{V^) — cf. the remark at the end of §1.3. This constant is equal to 

= -(A + 2p| A), (2.7) 



20 



PIERRE-LOiC MELIOT 



see [Appll, Equation (3.4)] and the references therein, or [Levll] and [Far08, Chapterl2] for a case-by- 
case computation. More precisely, denoting Eij the elementary matrices in M(n, k) with fc = M, C or H, 
the Casimir elements with respect to the scalar product (1.6) write as 

(g)2 



C'su(ji) 



l<i<j<n 

1 " 

- ^ lE,, ® lEu + ]E,^ ® ]E,i + VEu ® "iiEi, 



usp(ji) 



2n 



^ / E,j - Ej, ^ f + '^Eji \'^^ ^ f i^ij + j-Eji \ "^^ I ^Ejj + kEj 

see [Lcvll, Lemma 1.2] — beware that the tensor product is over M, since we deal with real Lie algebras. 
One can then use the relations of the universal Lie algebra U{i^) to rewrite C{ as the sum T + N of 

• an element T in U{t) that acts on by the coefficient of Formula (2.7); 

• an element N that vanishes on the highest weight vector of y , because part of the left ideal of 
[/(t*') generated by the elements of n, n being the nilpotent subalgebra of i"^ associated to the 
choice of the Weyl chamber; 

see the argument of [Far08, p. 287] for more details. In the case of a compact symmetric space, the same 
Formula (2.7) gives the action of ISP^^ on '^^{G / K). Indeed, remember that the Riemannian structures 
on G and G/K are chosen in such a way that for any / S '^°"{G) that is right if- in variant, 

A«/^(/)(gX) ^ AG(/)(.g). 

Consider then a function in '^^{G/K), viewed as a function on G. In Definition (2.3), / appears clearly 
as a linear combination of matrix coefficients of the spherical representation A G G^', so the previous 
discussion holds. 



2.3. Densities of a Brownian motion with values in a compact symmetric space. Let us now 

see how the previous results can be used to compute the density (fc) or (x) of a Brownian motion 
on a compact Lie group or symmetric space. These densities are in both cases X-invariant, so they can 
be written as 

pf{k)=Y,^a,it)x\k) or pfix)= J2 «a(O0'(^) 

by using either Peter- Weyl's theorem in the case of conjugacy-invariant functions on K, or Cartan's 
theorem in the case of left if -invariant functions on G/K. We then apply y = ^ | j_g to these formulas: 



Apf (fc) >^ KA .,^^A.,^ dpfjt) ^ daxjt) 



and similarly in the case of a compact symmetric space. Thus, ^^^^^^ — ^ax{t) and a\{t) = a\{0)e 
for every class A. The coefficient a\{Q) is given in the group case by 



J K 

and in the case of a compact symmetric space of type non-group by 



G 
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Combining this with the highest weight theory exposed in the previous paragraph, one obtains 

Proposition 14. The density of the law fit of the Brownian motion traced on a classical simple compact 
Lie group K is 

and the density of the Brownian motion traced on a classical simple compact symmetric space G/ K is 



Let us now apply this proposition in each classical case (see [BD85], [FH91, Chapter 24] and [Hcl78, 
Chapter 10] for precisions on these computations and on the classical root systems). For convenience, we 
shall suppose: 

• n>2 when considering SU(n), SU(n)/SO(n), SU(2n)/USp(n) or SU(n)/S(U(n - q) x U((j)); 

• n > 3 when considering USp(7i), USp(n)/U(n) or USp(n)/(USp(n - g) x USp(g)); 

• n > 10 when considering'' SO(n), S0(2n)/U(n) or SO(n)/(SO(n - q) x SO(g)). 

This ensures in particular that the root systems and the Schur functions of type B, C and D are not 
degenerate, and later this will ease certain computations. For Grassmanian varieties, we also suppose by 
symmetry that q< [^J. 

• special unitary groups SU(n,C): a maximal torus is 



T = I diag(zi,Z2, . . . ,2;„) 



Vi, Zi g T and J| = 1 I = T"/T (2.8) 



and the Weyl group is the symmetric group 6„. The simple roots and the fundamental weights, 
viewed as elements of i* , are = — 6*+^ and tUi = + • • • + for i e |1, n — 1], where e* 
is the coordinate form on t = iR" defined by e*(diag(iti, it2, • ■ • , i^n)) ~ ti- Here there is a small 
subtlety as in t*, t37„ = + • • • + e" =0. Therefore, the fundamental weights can also be written 
as 

m^ = nji~ -run^ - — -{e^ + he*) - -{e'+^ + h e"); 

n n n 

this is actually important in order to compute correctly the Casimir coefficients. The dominant 

weights are then the 

(Ai - A2)roi H h Xn-iUJn-i = Aie^ H h A„_ie""^ - |A| — , 

n 

where A = (Ai > A2 > • • • > A„_i) is any partition (non-increasing sequence of non-negative 
integers) of length (n — 1); it is then convenient to set A„ = 0. The half-sum of positive roots is 
given by 2p = 2{mi + • • • + Wn-i) = X]"=i(" + 1 — 2i)e*, and the scalar product on t* is ^ times 
the usual euclidian scalar product ^e* | e-'^ = Sij. So, 

^A ^ -Q - Xj+j-i 



l<2<j<n 



6 



For SU(2n)/USp(n) and SO(2n)/U(n), the restriction holds on the "2n" parameter of the group of isometries. So for 
instance, we shall consider the spaces SO(2n)/U(n) with 2n > 10. 
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and the Casimir coefficient is kx — + ^ A| + (n + f — 2i)Xi. The character of the 

representation of highest weight A is the Schiir function 

det(2;- ■')i<ij<„ 

where zi, . . . , z„ are the eigenvalues of k. 

unitary groups U(n, C): though we have chosen to examine only the Brownian motions on simple 
Lie groups or symmetric spaces, the same work can be performed over the unitary groups, which 
have a non-trivial and one-dimensional center, and therefore are not semi-simple but reducible. 
The theory is in fact almost the same as for SU(n, C). Irreducible representations of U(n, C) are 
labelled by sequences A = (Ai > • • • > A„) in Z", the action of the torus T" on a corresponding 
highest weight vector being given by the morphism A(zi, . . . , z„) = Zj^^ • • • z^". The dimensions 
and characters are given by the same formulas as before, and the Casimir coefficient is this time 

spaces of real structures SU(n, C) /SO(n, M): the involutive automorphism defining the symmetric 
pair is dig) = g. The intersection of the torus (2.8) with SO(n,M) is (Z/2Z)"/(Z/2Z), the set 
of diagonal matrices diag(a;i, . . . , Xn) with coefficients Xi in {±1} and determinant xi ■ ■ ■ x^ ^ 1. 
The dominant weights A = (Ai, . . . , A„) such that \{t) = A(xi, . . . , Xn) — ■ ■ ■ x^^" = 1 for 
every t E (Z/2Z)"/(2/2Z) are those with even parts: 

\/i e , X^ = mod 2. 



• spaces of quaternionic structures SU(2n, C)/USp(n,EII): 6{g) — J2n9 J2ni where J2n is the (2n x 
2n)-skew symmetric matrix 

/O 1 \ 



J2n — 



-1 



V 



1 

0/ 



The intersection of the torus (2.8) (with 2n entries) with USp(n, H) is the set of diagonal matrices 
diag(zi, , . . . , Zn, with all the ZiS in T. The dominant weights A trivial on T H X are 

then characterized by the condition 

Vi e |l,n] , \2i-1 = \2f 

In other words, all the parts are doubled. 

complex Grassmannian varieties SU(n, C)/S(U(n — q, C) x 11(9, C)): this time, it is a little sim- 
pler to work with U(n, C)/(U(n — C) x U(g,C)), which is the same space. The involutive 
automorphism defining the symmetric pair is then 9{g) = Kn^q g Kn.q, where 



T 



K„ 



in-2q 



and Tq is the {q x g)-anti-diagonal matrix with entries 1 on the anti-diagonal. Notice that this 
does not correspond to the usual embedding' of U(7i — g, C) x U((?,C) into U(n,C); but this is 
the good choice for so that the usual maximal torus T = T" satisfies the crucial hypothesis 
"TOP maximal torus in P". The intersection of the torus (2.8) with K is then the set of diagonal 



^The usual embedding corresponds to the matrix ^ = diag(l, . 
; it is conjugated to the matrix K„^q proposed. 



,1,-1, 



-1) with {n — q) entries 1 and q entries 
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matrices diag(zi, . . . , Zg, Zg+i, . . . , Zn^q, Zq, . . . , zi) with the z^'s in T. The dominant weights A 
trivial on T n X are then 

A = (Ai, . . . , Aq, 0, . . . , 0, — Ag, . . . , — Ai). 

Therefore, they are given by a partition of length q. 

compact symplectic groups USp(n,IHI): considering USp(n,IHI) as a subgroup of SU(2n, C), a 
maximal torus is 



T = I diag(zi,Zi \ . . 



1 ^n. 1 1 



Vi, z,£T} (2.9) 



and the Weyl group is the hyperoctahedral group f)„ — (Z/2Z) I (5„, acting by permutation of the 
Zi's and by exchange of z^ and z^^. The simple roots, viewed as elements of t*, are = e* — e'+^ 
for z £ |1, n — 1] and Q!„ — 2e"; and the fundamental weights are Wi = + ■ ■ ■ + ioi i G |1, «]. 
Here, e'(diag(iti, — iti, . . . , ii„, — ii„)) = U. The dominant weights can therefore be written as 
Aie^ + • • • + A„e", where A — (Ai > A2 > • • • > A„) is any partition of length n. The sum of 
positive roots is 2p — J^^^ii'^^ + 2 — 2«)e*, and the scalar product on t* is times the usual 
euclidian scalar product. Hence, 

j^x _ TT K - Xj+j - i T-r At + Aj + 2n + 2 - z - j 

l<i<j<n ■' i<i<j<n 

and the Casimir coefficient is ka = 5^ S"=i Af +(2n+2— 2i)Ai. The character of the representation 
of highest weight A is the Schur function of type C 

X l/Cj -SCAiZl,Zi ,...,Z„,Z„ j - , „-j + l _ („-j + l). 

UtLl^^j Zj )l<i,j<n 

where zf^,...,z^^ are the eigenvalues of k viewed as a matrix in SU(2n, C). 

spaces of complex structures USp(n, IHI)/U(n, C): the involutive automorphism associated to the 
symmetric pair is 0{g) = g (again, inside SU(2n,C)). This corresponds to the choice C = M® jM 
for the inclusion C C H and then U(n,C) C USp(7i,H). The intersection of the torus (2.9) with 
U(n,C) is the set of diagonal matrices diag(xi,a;i, . . . ,a:„,x„) with coefficients Xi in {±1}. The 
dominant weights A = (Ai, . . . , A„) such that A(t) — A(xi, Xi, . . . , a;„, a;„) = x^^ ■ ■ ■ x^" ~ 1 for 
every t G (Z/2Z)" are those with even parts: 

Vie A, = mod 2. 



quaternionic Grassmannian varieties USp(n,IHI)/(USp(n—g,IHI)xUSp((7,IHI)): 9{g) — i2n,g5^2n,gi 
where 



-'2n,q 



\ 



T4 appearing q times (with all the computations made inside SU(2n,C)). Again, this peculiar 
choice is made so that Tn P is maximal in P (of rank q), and it does not correspond to the usual 
embedding of USp(n— (7, H) x USp(g, H) into USp(n, H). The intersection of K with the torus (2.9) 
is the set of diagonal matrices diag(zi, zj~^, zi, . . . , Zg, z~^, Zq-^, Zq, Z2q+i, Z2q^i, ■ ■ ■ , z„, z~^) 



with the Zi's in T. The dominant weights (Ai, 

A = (Ai, Ai, . . . , Xq, Xq, 0, 



qi "q T "q 

, A„) trivial on K HT write then as 



,0). 



Thus, they are given by a partition of length q (with parts doubled). 
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special orthogonal groups SO(n,IR): odd and even special orthogonal groups do not have the 
same kind of root system, and on the other hand, SO(n,M) is not simply connected and has for 
fundamental group Z/2Z for n > 3. So in theory, the arguments previously recalled apply only 
for the universal cover Spin(n). Nonetheless, most of the results will stay true, and in particular 
the labeling of the irreducible representations; see the end of [BD85, Chapter 5] for details on 
this question. In the odd maximal torus in S0(2n + 1, M) is 

Re. = & -r«f ) e S0(2,M)| (2.10) 



T= |diag(i?e,,...,i?e„,l) 

and the Weyl group is again the hyperoctahedral group io„. The simple roots are ai — e* — e*+^ for 
i € [1, n — 1], and = e"; and the fundamental weights are Wi = + • • • + for i £ \\,n — 1], 
and Wn = \{e'^ + ■ ■ ■ + e"). Here, 

(diag((i',70>---'(a"ro")>0))=a. 

and it corresponds to the morphism diag(i?0j , . . . , 1 1) '"^ 6'^' . The dominant weights are then 
Aie^ + ■ • • + A„e", where A — (Ai > A2 > ■ ■ • > A„) is either a partition of length n, or an 
half-partition of length n, where by half-partition we mean a non-increasing sequence sequence 
of half-integers in N' = N + 1/2. The sum of positive roots is 2p = X]"=i(2't- + 1 — 2i)e^ and the 
scalar product on t* is ^2n+i) times the usual euclidian scalar product (e* | S"*) — ^ij- 

j^x _ -Q - Aj + j - i -Q \i + \j +2n + l - i - j 



^ - 2n + 1 - i - 7 

l<i<j<n l<i<j<n 

and the Casimir coefficient is k\ = 2n+i Tll^i ^1 + (2"- + 1 ^ 2i)Xi. The character of the repre- 
sentation of highest weight A is the Schur function of type B 

X {k)-Sbx{Z,,Z, ,...,Z„,Z„ ,1)- . „-, + l/2 _ -(„-j + l/2). ' 

where , . . . , z^^, 1 are the eigenvalues of k. In the even case, a maximal torus in S0(2n, M) is 



T= \ diag{Rg^,...,ReJ 



V*. Re. = & "cr.f) e S0(2,M)| (2.11) 



and the Weyl group is the subgroup of jo„ of index 2 consisting in signed permutations 
with an even number of signs —1. The simple roots are ai = e' — e'^^ for i e |l,n — 1] and 
a„ = e"^^ -I- e"; and the fundamental weights are Wi = -\- ■ ■ ■ + for i e |l,n — 2] and 
'^n-i.n — |(e^ -I- • • • -I- e"""'^ ± e"). The dominant weights are then Aie^ + • • • + A„_ie"^^ +eA„e", 
where e is a sign and (Ai > • • • > A„) is either a partition or an half-partition of length n. The 
sum of positive roots is 2p = ^^^i{2n — 2z)e', so 

_ TT A; - Aj + j - ^ A., + Aj + 2n - ^ - j 
i — i 2n — i — j ' 

l<i<j<n •' 

and K\ = 1^ X^iLi + ('^^ ^ 2i)A,j. The irreducible characters are the Schur functions of type C 
X^(fc) = sdx{zi,z^^, . . . , z„, z^^) 

^ det(z;^-+-^ - ^^^'^^"~■'^<^.<» + det(z^+"-^ + zr''^^"'-")i<^.<n 
det(zr'+^r^"~'^)i<^.<n 
where zj*^^, . . . , are the eigenvalues of k. 

real Grassmannian varieties SO(n, M)/(SO(n — (7,M) x SO(q,M)): in this case and for the next 
one, one cannot directly apply Theorem 13, since SO(n,IR) is not simply connected. A rigorous 
way to deal with this problem is to first look at quotients of the spin group Spin(n); of course 
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this will not be very practical, as Spin(n) has no "natural" matrix representation. Nonetheless, 
consider the Grassmannian variety of non- oriented vector spaces 

Gr*(n, (7,M) ~ Spin(n)/(Spin(?i — q) x Spin(g)); 

Gr(n, M) is a 2-fold covering of Gr* (n, q, M) . The defining involution of Gr^ {n,q,R) corresponds 
to the involution of SO(n,M) given by 9{g) = Nn^q g Nn^q, where 

/T2 \ 



n,q 



To 



V 



In-2q/ 

with q blocks T2. With T as in Equations (2.10) and (2.11), the intersection K (1 T is then 
(Z/2Z)« X (SO(2,K))LtJ-9, so the dominant weights trivial on KnT write as 

A= (Ai,...,A„0,...,0) ; Vze |1,(7], A, = 0niod2. 

They are therefore given by an integer partition of length q, with all parts even. Now, for 
the simply connected Grassmannian variety Gr(n, g,IR), there are twice as many spherical rep- 
resentations, as n T is in this case isomorphic to ((Z/2Z)'?/(Z/2Z)) x TLtJ^?, instead of 
(Z/2Z)« X TLtJ-9. Therefore, the condition of parity IS now 

Vi,jG |l,g], A, EE Aj mod 2. 

Thus, spherical weights correspond to partitions of length q with all parts odd, or all parts even. 



spaces of complex structures S0(2n,M)/U(n, C): denoting — ), a convenient choice for 

e and therefore the embedding U(n,C) ^ S0(2n,M) is 6{g) = M^gM'^, where 



^U4 



if n is even 



M„ 



C/4 



C/4 



V 



if n is odd 



I2J 



with each time m — [^J blocks U4. We skip here the technicalities related to the non simple- 
connectedness of S0(2?i,E); again they can be solved by considering the space as a quotient of 
Spin(2n). In the end, the dominant weights A trivial on the intersection KnT are given by 

A = (Ai, Ai, . . . , A„i, Am) or A = (Ai, Ai, . . . , A„i, A„i, 0) 

that is to say a partition with all non-zero parts that are doubled. 

Remark. In each case, the "degree of freedom" in the choice of partitions labeling the irreducible or 
spherical representations is exactly the rank of the Riemannian variety under consideration, that is to 
say the maximal dimension of fiat totally geodesic sub-manifolds. 

Remark. The special embeddings detailed in this paragraph were only useful to find out which highest 
weights were corresponding to spherical representations (via Cartan-Helgason's Theorem V,i, which re- 
quires K to be chosen so that P H T is maximal in T). However, from now on, one can come back to 
more traditional embeddings, in particular in the cumbersome case of Grassmannian varieties. 



In the following, we denote: 2)„ the set of partitions of length n; 3n the set of non-decreasing sequences 
of (possibly negative) integers; i2)„ the set of partitions and half-partitions of length n; 22)„ the set of 
partitions of length n with even parts; 22)„ ffl 1 the set of partitions of length n with odd parts; and 2)2)n 
the set of partitions of length n and with all non-zero parts doubled. It is understood that if i is too big, 
then Xi = for a partition or an half-partition A of prescribed length. 
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Theorem 15. The density of the law lit of the Brownian motion traced on a classical simple compact 
Lie group writes explicitly as follows: 

Ae2)„-i \i<«<j<« 

Ae3„ \l<i<J<n * / 

5; e-As„.w.»«-»).,(' n ^^^^^ n 

Aei«g„ \i<»<j<n •^'^ / 

respectively for special unitary groups SU(n, C), unitary groups U(n,C), symplectic groups USp(n, H), 
odd special orthogonal groups S0(2n + and even special orthogonal groups S0(2n,K). In this last 

case, sX = (Ai, . . . , A„_i, — A„), and it is agreed that sdx + sd^x stands for sdx i/ A„ = and therefore 
sX = X. 



We denote generically <px (x) a zonal spherical function associated to a spherical representation ( the 
function depends of course of the implicit type of the space considered). The density of the law /it of the 
Brownian motion traced on a classical simple compact symmetric space writes then as follows: 



V- -jI; (-^l^+Er=i'A?+(n+i-2i)A.) / -p^ Xi-Xj+j-i\ , . , 

Ae2a)„_i \l<i<j<n / 

^ ^-^(-l^+E?.r^A?+(2„..-2,A.)/ ^ A.-A,^+i-A ^^^^^ 
Ae2)2)2„-i \l<i<j<2n ^ * / 

Aeff), \l<i<j<" ^ * / 

A622)„ \1<»<J<" l<i<j<n 

^ e-sfe ^?=i^'+(2n+2-2i)A, I -Q " + J " ^ "Q A, + A, + 2n + 2 - i - j I 

Ae2)2)2, \l<i<i<n l<i<i<n 
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y e-^2T.L^^l+i2n+l-2^)xJ TT - + J - i TT X^ + X, + 2n + 1 - I - j \ 

/ ^ lAX n — ^ -L-L 2/7, -1-1 — - % — 1 j 

Ae22),U2'21,ffll \l<i<j<n l<i<j<ri J 

y e-a^^?=i^'+(2„-2.)A. I TT (A, - A, + J - z)(A, + A, + 2n - ^ - j) \ ^ 

1 (j — i)(2n — i — j) I 

Ae22),u2a),ffll yi<j<i<ji V /V J/ y 

for real structures SU(n, C)/SO(?i, M), quaternionic structures SU(2n, C)/USp(rt, H), complex Grassman- 
nian varieties Gr(n, g, C), complex structures USp(n, H)/U(n, C), quaternionic Grassmannian varieties 
Gr(n, g,IHI), complex structures S0(2n. M)/U(n, C), odd real Grassmannian varieties Gr(2n + 1, g, M) and 
even real Grassmannian varieties Gr(2n, 

Remark. In the case of complex Grassmannian varieties, it is understood that A„+i_i = — A^ as explained 
before. We have not tried to reduce the expressions in the previous formulas, so some simplifications can 
be made by replacing the indexing sets of type 22)p or 2)2)p by 2}p. 

Example (Brownian motions on spheres and projective spaces). Let us examine the case q = 1 for 
Grassmannian varieties: it corresponds to real spheres §"(M) = SO(n + 1, ]R)/SO(n, M), to complex 
projective spaces P"(C) = SU(n + 1, C)/S(U(n, C) x U(1,C)) and to quaternionic projective spaces 
P"(IHI) = USp(n-|- 1, Il)/(USp(n, H) x USp(l,IHI)). In each case, spherical representations are labelled by 
a single integer fc S N. So, the densities are: 

,r«(.) = fe-^^ fcl±f (2fc + n- 1) (2.12) 



Pt 



"(C) 



oo 



(n- l)!fc! 

Hk+^)t ((n - 1 + fc)!)2 



(,) . ^ e--^ + ") (2-13) 



k=0 

oo 



r^(m, N fc(fc+2„+i)t (2n + fc)! (2n - 1 + fc)! , , „ ^ ,h . n , 

pf^°'(-)=Ee (2.Vl)!(2n-l)!(fc + l)!fc! (^^ + ^-+1) <t>lM- (^-l^) 

In particular, one recovers the well-known fact that, up to the aforementioned normalization factor 
(n -1- 1), the eigenvalues of the Laplacian on the n-sphere are the k{k -I- n — 1), each with multiplicity 
^^^±f (2fc + n - 1); see e.g. [SC94, §3.3]. 



Example (Torus and Fourier analysis). Take the circle T — U(1,C) = S^(M). The Brownian motion on T 
is the projection of the real Brownian motion of density i^{9) = by the map 6 i~> e'^. Thus, 

pW') = '^- E Pf(0 + 2m.) = y^ g e-^^^=y^5(0,t). 

m— — CO m— — oo 

The series S{9, t) is smooth and 27r-periodic, so it is equal to its Fourier series X^i^-oo '^kiS{t)) e*^'^, with 
Thus, the density of the Brownian motion on the circle with respect to the Haar measure ^ is 

oo oo 

pj{e">)^ e-'*e^"^ = l + 2^e-'^' cosfc^, 

A:— — OO k—1 

Since S(fc)(e^^) — c^^^, this is indeed a speciahzation of the second formula of Theorem 15, for U(l, C). 
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Example (Brownian motion on the 3-dimensional sphere). Consider the Brownian motion on USp(l,H), 
which is also SU(2,C) by one of the exceptional isomorphisms. The specialization of the first formula of 
Theorem 15 for SU(2,C) gives 



if e^'^ are the eigenvalues of g e SU(2,C). It agrees with the example'^ of [Lia04a, §4]. and also with 
Formula (2.12) when n — 3, since the group of unit quaternions is topologically a 3-sphere. 

Remark. The previous examples show that the restrictions n > no are not entirely necessary for the 
formulas of Theorem 15 to hold. One should only beware that the root systems of type Bi, Ci, Di and 
D2 are somewhat degenerated, and that the dominant weights do not have the same indexing set as for 
Bti>2 or C„>2 or D„>3. For instance, for the special orthogonal group S0(3,M), the only positive root 
is e^, and the only fundamental weight is also e^. Consequently, irreducible representations have highest 
weights ke^ with k £ N; the dimension of the representation of label fc is 2fc + 1, and the corresponding 
character is again jf ^^^^ g-ie ^^j^g ^j^g non-trivial eigenvalues of the considered rotation. So 



if (7 is a rotation of angle 9 around some axis. 



^Notice that the space is dilated in Liao's paper by 2it, whence the factor (27r)^ = 47r^ between our time exponents and 
the ones given by Liao. 
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3. Upper bounds after the cut-off time 

Let be a probability measure on a compact Lie group K or compact symmetric space G/K, that is 
absolutely continuous with respect to the Haar measure 77, and with density p. Cauchy-Schwarz inequality 
ensures that 

4dTv(A^, = \p{x) - 1| dx^ < \pix) - dx=\\p- 

The discussion of Section 2 allows now to relate the right-hand side of this inequality with the harmonic 
analysis on X. Let us first treat the case of a compact Lie group K. If one supposes that p is invariant by 
conjugacy, then Parseval's identity (2.2) shows that the right-hand side is J^xeR \x^{P ~ However, 
by orthogonality of characters, for any non-trivial irreducible representation of _ftr — i.e., not equal to 
Ik k G K i-^ 1 — one has 

x'(l)= / x\k)dk^ [ x\k)x^Hk-')dk = 0. 

J K JK 

On the other hand, for any measure /i on the group, X^'^{i^) — /kX"'"!^) I^(dk) — fi(dk) — 1. Hence, 
the inequality now takes the form 

where the / indicates that we remove the trivial representation from the summation. Similarly, on a com- 
pact symmetric space G/K, supposing that p is A'-invariant, Parseval's identity reads \\p — M\%'2(^q/k:) ~ 
'ExeG'< \^^(P - 1)1^- However, for any non-trivial representation A, 



a 



Indeed, using only elementary properties of the Haar measure, one sees that 1(A) — jQP^{g)dg = 0, 
because it is a projector and it has trace x^(l) — 0- So again, the previous inequality can be simplified 
and it becomes 

xeG^ 

In the setting and with the notations of Proposition 14, a bound at time t on 4dxv(Mtj Vk)'^ (respectively, 
on -idTvipt^Vo/K)^) is then 

^g-t(A+2p|A)(^A^2 . respectively, ^ e'^^^+^H ^> 

Proposition 16. In every classical case, 4c?Tv(Mii Haar)^ is bounded by X^Aeiv ^n('^) e^*^"'^^, where 
the indexing sets Wn and the constants yl„(A) and i?„(A) are recalled in the table on page 30. 

Remark. One can shorten all these formulas by introducing the modified weight coordinates 

• Pi = Xi + n+^z^ in type A„_i, i.e., for SU(n, C); 

• Pi — Xi + n + ^ ^ i in type B„, i.e., for S0(2n + 1, M); 

• Pi ~ Xi + n + 1 ~ i in type C„, i.e., for USp(n, H); 

• Pi — Xi + n — i in type D„, i.e., for S0(2ri, M). 
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Then, 



is equal to 



V(^l,^2,---,M«) 

V(l,2,...,n) 



in type A, 



in type B„; 



n— 1 : 



(2n- 1)!! j V((2n- 1)2,...,12) 



nl J V(n2,...,l2) 
V(/zf,...,MD 



in type C„; 



in type D„; 



V((n-1)2,...,02) 



where Y{xi, . . . , Xn) denotes the Vandermonde determinant. Also, i?„(A) is always an afRne function 
(depending on n) of J27=i A^l ^YP^ C and D; and of J2i=i - iJ2i=i ^-i) i'^ tyP^ ^■ 

This section is now organized as follows. In §3.1, we compute the weights that minimize i3„(A); they 
will give the correct order of decay of the whole series after cut-off time. In §3.2, we then show case-by-case 
that all the other terms of the series Sn{t) of Proposition 16 can be controlled uniformly. Essentially, we 
adapt the arguments of [Ros94, Por96a, Por96b], though we also introduce new computational tricks. As 
explained in the introduction, the main reason why one has a good control over Sn (t) after cut-off time is 
that each term Tn{X,t) — yl„(A) e^* of the series Sn{t) stays bounded when t = tcut-off', for every n, 
every class A and in every case. We have unfortunately not found a way to factorize all the computations 
needed to prove this, so each case will have to be treated separately. However, the scheme of the proof 
will always be the same, and the reader will find the main arguments in §3.2.1 (for symplectic groups 
and their quotients), so he can safely skip §3.2.2-3.2.4 if he does not want to see the minor modifications 
needed to treat the other cases. 



Remark. The bounds obtained on the series Sn{t) in §3.2 are not at all optimal, and in particular one 
can conjecture that the exponent of n in these bounds can at least be multiplied by a factor 2. A possible 
way to improve our bounds would be to use the following alternative approach to the control of the series 
Sn{t); it is related to the problem of minimization of the logarithmic potential 



of a probability measure v on the real line (c/. [ST97]). Set t = {1 + e) icut-off, and consider only the 
case of compact groups. The crucial remark is that up to some explicit constant fnit), the general term 
Tn{X,t) of the series Sn{t) is equal to the exponential of a logarithmic potential of either the discrete 
measure X]"=i '^m; ca,se of unitary groups, or of the discrete measure Y^^=i ^it^i)^ case of 

special orthogonal groups or compact symplectic groups (the /i'^s are the modified weight coordinates): 




T„(A,t) cx l/(/xi,...,^„)2e'"^'"=i^? when K = \J{n); 




r„(A, t) cx V{pl, . . . , nl)^ e-5^ ^«"=i when K = S0(2n). 
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Suppose for a moment that the /Xi's are arbitrary (ordered) real numbers, instead of integers of half- 
integers. Then one can use the following well-known results on discrete minimizers of logarithmic poten- 
tials, see [Sze39, Chapters 5 and 6] and [MOT53, Volume 2, Chapter 10]. Denote 

the Hermite polynomials, and 



(x)^ — — fe"^ a;"-i) • L(°'> (x) = — — 



the generalized Laguerre polynomials of weig ht -1 and 0. If ^ E^^i < L, then the maximum of 
A(/ii, . . . , /in)^ is attained if and only if the /ij's are the n zeroes of the polynomial 




Similarly, under the same assumption ^ X^ILi /^l — maximum of (UlLi /^O ■ ■ • ' Mn) at- 

tained if and only if the are the n zeroes of the polynomial Ln \nx/L), and the maximum of 
A(/i^, . . . , /i^J is attained if and only if the /x^'s are the n zeroes of the polynomial Ln ^\in — l)x/L). In 
each case, these values can be calculated, see [Szc39, Theorem 6.71]. 



This enables one to have an explicit bound on r„(A,t) assuming that X^ILi contained in some 
interval [a,/3]. Multiplying these bounds by ^ times the number of integer of half-integer points in the 
"orange peel" 



(xi , . . . , Xji ) G 



E 



a;, 



and then summing these bounds over k provides an upper bound n^^'^ on the part of S'„(t) correspond- 
ing to weights such that X^iLi l^i enough, say bigger than 10 it^ — we skip the details of these 
computations, as they are quite similar to what will be done in §3.2. Unfortunately, we were not able 
to use this method to bound the whole series Sn{t): indeed, for "small" values of X]"=iMii points 
of the lattice of weights cannot approximate sufficiently the roots of the aforementioned renormalized 
Hermite or Laguerre polynomials, and therefore the corresponding bound was too crude for most of the 
corresponding terms in Sn{t). Nonetheless, this alternative method hints at possible better bounds on 

Sn{t). 



3.1. Guessing the order of decay of the dominating series. Remember the restriction n > 2 
(respectively, n > 3 and n > 10) when studying special unitary groups (resp., compact symplectic groups 
and special orthogonal groups) and their quotients. In what follows, we use the superscript ★ to denote 
a set of partitions or half-partitions minus the trivial partition (0, 0, . . . , 0). The lemma hereafter allows 
to guess the correct order of decay of the series under study. 

Lemma 17. Each weight Xmin indicated in the following table corresponds to an irreducible representation 
in the case of compact groups, and to a spherical irreducible representation in the case of symmetric 
spaces of type non-group. The table also gives the corresponding values of An and Bn- In the group case, 
Bni^min) is minimal among {_B„(A), A £ W*}. 
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K or G/K 


-^min 


('^min) 


('^min) 


S0(2n + 1,M) 


(1,0,. 


.,o)„ 


271 

2n+l 


(2n + l)2 


S0(2n,]R) 


(1,0,. 


.,o)„ 


2n-l 
2ri 


At? 


SU(n,C) 


(1,0,.. 


,0)„_i 




2 


USp(n, H) 


(1,0,. 


.,o)„ 


2n+l 
2n 




Gi{2n + l,q,R) 


(2,0,. 


..,0), 


2 


2n2 + 3n 


Gr(2n,g,E) 


(2,0,. 


..,0), 


2 


2n2 + n - 1 


Gr(n, g, C) 


(1,0,. 


..,0), 


2 


- 1 


Gr(n,g,H) 


(1,1,0, 


• •,0)29 


2 


(n- l)(2n + 1) 


S0f2n R]/l](n C) 


(1,1,0, 


..,o)„ 


2(n-l) 
n 


n(2n — 1) 


SU(n,C)/SO(n,M) 


(2,0,.. 


,0)„_i 


2(n-l)(n+2) 


n(n+l) 
2 


SU(2n,C)/USp(n,H) 


(1,1,0,.. 


• , 0)2n-l 


(ri-l)(2n+l) 


n(2n- 1) 


USp(n,H)/U(n,C) 


(2,0,. 


■,0)„ 


2(n+l) 
n 


n(2n + 1) 



Remark. For symmetric spaces of type non-group, one can also check the minimahty of -B„(A,nin), except 
for certam real Grassmannian varieties Gr(n, (7,IR). For instance, if q = 1, then (1)^ labels the geometric 
representation of SO(n,M) on C", which has indeed an invariant vector by SO(n— 1,M) x S0(1,M); and 
the corresponding value of B{X) is (n — l)/n < 2. Fortunately, Amin, though not minimal, will still yield 
in this case the correct order of decay of the series S{t). 



Remark. To each "minimal" weight Amin corresponds a very natural representation. Namely, for a 
special orthogonal group SO(n, M) (respectively, a compact symplectic group USp(n,IHI)), the mini- 
mizer is the "geometric" representation over C" (respectively C^") corresponding to the embedding 
SO(n,R) ^ SO(n,C) ^ GL(n,C) (respectively USp(n,Il) ^ SU(2n,C) ^ GL(2n,C)). For a spe- 
cial unitary group SU(rt,C), one has again the geometric representation over C", and its compose with 
the involution k i— > (fc*)^^ corresponds to the label (1, . . . , l)„_i, which also minimizes B„{X). The case 
of spherical minimizers is more involved but still workable: we shall detail it in Section 4. 

Proof. To avoid any ambiguity, we shall use in the following indices to precise the length of a partition 
or half-partition. Let us first find the minimizers of -Bn(A) in the group case: 



SU(7i): one has to minimize 



n 



J2X^ + {n + l-2^)X,^U ^ (A,-A,)M 

i=l \l<i<j<n I 



^i{n-i){X,-X,+i) \ =A + B 



over 2)*_2- In B, at least one term is non-zero, so 

B > min i(n — i) 
Ug[i,„_i] 
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with equality if and only if A = (1,0,..., 0)„_i or A = (1, . . . , l)„_i. In both cases, A is then 
equal to However, is also the minimum value of A over 2)*_]^. Indeed, there is at least 
one index I G |1, n — 1] such that A; > A/+i. Then all the (Ai — A^)^ with i < I and j > I + 1 give 
a contribution at least equal to 1, and there are l{n — I) such contributions. Thus 

l(n — I) n — 1 

A > > , 

n n 

and one concludes that mini?„(A) is obtained only for the two aforementioned partitions, and is 

equal to ^(A^in + S,nin) = 1 - 

S0(2n): the quantity to minimize over ^2}^ is 

^ i{2n - 1 - i){Xi - A,;+i) j + n(n - l)A„_i = A + B + C, 
1=1 J 
again with A, B and C non-negative in each case. Only A involves A„, so a minimizer satisfies 
necessarily A„ — (partitions) or A„ = ^ (half-partitions). In the case of partitions, a minimizer 
of -|- C is (1,0,..., 0)„, which gives the value min^gii i{2n — 1 — i) = 2rt — 2. The same 
sequence minimizes A over SI)*, so the minimal value of A + B + C over non-trivial partitions 
is 2n — 1 and it is obtained only for (1,0,..., 0)„. On the other hand, over half-partitions, the 
minimizer is (^j---!^)^! giving the value 

n n(n — 1) n(2n — 1) 

4 ^ 2 ^ 4 ' 
Since we assume 2n > 10 and therefore n > 5, this value is strictly bigger than 2n — 1, so the 
only minimizer of i?„(A) in ^2}* is (1,0,..., 0)„. 

• S0(2n + 1): exactly the same reasoning gives the unique minimizer (1, 0, . . . , 0)„, with corre- 
sponding value 2nfoT A + B + C = {2n + 1) i3„(A). 

• USp(ri): here one has only to look at partitions, and the same reasoning as for S0(2n) and 
S0(2n -|- 1) yields the unique minimizer (1, 0, . . . , 0)„, corresponding to the value 2n + 1 for 
2nB„(A). 

The spherical minimizers are obtained by the same techniques; however, some cases (with n or g too 
small) are exceptional, so we have only retained in the statement of our Lemma the "generic" minimizer. 
The corresponding values of An{X) and i3„(A) are easy calculations. □ 



Suppose for a moment that the series S'„(t) of Proposition IG has the same behavior as its "largest 
term" A„(Aniin) e^*'^"'-'^"""'. We shall show in a moment that this is indeed true just after cut-ofF time 
(for n big enough). Then, Sn{t) is a O(-) of 

n^e"' for SU(n), SO(n), USp(n); 

n^e-^* for SU(n)/S0(7i), SU(2n)/USp(n), S0(2n)/U(n), USp(7i)/U(n), 

SU(n)/S(U(n - g) x \J{q)), S0(n)/(S0(n ~ q) x SO{q)), USp(n)/(USp(n - q) x USp(g)). 

Set then tn,e = a (1 + e) logn, with a = 2 in the case e^*, and a = 1 in the case e^^*. Under the 
assumption Snit) ^ An{\min) s^* '^"^^"""^ , one has 

Sn{tn,e) = 0{n-^'). 

Thus, the previous computations lead to the following guess: the cut-ofF time is 

• 2 log n for classical simple compact Lie groups; 

• logn for classical simple compact symmetric spaces of type non-group. 
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3.2. Growth of the dimensions versus decay of the Laplace-Beltrami eigenvalues. The estimate 
Snitn,E) ~ ^n(Amin) e~*" = ^'>('^) — 0{n^^'^) might seem very optimistic; nonetheless, we are going to prove 
that the sum of all the other terms v4.„(A) e^*" ' ^"^■^') in Sn{t) does not change too much this bound, and 
that one still has at least 

We actually believe that at least in the group case, the exponent 2e is good, cf. the remark before §3.1 
— the previous discussion shows that it is then optimal. 

Suppose that one can bound v4„(A) e^*" ' by c(n)l'^l, where |A| is the size of the partition (sum of 
all parts) and c(n) is some function of n that goes to as n goes to infinity (say, Cn^^'^). We can then 
use the following: 

Lemma 18. Suppose x < ^. Then, the sum over all partitions X^A'^''^'' "which is convergent, is smaller 
than 1 + 5a;. Consequently, 

J2 < 5x. 

A#(0,...,0) 

Proof. The power series P{x) — x'^' — YliLi rrjrr = ^ + x + 2x^ + Sx"^ + Sx"* + • • • has radius of 
convergence 1, and it is obviously convex on M+. Thus, it suffices to verify the bound at x = and x — ^. 
However, 

P(0) = 1 = 1 + (5 X 0) ; P (^^ < 3.463 < 1 + (^5 x 

□ 



With this in mind, the idea is then to control the growth*^ of the coefficients yl„(A), starting from the 
trivial partition (0, . . . , 0). Of course, the arguments will be slightly different in each case, and we will 
have to adapt them to the case of half-partitions, even partitions, partitions with parts doubled, etc.; 
however, the reasoning will be the same each time. The simplest cases to treat in this perspective are 
the compact symplectic groups and their quotients. 

3.2.1. Symplectic groups and their quotients. Set tn,e = 2(1 + £)logri; in particular, i„,o = 21ogn. We 
fix a partition A € and for fc < A ji. W6 denote pk,7i 

the quotient of the dimensions associated to 

the two rectangular partitions 

{k,...,k)n and (fc - 1, . . . , fc - 1)„. (3.1) 
Using the formula given in §2.3 in the case of compact symplectic groups, one obtains: 

_ jj 2k + 2n + 2- i- j _ yr 2 ( \^ ^ 

2fc + 2n - i - 7 ^ 2k + 2n-i-j " \ ^ 2k + 2n - i ~ j 

l<i<j<n ■' l<i<j<n ■' \\<.i<j<n 

The double sum can be estimated by standard comparison techniques between sums and integrals. 
Namely, since x,y ^ ik+in-x-y convex on {(x,?/) | x > 0, y > 0, 2fc + 2n > x + j/}, one can 
bound each term by 

2 /■/■ 2 , , 

< / / — ax ay. 



2k + 2n-i~j i,i+i]x[j-i j+i] 2fc + 2n - x - y 

We use this bound for non-diagonal terms with indices i < j, and for diagonal terms with i = j, we use 
the simpler bound 

n ^ n—1 ^ ^ 

E .TT^ -T^kTu- ^^--^"^ - ^ fc + + - - 1) - log 

i=l u=0 



This is also what is done in [Por96a, Por96b], but the way we make our partitions grow is different. 
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where Hn denotes the n-th harmonic sum. So, 

l0gPfc,„< V ^. , ^ ^ : : < Hk+n-l - Hk-l + 1 1 ^ dx dy 

l<i<j<n ■^■^[2-'"+2\ " 

< — + logffc + 7T, — 1) — log k 

k 

+ {2k + 2n - 1) log(2fc + 2n - 1) + {2k - 1) log(2A: - 1) - 2(2A: + n - 1) log(2fc + n-l). 

On the other hand, the same transformation on partitions makes — t„.o Bn{X) evolve by —{2k + n) logn. 
So, if rj^ ^ is the quotient of the quantities e"*" ^ ^"(-*') with A as in Equation (3.1), then 

2k -\- n 1 
log 77fc,n < ^ — log " + ^ + log(^ + n - 1) - log fc 

+ {2k + 2n-l) \og{2k + 2n - 1) + {2k - 1) log(2fc - 1) - 2(2fc + n - 1) log(2fc + n-l). 

Suppose k > 2. Then, one can fix n > 3 and study the previous expression as a function of k. Its 
derivative is then always negative, so logryfe.„ < log 772, „, which is also always negative. From this, one 
deduces that 

for any rectangular partition (A„, . . . , A„)„; indeed, the left-hand side is the product of the contributions 
rjk.ii for k in |1,A„]. However, iji^n is also smaller than 1: in this case, the dimension is given by the 
exact formula 



so rji^n = Cat„+i e-T^ which can be checked to be smaller than 1 for every n > 3. So in fact, 

for any rectangular partition (A„, . . . , A„)„. 

The previous discussion hints at the following more general result: 
Proposition 19. In the case of compact symplectic groups, at cut-off time, 

- 3 

for any integer partition A of length n (not only the rectangular partitions). 

Proof. We fix Z S |1, 71 — 1], and the idea is again to study the quotient pk.i of the dimensions associated 
to the two partitions 

(fc + A;+i, . . . , A; + A;+i, A;+i, . . . , A„)„ and (A; - 1 + A;+i, . . . , fc - 1 + A;+i, A;+i, . . . , A„)„, (3.2) 

where k is some integer smaller than A; — A/+i — in other words, the n — I last parts of our partition 
have already been constructed, and one adds k to the I first parts, until k — Xi ~ Xi+i; see Figure 4. 

The transformation on partitions described by Equation (3.2) makes the quantity —tn QBn{X) change 

by — ^^^'^ "tt^""''' logn. We shall prove in the following that this variation plus logpfc.; is almost always 
negative. For convenience, we will treat separately the cases ^ = 1 or 2 and the case I > 3; hence, suppose 
first that I E |3,ri — 1|. The quotients of Vandermonde determinants can be simplified as follows: 

- TT k + j-l + Xi+i-Xj k + Xi+i + Xj+2n+l-j yr 2k + 2Xi+i + 2n + 2 - i - j 



.^^^^ k+j-l-l + Xi+i -Xj k + Xi+i + Xj+2n + l-j-l 2/c + 2A/+i + 2n - i - j 
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A/+1 + k 



A;+i 



12 3 ••• I l + l ■■■ n 



Figure 4. One makes the partitions grow layer by layer, starting from the bottom. 



Notice that the second product Pk.i,{2) this formula is very similar to Pk,n] the main difference is that 
indices i,j are now smaller than I (instead of n). Hence, by adapting the arguments, one obtains 

\<i<i<i ■' •'"'[s-'+sJ " 

< ^ h log(fc' + n - 1) - log(fc' + n-l) + (2k' + 2n - 1) log(2fc' + 2n - 1) 

k' + n — I 

+ {2k' + 2n - 2/ - 1) log(2fc' + 2n~2l-l)- 2{2k' + 2n-l-l) \og{2k' + 2n-l-l) 

where k' stands for k + A;+i. So, if (rik.i)^ is the quotient of the quantities (D^)^ g-t>i,o B„(x) ^j^^]^ ^ g^g 
Equation (3.2), then log?7fej < log%^; + logp^ where log%j is given by 

l(2k' + 2n-l)^ 1 , s , 

- ^ log 71 + — + log(fc' + n - 1) - log(fc' +71-0 

+ {2k' + 2n - 1) log(2fc' + 2n - 1) + {2k' + 2n-2l-l) log(2fc' + 271-21-1) 

- 2{2k' + 2n-l-l) log(2fc' + 2n - ? - 1), 

and Pk,i.(i) is the first product in the expansion of pk^i- Let us analyze these two quantities separately. 

• logTT/j /: here the technique is really the same as for log77fc „. Namely, with n and I fixed, logTy^ ; 
appears as a decreasing function of a; = fc', because its derivative with respect to x is 

llogn 1 1 1 

2 _ _| 

n {x + n — ly^ X + n — 1 x + n — I 
+ 2(log(2a: + 2n - 1) + log(2a; + 2ji - 21 - 1) - 21og(2x + 27i - / - 1)). 

A upper bound on the first line is —U—lllSMl < g (remember that n > 3 and therefore logn > 1), 
and the second line is negative by concavity of the logarithm. From this, one deduces that 
logrjk^i < log^i^j, and we shall use this estimate in the following in order to compensate the other 
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part of logryfe,;: 

, „ ^ l{2v + 2 + 2n - Z) , 1 , , . , , , n 

logr?fe,; < logn + — — — , + log(w + n) - log(v + n + l-l) 

+ {2v + 2n + 1) log(2u + 2n + 1) + {2v + 2n - 2Z + 1) log(2u + 2n - 2/ + 1) 
- 2(2t; + 2n-l + l) \og{2v + 2n-l + l) 
where v stands for A;+i. 

logPfe,;,(i): in the product Pk,i,{i), each term of index j writes as 

(k' + nf - (Aj + n + 1 - if ^ (k' + nf - (A,+i + n + 1 - j)^ 



(fc' + n - ly - (Aj + n + 1 - j)2 - (/c' + n - 0^ - (A,+i + n + 1 - j)^ 

^ + j - 1 fc" + 2n + 1 - j 
~ fc + j - Z - 1 k" -\-2n+\- 3 - I 
with fc" = A; + 2Ai+i = fc + 2v. and multiplying all these bounds together, one gets 
(fc + n-1)! ik-\)\ {k" + 2n-l)\ (k" + n - 1)1 
- {k + l- 1)! (A; + n-Z- 1)! (A;"+n)! (/c" + 2n-20!' 
Again, this is decreasing in k, so 

n! (2t; + 2n - Z + 1)! {2v + n-l + 1)! 
- Z! {n - /)! (2i; + n + 1)! {2v + 2n-2l + 1)! " 
Recall the classical Stirling estimates: for to > 1, 

log to! = TOlogTO+ ^logm- TO + log\/27r + r„, with < r„ < 

2 12to SdOto-^ 

It enables us to bound log Pk,i,{i) by the sum of the following quantities: 

★ A={2v + 2n-l + l) log(2?; + 2n - I + I) + {2v + n - I + 1) \og{2v + n - I + 1) 

-{2v + n + l) log(2w + n+l) - {2v + 2n-2l + l) log(2i; + 2n - 2/ + 1). 

-k B= l{\og{2v + 2n-l + l)+ log(2v + n-l + l)- \og{2v + n + 1) - log(2?; + 2n - 2/ + 1)), 
which is non-positive by concavity of the logarithm. 

■k C = n log n — I log I — [n—l) log(n — I). 

■k D = 2 (logn — logZ — log(n — I)). This is non-positive unless n = I + 1 — recall that we 
assume for the moment / € |3, n — 1]. In that case, it is smaller than ^(n-i) ■ 

. p _ X f k _ 1 5_ _i 1 I 1 1 1 \ 

12 ; n-l 2v+2n-l+l ^ 2u+n-J+l 2v+n+l 2v+2n-2l + l J' 

* ~ 3^0 (f + (n-;)3 + (2v+n+l)3 + (2w+2n-2i+l)3 ) ' 

The sum of the two last terms EF = E + F happens to be negative. Indeed, E and F are 
decreasing in v (we use the convexity of a; ^ to show that ^ < 0), so it suffices to check the 
result when v = 0. Then, with I fixed, 

,x 1 /I 1 1 1 1 1 1 

EF{n, I) = 



12 \n I n-l 2n-l + l n-l + 1 n + 1 2n-2l + l 
1/11 1 1 



360 (n-Z)3 (n+l)3 (2n-2Z + l)3, 

is decreasing in n, hence smaller than its value when n = l + l. So, it suffices to look at EF{1 + 1, 1), 
which is now increasing in /, but still negative. Thus, in the following, we shall use the bound 

logpk,i,(i)<A + C + D<A + C+^J—^. 
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Adding together the bounds previously demonstrated, we get 

^ l{2v + 2 + 2n-l) ^ 1 1 , / ^ , / 1 n 

log r]k,i < log n + H ■ — + log u + n) ~ log i; + n + 1 - I) 

2n 2n — 2 v + n + 1 — I 

+ {2v + 2n + l) log(2w + 2n + 1) - (2u + 2n - / + 1) log(2w + 2n-l + l) 
+ {2v + n-l + l) log(2u + n - I + 1) - {2v + n + 1) \og{2v + n + 1) 
+ n log n — I log I — {n — I) log(n — I). 

By concavity of xlogx, the sum of the second and third rows is non-positive. What remains is decreasing 
in / and in v, and when I ~ 3 and w = 0, we get 

^^<^Zn+ + -^+\og(^^] +(n-3)log (^^) -SlogS 
2n 2n~2 n~2 \n — 2J \n — 3 J 

which is maximal for n — 5, and still (barely) negative at this value. Thus, we have shown so far that 
'>lk,i < 1 for any k, any Z g |3, n — 1], and any partition A that we fill as in Figure 4. 



When Z = 1 or / = 2, the approximations on log% ; that we were using before are not good enough, 
but we can treat these cases separately. When 1 = 1, 



Pk,l 



X2 + k + n -Q fc + j - 1 + Aa - k + X2 + Xj + 2n + 1 - j 

J=2 



< 



A2 + A: + n - 1 J- fc + ?■ - 2 + A2 - A,- fc + A2 + A,- + 2n - 7 

j=2 

k + n -A- fc + j - 1 fc + 2n + 1 - j k + 2n - 1 



k + j - 2 k + 2n- j 

j=2 



k 

If /c = 1, which only happens once when one makes the partition grow, then the bound above is 

2ti + 1 1 

2ne fs-'°sn < 2. On the other hand, if fc > 2, then the bound is decreasing in k and therefore 
smaller than (n -I- i) e"^ < 1. So, one also has rj^^i < 1 for any k but k = 1, where a correct 
bound is 2. Similarly, when 1 = 2, 

_ Xs + k + n 2A3 + 2fc + 2n - 1 -A- fc + j - 1 + A3 - Aj fc + A3 + Aj + 2n + 1 - j 
^'''^ ^ A3 + fc + n - 2 2A3 + 2fc + 2n - 3 fc + j - 3 + A3 - Aj fc + A3 + A^ + 2n - 1 - j 

^ k + n 2fc + 2n - 1 -A- fc + j - 1 fc + 2ri + 1 - j _ fc + 2n - 2 fc + 2n - 1 2fc + 2n - 1 
~ fc + n- 2 2fc + 2n - 3 fc + 7 - 3 fc + 2n - 1 - 7 ^ fc fc+ 1 2fc + 2n-3' 

fc + 2n - 2 fc + 2n - 1 2fc + 2n - 1 2,1+2^-2 

Vk 7 < e " . 

fc fc + 1 2fc + 2ri-3 

Again, the last bound is decreasing in fc, smaller than 2 -I- ^ < | when fc = 1 and smaller than 1 when 
k = 2. Hence, rik,2 < 1 unless fc = 1, where a correct bound is | (and again this situation occurs at most 
once whence making the partition grow). 



Conclusion: every quotient rjkj satisfies rjk^i < 1, but the two following exceptions: fc = 1 and Z = 1 or 2. 
The product of the bounds on these two exceptions is 2 x | = so for every partition A, one has indeed 

1=1 fe=l 



□ 
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Remark. A small refinement of the previous proof shows that the worst case is in fact the partition 
(2, 1, 0, . . . , 0)„ — by that we mean that any other partition has quotients pk^i that are smaller. Its 
dimension is provided by the exact formula 

. ^ Snjn^ - 1) 
3 

so one can replace the bound -y of Proposition 19 by |. 

Remark. The main trick in the proof of Proposition 19 is the way we make our partitions grow. One 
could have tried something simpler, namely, make the first part Ai grow box by box, then the second 
part A2, etc. However, with the same technique of estimation of quotients of dimensions, we were then 
not able to prove something better than 

instead of Ce 2^-^"(^) — this is in accordance with the bound 2a log n mentioned in SaloS-Coste's 
Theorem 7. By making the partitions grow layer by layer, we use in a much better way the fact that the 
parts Ai, A2, . . . , A„ are ordered, and the "compensations" in the growth of given by the determinantal 
structure of its formula. 



The upper bound (1.11) is now an easy consequence of Lemma 18 and Proposition 19. For any partition 
A, notice that 



^ n ^ n 1 ^ 

Bn{X) > — ^(2n + 2-2i)A, = — ^i(2n + l-i)(A, -A,+i) > -J2i{X^ ~ X^+l) 



i—1 i — 1 i—1 

From this, one deduces that in the case of compact symplectic groups, 

if one assumes that < | (in order to apply Lemma 18). By Proposition 16, one concludes that 

,USp(n,H). -rj \ ^ 3 

"TV lM2(i+e)iogn, naarj S — ■ 



Here one can remove the assumption ^ < otherwise, the right-hand side is bigger than 1 and therefore 
the inequality is trivially satisfied. This ends the proof of the upper bound in the case of compact 
symplectic groups. For their quotients, one can still use Proposition 19, as follows. For quaternionic 
Grassmannians, 

assuming -\ < \. This implies that 

,Gr(n,g.H)/ „ \ / ^ 

"TV (/^(i+e)iogn,Haar) < — . 

Again, the assumption on ni is superfiuous, since otherwise the right-hand side is bigger than 1. Exactly 
the same proof works for the spaces USp(7i)/U(n), with the same bound (it may be improved by using 
the fact that one looks only at even partitions). 
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3.2.2. Odd special orthogonal groups and their quotients. Though the same reasoning holds in every case, 
we unfortunately have to check case by case that everything works. For odd special orthogonal groups 
S0(2n + 1, M), set tn^e = 2 (1 + e) log(2n + 1), with in particular tnfi = 2 log(2n + 1). The main difference 
between S0(2n + 1) and USp(n) is the appearance of half-partitions, which is solved by the following 

Lemma: 

Lemma 20. For any integer partition X, denote A ffl | the half-partition Ai + |, A2 + |, . . . , A„ + i. 



V(Sn(Affli)-S„(A)) < gn(log2-'25(^) < ^ 



Proof. The quotient of dimensions is 

TT A, + A, +2/^ + 2-^-j ^ yr 2n + '^-i-J 
.,4 , \i + \j+2n+l-i-j- J-} 2n+l-i-j 

and the difference (i?n(A ffl |) — B„(A)) is equal to 

log(2n + l) /-A 1 2n+ 1 - 2A log(2n+ 1) /n n^X nlog(2n+l) 



+ 1) , 1 2n+l-2i\ 



2n+l l^' 4 2 /- 2n + l V4 2 

\i=l / 

This yields the first part of the inequality, and the second part is an easy analysis of the variations of the 
bound with respect to n. □ 



Then, for any integer partition A, one can as before prove a uniform bound on e~ '°8(2n+i) 
the differences are tiny, e.g., in many formulas, 2n + 2 is replaced by 2n + 1, or ^ is replaced by ^n+i ■ 
So more precisely, with the same schema of growth of partitions, one has the following bounds: 

• ?7i^„: it is given by the exact formula {^^^-l) e~ ^n+i^ iog(2Ti+i)^ which is indeed smaller than 1 for 
n > 5. 

• 'nk>2,n' the comparison techniques between sums and integrals give 

n(2k — 1 + n) 2 1 

logr/fc,„ < 2n + l — ^"^^^"^ ^^^^ 2k^ ^ k ~^ + n - 2) - log fc 

+ {2k + 2n-2) log(2/c + 2n - 2) + {2k - 2) log(2fc - 2) - 2(2fc + n - 2) log(2A; + n - 2). 

This bound is decreasing in k, whence smaller than its value when k = 2, which is negative for 
every value of n > 5. 

• Vk,lel3,n-iY- as before, pk,i splits into Pk,i,{i) and Pfe,;,(2): 

n /v +.;-! + A;+i - A^- /.: + A/+1 + A^ +2/) -j y-r 2k + 2A;^i + 2// + 1 - < - j 
.^^^^ k + j-l-l + Xi+i-Xj k + Xi+i + Xj+2n-j-l i<j-<,.<, 2fc + 2A,+i + 2n - 1 - i - j ' 

The bound on log?7fc,/, the sum of logpk,i,{2) and of the variation of — ^^B„(A), is 

, _ Z(2fc'- l + 2n-0 , 1 , „^ , , 

log%,i < ^— ^ log(2n + 1) + — + log(fc' + n - 2) - log(fc' + n - I - 1) 

+ {2k' + 2n-2) log(2fc' + 2n - 2) + {2k' + 2n - 2Z - 2) log(2A:' + 2n - 2/ - 2) 

- 2(2A:' + 2n-l-2) log(2/c' + 2n-l-2) 

l{2v + 2n+l-l) , 1 , , ,x , . 

< — ^ , log(2n + 1) + , + log(t; + n - 1) - log(w + n-l) 

2n + 1 V + n— I 

+ {2v + 2n) \og{2v + 2n) + {2v + 2n- 21) \og{2v + 2n- 21) 

- 2{2v + 2n-l) log(2t; + 2n-l) 
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with k' = k + A;+i = k + V. On the other hand, in the product pk,i,{i), each term of index j 
writes as 

(A:' + n - 1/2)2 - (Aj + n + 1/2 - j)2 ^ ^j.^ + n - l/2f - {Xi+i + n + 1/2 - jf 



{k' + n- 1/2 - 0^ - (A,- + n + 1/2 - j)2 - {k'+n- 1/2 - l)^ - (A,+i + n + 1/2 - 

k + j-1 k" + 2n- j 



< 



k + j-l-lk" + 2n-j-V 
so the quantity pk,i,(i) is bounded by 

(fc + n-1)! (A:" + 2n - / - 1)! {k" + n - I - 1)\ 



{k + l-iy. {k + n-l-iy. (fc"+n-l)! (&" + 2n-2/-l)! 
^ n!(2i; + 2n-/)!(2t; + n-Z)! 



Z! (n - /)! (2u + n)! (2t; + 2n- 2l)\ ' 
Again, Stirling approximation leads to 



log Pk,i,{i) < (2v + 2n-l) log(2D + 2n-l) + {2v + n-l) log(2u + n-l)-{2v + n) log(2t; + n) 

- (2t; + 2n-201og(2t; + 2n-2/)+nlogn-Zlog/- (n - /) log(n - /) + — ^ — , 

2n — 2 

and therefore 

;(2u + 2n + l-n , 1 , , ,^ , . 

logr?fe,i < log(2n + 1) + — + log(^; + n - 1) - log{v + n-l) 

+ {2v + 2v) log(2?; + 2n) + {2v + n - I) \og{2v + 71-1) 
- {2v + n) \og{2v + n)- {2v + 2n-l) log(2i; + 2n-l) 

1 



+ n log n — I log I — [n — l) log(n — /) + 



2n-2 



< -M2!i±l_0 iog(2n + 1) + ^ + log(n - 1) - log(n - 

2n + 1 n — l 

+ n log n — / log I — {n — l) log(n — + 



2n-2 

The last bound is decreasing in I, so it suffices to look at the case Z = 3; then the bound is 
decreasing in n, so it suffices to check that the bound is negative when n = 5, which is just a 
computation. We conclude that \ogr]k,i < for any k and any i e |3, n — 1]. 

77fe,i: a bound on pk,i is ^±f^ fkXH-l ^ ^° 



nu^ < k + n-22k + 2n-l gj^±jp^ iog(2n+i) < (n-l)(2n+l) ^^iog(2n+i) 
- fc 2fc + 2n - 3 - 2n - 1 



n—l log(2n + l) 1 log 11 

< e 2n+i < _ e 11 < 1. 

- 2n - 1 - 2 



• 77fc,2. a bound on pk,2 is ^^^^ 2fc+2n-5 fe+;^' 

fc + 2n — 4 fc + 2n— 3 2k + 2n — 1 k + n— 1 4fc+4,i-6 1„„/9„ , i n 4iog(2,i+i) 

' ' - k k + 1 2k + 2n-5 k + n-2 - 2n+l 

The last bound is bigger than 1 only when n = 5 or 6. The maximal value is obtained for n = 5, 
and is smaller than 1.09 < j^. More( 
smaller than 1 even when n = 5 or 6. 



and is smaller than 1.09 < j^. Moreover, if fc > 2, then one has a much better bound, that is 



Putting all together, one sees that at most one quotient rik,i may be bigger than 1 (and actually only 
when n = 5 or 6). Thus, we have proved the following: 
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Proposition 21. In the case of odd special orthogonal groups, at cut-off time, 

11 

To 



for any integer partition A of length n. For half-integer partitions, the bound is replaced by ^ . 



There is one last computation that needs to be done, namely, the special case^" A — (i, . . . , i), 
(0, . . . , 0)„ ffl ^. The value of Bn{\) is then ^, and — 2". Thus, in this special case. 



4 (2n + l) 
for every n > 5. On the other hand, 

1 " I A I " 



(» + l)|A| > N 

- - . - - - . - 2n + 1 - 2 ' 

1=1 1=1 



so we can now write: 



4 (2n + l)^ ^^^^ 

< 11 1 ^ |^^^A-)2g-t„,oB,.(A) _^ pAffli^2g-t„.oB„(Affli)~j g-2elog(2n+l)S„(A) 



11 




1 


T 


(2n 


+ 1)^ 


11 




1 


T 


(2n 


+ IY 


11 




1 


T 


(2n 


+ iY 


11 




1 


T 


(2n 


+ 1)^ 




1 


<i 


(2n 


+1)= 



< i^!: ^ I V f — + — ^ e-^l^l iog(2n+i) 



^ V 100 25 / 

121 1 33 144 

< \ > — < < 

" 4f (2n + l) 



Assg 



,SO(2ri+l,K1/ X , 6 



"TV (,M2(l+£)log(2n+l):-naarJ s 



(2n + 1)5 



and again we can now remove the assumption (2w+i)g — 5' 1"-^^^ same technique applies to odd real 
Grassmannians, with 

S ( \ - V n^c-^^"(^) < — V c-5l^l'°g(2«+i) < ^1 < 1^ 

"12;,^^, -10^ - 10(2n + l)f - (2n+l)i' 

^ ^ Ae(2ig,u22),ffli)* Ae2)* \ ' ) \ < j 

and therefore 

2 



jGr(2n+l,(3,R) / tt \ ^ 

"TV vA^(i+e) iog(2n+i)i -tiaarj s 



(2n + l)t 



In the next paragraph, we shall prove exactly the same bounds for even special orthogonal groups and 
even real Grassmannians. 



It corresponds to the spin representation of the odd orthogonal group SO(2n + 1, M). 
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3.2.3. Even special orthogonal groups and their quotients. For even special orthogonal groups, denote 
tn,e = 2 (1 + e)log(2n). The possibility of a sign ± for the last part A„ of the partitions leads to a 
coefficient 2 in the series Sn (t) , and on the other hand, the case of half-partitions is reduced to the case 
of partitions by way of an analogue of Lemma 20. Indeed, 

- - 5 

for any n > 5 and any partition. That said, let us analyze as before the various quotients pk,i and r]k^i 
corresponding to the growth of partition described by Equation (3.2): 

• 77fc „: the general formula is 

/n-l . 

2k + n-i-l 2fc + n-i-2 



_ I -p-r 2fc + 2n - 2i - 1 2fc + 2n - 2i - 2 \ _ 2fc+r.-2 iog(2„) 



\ 2 = 1 

which is decreasing in k and reduces to (^"^^) e ^ ^ when k = 1. This latter bound is always 
smaller than 1. 

%,;el2,«-i]: the quotient of dimensions pk.i = Pk,L(i} Pkd,{2) is equal to 

fc + j - 1 + A;+i - Aj fc + Ai+i +Aj +2n- 1 - j -p-r 2fc + 2A;+i + 2n - i - j 



fc + j - / - 1 + Ai+1 - Aj fc + Ai+1 + Aj + 2n - 1 - j - ? ^<|^^.<;2fc + 2Az+i + 2n - 2 - i - j 

The main difference with the previous computations is that Pk,i,(2) is a product over distinct 
indices i < j, so we will not have to worry about diagonal terms in the corresponding sum (see 
the argument at the beginning of §3.2.1). Hence, with the same notations as before^^, 

log% / < + iog(2„) + (2fc' + 2n - 3) log(2fc' + 2n - 3) 

2n 

+ {2k' + 2n - 2/ - 3) log(2fc' + 2n - 2? - 3) - 2{2k' + 2n - ? - 3) log(2A:' + 2n - ^ - 3) 

< _ ^(2^' + 2n-0 Yog{2n) + {2v + 2n-l) \og[2v + 2n-l) 
2n 

+ {2v + 2n -21-1) log(2n - 2? - 1) - 2{2v + 2n-l-l) \og{2v + 2n-l- 1); 
log Pk,i.(i) <{2v + 2n-l-l) \og{2v + 2n - I - I) + {2v + n - I - 1) \og{2v + n - I - I) 
- (2u + n - 1) log(2w + n - 1) - (2u + 2n - 2? - 1) log(2u + 2n -21-1) 

1 



+ n log n — I log I — {n — I) \og{n — I) + 



2n-2' 



Adding together these bounds, using the concavity of xlogx and then the decreasing character 
with respect to v gives 

_ l{2n — l) 1 
log r]k,i = log 77fc,/ + log Pk,L(i) < T) log(2n) + nlogn - llogl - [n - I) \og[n - I) + 



2n °' ' ° o V y ' 2n-2 

which is decreasing in Z > 2. Then, 

2n — 2 1 

log(2n) + n log(n) - 2 log 2 - (n - 2) log(n - 2) H 

n 2n — 2 

is decreasing in n, and one can check that it is negative when n — 5. So, rj^^i < 1 for any k and 
any ;e |2,n-ll. 



^"'^We make the convention that log = 0, since 2v + n — l — 1 might be equal to in the following (when v = A;^_i = 
and I = n — 1). 
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?7fc,i: one has pk,i < ttn-l ' therefore 

fc + 2n — 3fc + n— 1 _ 2n+2fc-3 
VkA < r T— ^ e 2„ 



log(2n) 



Suppose k > 2; then the right-hand side is smaller than ^^2n^' Vk,! < 1- On the other 

hand, for k = 1, which happens only once, 

log(2n) log 10 4 

Vi,! < e 2„ < e 10 < -. 

o 

Thus, we have shown that: 

Proposition 22. In the case of even special orthogonal groups, at cut-off time, 

D^e 2-^'^W<_ (^respectively, — ) 
for any integer partition (resp. any half-partition) A of length n. 



Besides, the same proof as in the case of odd special orthogonal groups shows that Bn{\) > ^ for 
any partition. For the special half-partition A = (0, . . . , 0)„ ffl | that cannot be treated by combining 
Lemmas 18 and 20, one has Z?'*' — 2"~^ and i?„(A) = j, hence 

/£,AN2 ^t„,.S„(A) ^ (n^l)log4-^iiS|^ log(2„) ^ 1 

^ ^ - - (2n)^ 

for n > 5. We conclude that 



< ^ + V + — e-^^l^l < < 



(2n)'^ V 9 225 y " 45(2n)^ " {2n^) ' 



and therefore, by Proposition 16, 



d^°(^"'*'(^2(i+£)iog(2n),Haar) < 



(2n)5 

For even real Grassmannian varieties, 

2/ , "3,4^ - 3(2n)f - (2n)§' 

and again, the total variation distance is bounded by 2/(2n)3. So, the inequalities take the same form 
for even and odd special orthogonal groups or real Grassmannians, and the proof of the upper bound in 
this case is done. The same inequality holds also for the spaces of structures S0(2n)/U(n). 

3.2.4. Special unitary groups and their quotients. Set tn,e 2(1 + e)logn. For special unitary groups, 
Weyl's dimension formula fortunately takes a much simpler form than before, but on the other hand, the 
computations on Bn{X) are this time a little more subtle. We shall still prove that almost every quotient 
Vk.i of the quantities e-*'<^,o BnW -v^rith A going from 

(Ai+i + fc — 1, . . . , A;+i + fc — 1, Xi+i, . . . , A„_i)„_i to (A;+i + k, . . . , A/+1 + k, . . . , A„_i)„_i 

is smaller than 1; but in practice, what will happen is that the negative exponentials may be much larger 
than before, whereas the quotients of dimensions pk.i will be much smaller. Consider for a start rjk^n-i- 
One has 

-p-i- k -\- n — i k + n — 1 
Pk,n 1 IL jf_i^ri — i k ' 

i—l 
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(ra-l)(n+2fc-l) 



whereas i3„(A) is changed by — — '-^^ -. So, 

k + n-1 (^-i)(»+2fc-i) f ne"'^'°s" = e^ < 2i if fc = 1, 



Vk.n-l 



k I ii±ic-^^^a?^i°g « < iLLi < 1 iffc>2, 



by using the decreasing behavior with respect to k. Notice that is indeed much smaller than before 

(linear in n whereas before it grew exponentially in rt), but i?„(A) for fc = 1 is almost constant instead of 
linear in n. 



In the general case, 

n k' - Aj + j - 1 ^ -A- fc + j - 1 
k'-\^+j-l-l- J-J- k + j-l-l 

3=1 + 1 J + l ■> 

with the usual notation fc' = fc + A;-|_i. On the other hand, the transformation on partitions makes -B„(A) 
change by 

-l(n - l)(n + 2fc' - 1) + 2^|A|;+i,„ 
where |A|i+i.„ is the restricted size '^v Notice now that 



-(n-Ofc' + |A|/+i,„ = ^ Aj-A;+i-fc< ^ -k ^ -{n- l)k. 
j=i+i 



So,, 

n 

r,. , < TT 

k+j-l-l ~ \K 

J=l+1 ^ ' 

which can as usual be estimated by Stirling (this is the same kind of computations as before). Hence, 
with / > 3, the last bound is always smaller than 1, and also if Z = 2 unless n = 4. If n = 4 and I — 2, 
then 

(fc + 2)(fc + 3) 3^i„g2 iffc = l, 

fc(fc + l) -\\ iffc>2. 

Finally, when I — \, one has exactly the same bound as for / = n — 1, so 2i when fc = 1 and 1 for fc = 2, 
Multiplying together all the bounds (3/25 and twice 23), we obtain: 

Proposition 23. In the case of special unitary groups, at cut-off time, 

- 2 

for any integer partition A of length n — 1. 



An other big difference with the previous cases is that one cannot use Lemma 18 anymore. Indeed, 
for A — (fc, . . . , fc)„_i, Bn{\) = = so there is no hope to have an inequality of the type 

5n(A) > fc |A| for any partition. That said, set 5i — \i — Xi+i] then, 

B„(A) = ^ (A,-A,f + i5:^(n-0<5,>;^^^,. 

1 < t < J < n z— 1 i— 1 

This leads us to study the series 

n— 1 ^ 

(5i,...,(5„_i>o i=i -L — a: " 
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Clearly, each Tn{x) is convex on IR+, so if we can show for example that r„ (|) stays smaller than 
1 + ^ for every n, then we will also have the inequality Tn{x) < 1 + Kx for every < a; < g. Set 
Un{o^ = log(T„(a;)); one has 



= 5^ - log (l - ar^) < 5^ < 2 ^ ; 

i—1 i—1 i—1 



for < x < |. It follows that T„(x) <1 + Kx with K < 169. Suppose < |. Then, 



which leads to 



,SU(n,C)/ TT \ / 10 

"TV (M2(i+e) log n. -naar) < — . 



If ^ |, then this inequality is also trivially satisfied. Hence, the case of special unitary groups is 
done. For the quotients SU(n)/SO(n), one obtains 



and therefore 



,<'f^"f-)-ll<-<- 



,SU(n,C)/SO(Ti,M). TT \ ^ ° 

dry (M(i+e)iogn,Haar) < — . 

71/ 2 



The proof is exactly the same for SU(2n)/USp(n) and gives the same inequality, however with (2n)2 
instead of ni. 



For the complex Grassmannian varieties, we have seen that it was easier to see them as quotients 
of U(r?,) (instead of SU(n)), and this forces us to do some additional computations. Though the cut- 
off phenomenon also holds in the case of U(n), the set of irreducible representations is then labelled 
by sequences of possibly negative integers, which makes our scheme of growth of partitions a little 
bit more cumbersome to apply. Fortunately, for Grassmannians, the spherical representations can be 
labelled by true partitions, but then the dimensions are given by a different formula and we have to do 
once again the estimates of quotients pk,i and r]k,i- For a partition of size p = [^\, one has B„(A) = 
f ELi A? + (n + 1 - 2i)Xi and either 



An{X) 




Xi + Xj + n + 1 - i - j 
n + 1 — i — j 




if n = 2p, or 




n + 1 — i — j 

when n = 2p + 1. Let us give the details when n = 2p. Again, one looks at pk,l = -^nW / -^nif^) and 
Vk,l = pfc ,c-i°s"('^"(^)-'^"(^)), with /i and A equal to 

(A(+i + fc - 1, . . . , Xi+i + k - 1, Xi+i, . .. , Xp)p and (A;+i + k,.. ., Xi+i + k, A;+i, . . . , Xp)p. 

The quotient of dimensions is 



Pk,l 



{2k' + n-j){2k' + n-j-l) \ I -A- [k' - X, + j - l){k' + Xj + n - j) 




n (2A;' + n-j-0(2fc'+n-j-Z-l) I I .f,^-^ {k' - X^ + j - I - l){k' + X^ + n - j - I) 
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and a lower bound is then obtained by the usual replacement A;+i = = and then k = 1: 

n-2l + l + 
Pk,i < — j — , 

71+1 \ / 

This leads to the inequality 

The last quantity is decreasing in I, as the quotient of two consecutive terms of parameters n, I and n,l + l 
is smaller than 

'n + l- l 



So, 



e-^i°g"| <i. 

1j ~\~ 1 



ry../<'^(n + lfe-2i°*5" = ^<l. 
n + 1 rt^ 



It follows that An{X) e i°g"^"(^) is smaller than 1 for any partition (we leave to the reader the verification 
of the other case n = 2p + 1, which is very similar). Finally, one can compare directly Bn{X) to |A|: 

S„(A) - ^ E A? + (n + 1 - 2i)K > 2 ^^^^(A, - A,+i) > *(A, - A,+i) = |A|. 

i—1 i—1 i—1 

We conclude that 

5„ f%^) < Y: < ^ < ; 4f ■'''^HM(i+s)iog.,Haar) < 4 

\ 2 J ^ ^ / 6 ^2 

and this ends the proof of all the upper bounds of type (1.11). 
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4. Lower bounds before the cut-off time 

The proofs of the lower bounds before cut-off time rely on the following simple ideas. Denote Amin the 
(spherical) irreducible representation "of minimal eigenvalue" identified in Section 3.1. We then consider 
the following random variable: 

^ I x'^'"'"(fc) in the case of groups, 

^ I Vf^'"'" (j)^""" (gK) in the case of symmetric spaces of type non-group. 

In this equation, k or gK will be taken at random either under the Haar measure of the space, or under 
a marginal law fit of the Brownian motion; we shall denote Eqo and Ef the corresponding expectations. 
When f2 is real valued, we also denote Varoo and Varf the corresponding variances: 

VarP = E[n^] - E[n]^ = E[(rj - E[rj])2] . 

In the case of unitary groups and their quotients, 51 will be complex valued, and we shall use the notations 
Varoo and Var^ for the expectation of the square of the module of Q — E[Q]: 

Ya.T[n] = E[|f7p] - \E[n]\^ = E[|r2-E[r2]|2] . 

The normalization of Equation (4.1) is actually chosen so that f2 is in any case of mean and variance 1 
under the Haar measure. 

Remark. In fact, much more is known about the asymptotic distribution of these functions under Haar 
measure, when n goes to infinity; see [DS94]. For instance, over the unitary groups, the moments of 
order smaller than uq of (s) = trg agree with those of a standard complex gaussian variable as 

soon as n is bigger than tiq. In particular, if g is distributed according to the Haar measure of U(n, C), 
then tr g converges (without any normalization) towards a standard complex gaussian variable. One has 
similar results for orthogonal and symplectic groups, this time with standard real gaussian variables. 



One will also prove that under a marginal law fit, the variance of f2 stays small for every value of t, 
whereas its mean before cut-off time is big (not at all near zero). Standard methods of moments allow 
then to prove that the probability of a event 

E^ = {k\ \nik)\ > a} or {gK \ \n{gK)\ > a} 

is before cut-off time near 1 under fit, and near under Haar measure (for an adequate choice of a). This 
is sufficient to prove the lower bounds, see §4.2; in other words, $1 is a discriminating random variable 
for the cut-off phenomenon. 

The method presented above reduces the problem mainly to the expansion in irreducible characters 
or in spherical zonal functions of $7^ or of |r2p; cf. §4.1. In the case of compact groups, this amounts 
simply to understand the tensor product of y-**"" with itself, or with its conjugate when the character 

is complex valued. However, for compact symmetric spaces of type non- group, this is far less obvious. 
Notice that a zonal spherical function c/)^ of a space G/K can be uniquely characterized by the following 
properties: 

• it is a linear combination of matrix coefficients of the representation V^: 

cj^\gK)^J2J2^''pU9K)- 

• it is in ^'^{G/K)^ , i.e., it is ii'-bi-invariant; and it is normalized so that (j>^{eK) — 1. 

Consequently, if (\/^»in)«'2 = y^i q . . . q y^s q yei , , , q yet ^^^-i-^ ^-]^g yi', spherical irreducible 
representations and the V^^ non-spherical irreducible representations, then there exists an expansion 

(0^n-)2 ^ c,,(f>''^ + c^.J''^ + ■ ■ ■ + c,^^"'. (4.2) 
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Nonetheless, it seems difficult to guess at the same time the values of the coefficients Cu in this expansion^^. 

Example. Let us examine in detail the case of the complex projective line 

P^(C) = SU(2,C)/S(U(1,C) X U(1,C)) =U(2,C)/(U(1,C) x U(1,C)). 

As we shall see later, the label A,nin = (1,0,. .. ,0,-1) „ of the discriminating representation in the 
case of a complex Grassmannian varieties Gr(n, g,C) = U(n, C)/(U(p, C) x U((7,C)) corresponds to the 
adjoint representation of SU(n,C) on s[(n,C). In the following, we use the traditional embedding of 
U(p, C) X U(5, C) into U(n, C) by block diagonal matrices of sizes p — n — q and q. A scalar product 
on s[(n,C) for which SU(n, C) acts by isometry is (M | TV) = trMA^^^; and a spherical vector for the 
subgroup S(U(p,C) x U(g,C)) is 

y/nm V Piq 

As a consequence, if {gij)i<i.j<n are the coefficients of a matrix g e SU(n, C), then the spherical function 

(k(i,0;---,0;~i)n "v^rites as 



(1,0,...,0,-1)„ 



(g) = tiiMpJgMp^gg-y) = tiiMp^qgMp^qg-^) 




i=p+l i=p+i 



Me E 1^.1^ + E T.\9.n 
j=p+i i=p+ij=i J 

-\{tt\m?\+\i± E wA-i (4.3) 
^ \i=i j=i J ^ \i=p+ij=p+i J 

by using on the third line the fact that rows and columns of a unitary matrix are of norm 1. In particular, 
the random variable f2 is real-valued. Now, it can be shown (by a calculation with Schur functions) that 

(^(l,0,...,0,-l)„)®2 ^ ^(0,...,0)„ ^ ^(2,0,....0,-2)„ ^ ;^^^^-j^(1.0,...,0,-l)„ ^ i^^^^(l,l,0....,0.-l,-l)„ 
© y(1.0,....0,-l)„ ^ ;^^^^p.(2,0,...A-l,-l)„ e i^^3^(l.l,0,....0,-2)„^ 

Actually, the first line is the decomposition in irreducibles of the symmetric square 5^(sl(n,C)), whereas 
the second line is the decomposition of the skew-symmetric square v4^(s[(n, C)). The two last irreducible 
representations are not spherical, so there should be an expansion of the form 

(^(l,0,...A-l)„)2 ^ _J_ ^ ^^(1,0,...A-1)„ + J^(2,0,...,0,-2)„ ^ ^^(1,1,0,.. 

— 1 

But since one does not know a priori what are the spherical functions (^(2.o,. -,o,-2)„ 0(i,i,o,...,o,-i,-i)„^ 
it seems really hard to find by this direct algebraic approach the coefficients a, &, c. 

In the case of P^(C), one can proceed as follows. The matrices of SU(2,C) write uniquely as ('^ 
with |wp + |zp = 1, and the zonal spherical function (j)^'^^ is then 

<t)^^-^{w,z) = 2|w;|2 - 1. 

Denote H = En — E22, X = E12 and Y = E21 the generators of the complex Lie algebra s[(2,C); 
the decomposition (51(2, C))®^ = l/"'*^ ® V'^'^^ ® V^'~^ can be seen as a decomposition in irreducible 
s[(2, C)-modules. The first space is generated by the Casimir element 

2 2 
C = 2^ E,j(E) Eji - Eu (E) E^j 

i,j = l = l 



^^The only "easy" computation is the coefficient of the constant function in {4>^)^ , or more generally in a product </<^ (j>f. 
Indeed, 

^ J X otherwise. 



THE CUT-OFF PHENOMENON FOR BROWNIAN MOTIONS ON SYMMETRIC SPACES OF COMPACT TYPE 51 



and it corresponds to the constant function. The second space is generated by the five symmetric tensors 

Sa^3SiX,Y)-C ■ Sb=SiH,X) ; S, = SiH,Y) ; 
Sd^S{X,X) ; Se^S{Y,Y). 

where S{V, W) — ViS^W + Wi^V. The action of H reads then as follows: 

H ■ Sa = ; H ■ Sb ~ 2Sb ; H ■ Sc — -~2Sc ; H ■ Sd ^ 4:Sd ; H ■ Se ^ — 4S'e. 

As H generates linearly the subalgebra s(g[(l, C) x C)), it follows that the spherical vector associated 
to V^'~^ is up to a scalar constant 

Sa = Ell (E) E22 + E22 ® Ell + E12 (g) £^21 + £^21 ® £^12 - £^11 «) £^11 - £^22 ® £^22 = S{X, Y) — '- 

and therefore, (/)^'^'^{w, z) = 1 — 6|wzp = 1 — 6\w\'^ + 6|u;|'*. So, one sees that 

and in particular the spherical function does not appear in the right-hand side. Now, it seems quite 

hard to generalize this argument to the general case of sl{n, C) acting on (s[(n, C))*^^. Indeed, one would 
first need a description of each irreducible submodule^'^ inside (s[(n, C))®^, and then to find adequate 
s(g[(p, C) (E) gl(g, C))-spherical vectors in these spaces — it seems to us the only direct algebraic way to 
determine the spherical functions 0(2.0,. ..,o,-2)„ ^^^^ 0(i,i,o,...,o,-i,-i)„ ^ 



The previous example shows how difficult it can be to find the expansion (4.2); as far as we know, 
for a general zonal spherical function, there is a definitive answer only in the case of symmetric spaces 
of rank 1, see [Gas70]. For our problem, one can fortunately give in every case a geometric description 
of the discriminating spherical representation and of the corresponding spherical vector. This yields an 
expression of (f>^"™{gK) as a polynomial of the matrix coefficients of g, similar to Equation (4.3). Now 
it turns out that the joint moments of these coefficients under fj,t and Hca — Haar can be calculated 
by mean of the stochastic differential equations defining the G-valued Brownian motion; see Lemma 27, 
which we reproduce from [Lcvll, Proposition 1.4]. As {(j>'^'^''^{gK))'^ or \(j)^"'^"{gK)\'^ is also a polynomial 
in the coefficients gij, one can therefore compute its expectation under /it, and this actually gives back 
the coefficients in the expansion (4.2). Thus, the algebraic difficulties raised in our proof of the lower 
bounds will be solved by arguments of stochastic analysis. 



^■'Even the identification of the irreducible component y'^'"' inside cS^(sl(n, C)) and ^^(s((n, C)) is not easy. 
Namely, one can show that the first space is linearly generated by the symmetric tensors 

n 2 ^ 

i=i " 1=1 

n 2 " 

SEki = ^ S{Eii,Eki) ^ S(Eki, Eii), 

i — l i — \ 

whereas the second space is generated by the skew-symmetric tensors 

n n 
i=l i = l 

with A{V, W) = V^W-W^V. Thus, the spherical vectors of label (2, 0, ... , 0, -2)„ and (1, 1, 0, . . . , 0, -1, -1)„ lie in 
the orthogonal in S'^{sl{n,C)) of the subspace generated by the SHi^'s, the SEi^i's, and the Casimir element; but this does 
not really help us to determine these vectors. 
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4.1. Expansion of the square of the discriminating zonal spherical functions. The orthogonality 
of characters or of zonal spherical functions ensures that for every non-trivial (spherical) irreducible 
representation A, 



= 0; 



= 0. 



On the other hand, the function corresponding to the trivial representation, which is just the constant 
function equal to 1, has of course mean 1 under the Haar measure, and also under ^j. On the other hand. 
Theorem 15 allows one to compute the mean of a non-trivial irreducible character of zonal spherical 
function under jit- 

Et[x^]- / pf(fc)x^(fc)rffc=[x'](pf)=^^c-4^"W = {A„(A)e-*^"(^)}' 



K 







= Vd^ f 






Jx 



X=G/K 



with the notations of Proposition 16, and where [x^]{f) or [(/>'^](/) denotes the coefficient of or (/>'^ in 
the expansion of /. In particular, the table of Lemma 17 yields immediately the value of Et[r2] in each 
case. 



In order to estimate Vari [Q] and Varoo [^] , we now need to find a representation-theoretic interpretation 
of either fl^ when fl is real-valued, or of when ft is complex- valued. Let us begin with the case of 
compact groups: 

Lemma 24. Suppose G = S0{2n,R) orSO{2n + l,R) or USp(n,H). T/ien f7 = x^^'"'-'""^" is real-valued, 
and 



^2 ^ (-^(l,0,...,0),.-)2 ^ ^(2,0,...,0)„ _^ ^(1,1,0,. ..,0)„ _^^(0,0,...,0)„ 

On the other hand, when G — SU(n,C), is complex-valued, and 

|^|2 ^ ^(l,0,...,0),.-i ^(l,...,l)„-i ^ ^(2,l,...,l)„-i _^ ^(0,0,...,0)„_ 



(4.4) 



(4.5) 



Proof. In each case, f2(fc) = tr/c, up to the map (1.2) in the symplectic case; this explains why Q, is 
real- valued in the orthogonal and symplectic case, and complex- valued in the unitary case. Then, the 
simplest way to prove the identities (4.4) and (4.5) is by manipulating the Schur functions of type A, 
B, C and D; indeed, these polynomials evaluated on the eigenvalues are known to be the irreducible 
characters of the corresponding groups, see §2.3. We start with the special orthogonal groups. In type 



B„, (zi 



z„ ^ + 1)^ is indeed equal to the sum of the three terms 



l<i<j<n 



Zj + Z^ ^Zj 1 



\l<i<j<n 

sb(o,...,Q){Z, ^~\l) = 1; 



ZiZj + Zj -\- Z^ Zj 



/ n 



n; 
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whereas in type D„, (zi + • • • + z„ + z^^ + • • • + z^^)^ is equal to the sum of the three terms 
srf(2,o,...,o)(^, ^ ^ ZiZ^ + z^zj^ + zr'^Zj + z^^zj^ - n - 1; 

\l<i<i<n / 
\l<i<j<n 

sd(o,.....o){Z, Z^^) = 1. 

For compact symplectic groups, hence in type C„, (zi + • • • + 2;„ + z^^ + ■ • • + z~^)^ is indeed equal to 
the sum of the three terms 

sc(2,o,...,o)(^, Z~^) = ^ ZiZj + ZiZ~^ + z'^Zj + z~^z~^^ 

\l<i<j<n 

sc(i,i.0,...,0)(^, ^~^) = X! ZiZj + ZiZj'^ + z~'^Zj + z~'^zj'^ \ + n- 1; 

\l<i<J<n / 

sc{o,...,o){Z,Z^'^) = 1; 

and this is also (sc(i,o,...,o) ^"^))^ = {x^^-°---'°'> (k))'^ = n{k)'^. Thus, Formula (4.4) is proved. 



In type A„_i, notice that for every character x^, X'*'(^') — X^{k ^) = X'^ (k), where A* is the sequence 
obtained from A by the simple transformation 

(Ai > A2 > • • • > Xn-i)n-i ^ (Ai > Ai - A„_i > • • • > Ai - A2)„_i. (4.6) 

Indeed, if zi, . . . , z„ are the eigenvalues of k, then 



X^(fc) = ■S(Ai,...,A„_i)„_i(% \ . . . ,Z„^) = S(Ai,...,A„_i,0)„(2;i \ ---,Z„^)^ S(0 -A„_i,... -Ai)„(^n,- • ■ , Zi) 
= S(Ai,Ai-A„_i,...,0)„(2:i, • ■ ■ ,^n) = S(Ai,Ai-A,._i,...,Ai-A2)„-i(2:i: • ■ • , ^ri) = X^ (k) 

Here, one uses the relation Z1Z2 • • • z„ = 1 for every element of the torus of SU(n,C), which enables one 
to transform a n-vector of possibly negative integers into a (n — l)-vector of non-negative integers^"'. In 
particular, 

\n{k)\' = |^(1A...,0)„-.(^)|2 ^ ^(l,0,...,0)„_.(fc)^(l.l,....l).-i(fc). 

Then, a simple calculation with symmetric functions yields Formula (4.5); 
^(i,o,...,o)„_,(fc) ^k)^{z, + --- + z,,){z^' + ■■■ + z-i) 

-l + ^z,z-i-f zriz, +1 

i<j / 

= S(1,0,...,0 -1)„(^) + S(0....,0)„(^) = S(2,l,...,l),i-i(^) + S(o^...^o)„_i(^) 
^^(2,l,...,l).-i(fc)+^(0...,0)„_i(^) 

where Z = {zi, . . . , z„} is the alphabet of the eigenvalues of k. □ 



The computation is probably much clearer on the unitary group U{n, C), where we do not have to do this transformation 
of weights. 
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4.1.1. Values of the zonal functions and abstract expansions of their squares. As explained in the intro- 
duction of this part, the case of compact symmetric spaces of type non-group is much more involved. We 
start by finding an expression of D,{gK) in terms of the matrix coefficients gij of the matrix g. 

Proposition 25. In terms of the matrix coefficients of g, ^^'"'"(yif) is given by the following formulas: 



G/K 






k 


Gr(n,q,M) 


SO-L(n,C) 


^ ^ ji — 1 ^ J n — 1 iQii') ~t~ 77 ^ Ji — n-Li ^ J ^ — iQii ) 1 

rj ■• i — X * J — J. \'-' "J / q / ^ i — Pi -L ' — Pi -'- J ' 


R 


Gr(n,q,C) 


s[(n,C) 


p sr=i s^=i ifijp + g x]j=p+i isijp ~ 1 


R 


Gr(n,9,H) 


sp-L(2n,C) 




M 


S0(2n,K)/U(n,C) 




n Z^i=l Z^j = lff(2i)(2j)ff(2j-l)(2j-l) - .g(2i)(2j-l)5(2i-l)(2j) 


M 


SU(n,C)/SO(n,M) 


52(C") 


^ Sj = l(.9y)^ 


C 


SU(2n,C)/USp(n,H) 




n 1^1=1 2^3 = 1 9{2i){2j)gi2t-l){2j-l) - .g(2i)(2j-l).g(2i-l)(2j) 


C 


USp(n,H)/U(n,C) 


52(C2") 


^ + am,)? mi9^,)r iM9^,)r 


M 



For rea/ Grassmannians, so^(n, C) denotes the orthogonal complement o/so(n,C) msl(n, C); and for 
quaternionic Grassmannians, sp^(2n,C) denotes the orthogonal complement o/sp(2n,C) ms[(2n,C). 

Proof. Each space described in the statement of our proposition is endowed with a natural action of 

G = SO(n) or SU(n) or USp(n), namely, the action by conjugation in the case of Grassmannian varieties, 
and the diagonal action on tensors in the case of spaces of structures. Then, to say that 

^soKEf'"' =^0^^^) ; y^;„''c)-''-^)"=sl(n,C) ; V^^g^Jj^f " = sp^(2n, C) ; 

^^(1,1,0,...,0)„ _ A2((rn\ . -j^(2,0,...,0)„_i _ c2/^n\ . -j^(l,l,0,...,0)2„-l _ ^2/^2ri\ 

''^SO(2n,K) — ./^ ) , ''sU(n,C) — I i ''^SU(2n,C) — ) , 

-,r(2,0,...,0)„ _ c2/fn2n\ 
^USp(n,H) ) 

is equivalent to the following statements: the trace of g G S0(7i,M) acting on 50^(71, C) is given by the 
Schur function of type B or D and label (2,0, . . . ,0)liij; the trace of g G U(n,C) acting on sl{n,C) is 
given by the Schur function of type A and label (1,0,..., 0, —l)n', etc. Let us detail for instance this last 
case. We have seen in the previous Lemma that 

S(1,0,...,0 -l)„(^) = (Zl H \- Zn)(Zi^ H h z~^) - 1. 

On the other hand, the module 0((n, C) on which SU(n,C) acts by conjugation is the tensor product of 
modules (C") (E) (C")*. It follows that the trace of the action by conjugation of g e SU(n, C) on gl{n, C) 
is 

Xia) = (trg) (tr(5^')*) = (zi + • • • + z„)(zf ^ + • • • + z"!) 

if zi, . . . , z„ are the eigenvalues of g. Subtracting 1 amounts to look at the irreducible submodule sl{n, C) 
inside gl{n, C). The other cases are entirely similar, and the corresponding values of the Schur functions 
have all been computed in Lemma 24. 

Once the discriminating representations have been given a geometric interpretation, it is easy to find the 
corresponding X-invariant (spherical) vectors. We endow each space of matrices with the invariant scalar 
product {M I N) = ti MN^ , and each space of 2-tensors with the scalar product {xi (E) X2\ yi <E) 1/2) — 
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{xi I yi) {x2 I 2/2), where {v \ w) is the usual Hermitian scalar product on C" or C^". We also denote {ei)i 
the canonical basis of C" or C^". Then, the iiT-spherical vectors write as: 



Lr 


A 




SO(n) 


SO{p) X SO{q) 


1 /-g/p \ 
^/npq \ pig J 




SU(n) 


S{U{p) X U{q)) 


1 ( -qip \ 

^/npq \ pig J 




USp(n) 


USp(p) X USp(g) 


1 ( -qhp 
y/2npq \ p/2g 


) 


S0(2n) 


U(n) 




SU(n) 


SO(n) 


Si=l ® 


SV{2n) 


USp(n) 


Yll^l e2i fXl 621-1 - 62i-l (Xl 62i 


USp(n) 


U(n) 


^ E -^1 <E> 6, 



In each case, e'^°"" belongs trivially to y-**"!- and is of norm 1, so the only thing to check then is the 
ii'-invariance. In the case of Grassmannian varieties, the matrix e'^™'- commutes indeed with G{p) x G{q), 
since it is also {p, g)-block-diagonal and with scalar multiples of the identity matrix in each diagonal block. 
For U(n) inside S0(2n), we use the well-known fact that inside^'' SL(2n,C), 

S0(2n,M) n Sp(2n,C) ~ U(n,C), (4.7) 

the isomorphism being given by the map (1.1). This implies in particular that U(rt) leaves invariant the 
skew-symmetric tensor X^iLi 'S^ G2i-i — S2i-i (S> e2i corresponding to the skew-symmetric form defining 
Sp(2n, C). The intersection formula (4.7) also proves that U(n) leaves invariant the symmetric tensor 
Tliili ^ Si, whence the value of the spherical vector for U(n) inside USp(n). Finally, for SO(?t.) inside 
SU(n) and USp(n) inside SU(2n), we use again the defining symmetric bilinear form or skew-symmetric 
bilinear form associated to the group K to construct a /(T-invariant vector. The value of is then 

given by the formula (l)'^{g) = (e"^ | p^{g)e^), that is to say 

tr(Mp,,gMp,,5*) ; ti {Mp g Mp^g g^) ; ^ tr(Mp,, g Mp,, 5^) 
for real, complex and quaternionic Grassmannians; 

for SU(n)/SO(n) and USp(n)/U(n); and 

^ n 71 

„ 3(2i)(2j)3(2i-l)(2j-l) - g(2»)(2j-l)5(2«-l)(2j) 

i=l j = l 

for S0(2n)/U(n) and SU(2ri)/USp(n). Here by g we mean the complex matrix of size 2n x 2n obtained 
from a quaternionic matrix of size n x n by the map (1.2). In this last case, the computations can in fact 



Here, Sp(2n, C) is the (non-compact) complex symplectic group, that is to say the set of complex matrices of size 
2n X 2n that leaves invariant the skew-symmetric form ^2i-i A 62^. The notation USp(n,,IHI) used in this paper was 

meant to avoid any confusion between Sp(2n, C) and its compact form, the compact symplectic group. 
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be done inside M(n,IHI): indeed, 

(ff(2i-l)(2j-l))^ + (5(2i-l)(2j))^ + (g{2i)(2j-l)f + {g(2i)(2])Y 

= 2(([l](.g,,))2 + - - ([k](ft,))'), 

whereas AI* — {M)^ and | trM — 5R(tr A/). Thus, the formulas for the discriminating spherical functions 
of the spaces of structures are entirely proved, whereas for Grassmannian varieties, it suffices to check 
that for any unitary quaternionic matrix N, 

^ p p ^ n n 

mtT M,,,NM,,,N*) --EE if^^- 1' + ^ E E iff^^- 1' - 1; 

^ i=l j = l ^ i=p+l j=p+l 

indeed the real and complex cases are specializations of this formula. This last calculation is identical to 
the one done in the example before §4.1, see Formula (4.3). □ 



Lemma 26. There exists coefficients a,b,c,... (different on each line, and depending on n and q) such 
that the following expansions hold: 

Gr(n,q,M) : (cb^'-°'--'^iiA' = — ^ + a</,(^'°-'°)itJ + 6</)(^'^'°--°^LtJ 



Gr(n,g,C): (^C^.i.-.D.-i) 
Gr(n,q,H): (0(iaA...,o)„^ 
S0(2n)/U(n) : (0(i4.o,...,o)„^ 
SU(n)/SO(n) : 
SU(2n)/USp(n) 
USp(n)/U(n) : 



^^^(2,2,0,...,0)l5j +d0(34A...,O)L5j ^^^(4,0,...,0)l.j. 
2 1 



(2,1, ...,!)„- 



4,2,...,2)„- 



1 +CC 



2,2,l,...,l,0)„_i . 



2n2 - n - 1 
2 1 



(2,0,...,0)„_i 

(1,1,0, ...,0)2,.-l 

2 



2n2 — n 
2 

+ 

+ n 

2 1 



2n2 — n 



+ a0(l^O,...,O)„ _^ ^^(l*,0,...,0)„ ^ ^^(2,2,0,...,0)„. 
a0(4,2,...,2)„_i. 

,(2,2,1,.. .,l,0)2„-i. 



(2,0, ...,0), 



1 



2n2 + n 



^ ^^(2,0,...,0)„ ^ ^^(2,2,0,...,0)„ + c^(4,0,...,0)„ 



In these formulas, it is understood that if the label A of the spherical function (j)^ does not make sense for 
a choice of n and q, then this term does not appear in the expansion. 

Proof. Each time, one computes the expansion in irreducible representations of ]/ ^ y^mia ^]^g (.g^^ 
of real- valued spherical functions, and of y-^m" (gj^-^min in the case of complex- valued spherical functions, 
where A i-> A* is the transformation of weights given by Equation (4.6). This expansion can be found 
with Schur functions; let us detail for instance the case of complex Grassmannian varieties Gr(n, (7,C). 
With an alphabet of eigenvalues Z — {2:1, . . . , 2;„} such that Z1Z2 ■ ■ ■ Zn = 1, one has 

S(0,...,0)„_i(^) = 1 
S(2,l,...,l)„_i(^) = S(i,o,...,0,-l)„(^) = E ^^^7^ ~ ^ 



,h3 = l 



S(4,2,...,2)„„i(-Z') — S(2,o,...,0,-2)„(-Z') — I E E ^«^J^/c 

l<i<i<n l<A;<;<n 



E- 

i,«,i=i 
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S(2,2,l,...,l,0)„_i(^) = ■S(l,l,0,...,0,-1,-1)„(-^) = I 

, l<i<j<n l<fe<;<n 



-1 ,-1 











)-l 













S(3,l,...,l,0)„_i(^) = S(2,0,...,0,-l,-l)n(^) = X ^i^j^k ^-^1 M + 1 

\l<i<j<nl<k<l<n J \i,j=l J 

S(3,3,2,...,2)„_i(^) = S(l,l,0,...,0,-2)„(^) = X X ^i^J^k^^i^ ~ ^ ^^^^"^ ^ ^" 

\l<i<j<nl<k<l<n J \i,j=l J 

Consequently, 

^^(2,l,...,l)„_i-j®2 ^ y(0,...,0)„_i ^ 2 V'^^'^'-'^^"-! y(4,2,...,2)„_i ^ y(2,2,l,...,l,0)„_i 
^ y(3,3,2,...,2)„_i ^ y(3,l,...,l,0)„_i^ 

because the same equality with Schur functions holds. The second line corresponds to non spherical 
representations, so only the terms of the first line can contribute to ((^(^'I' - 'i)"-!)^. Entirely similar 
calculations yield the following expansions: 

• Gr(n,g,M): 

j^y(2,0,...,0)L5j = y(0.-.0)LtJ e y(2,0,...,0)LtJ V^''''°-'°^^ii 

^y(2,2,0,...,0)L5j ^y(3,l,0,...,0)L5j ^y(4,0,...,0)L5j 

• Gr(n,g,H): 

^y(l,l,0,...,0)„^®^ ^ y(0,...,0)„ ^ y(l,l,0,...,0)„ ^ y(l,l,l,l,0,-..,0)„ ^ y(2,2,0,...,0)„ 
^ y(2,l,l,0,...,0)„^ 

Only the terms on the first line are spherical. 

• S0(2n,M)/U(n,C): 

^y(l,l,0,...,0)„~j®^ ^ y(0,...,0)„ ^ y(l,l,0,...,0)„ ^ y(l,l,l,l,0,...,0)„ ^ y(2,2,0,...,0)„ 
^ y(2,0,...,0)„ ^ y(2,l,l,0,...,0)„^ 

again with non-spherical representations gathered on the second line. 

• SU(n,C)/SO(n,K): 

y(2,0,...,0)„_i ^ y(2,2,...,2)„_i ^ y(0,...,0)„_i y(4,2,...,2)„_i ^ y (2,l,...,l)„-i ^ 

and the last term is not a spherical representation. 

• SU(2n,C)/USp(n,e): 

y(l,l,0,-.,0)2„-l ^ y(l,...,l,0)2n-l ^ y(0,...,0)2„-l ^ y(2,2,l,...,l,0)2n-l y(2,l,...,l)2„-l^ 

and again the last term is not spherical. 

• USp(n,H)/U(n,C): 

^y(2,0,...,0)„~j®^ ^ y(0,...,0)„ ^ y(2,0,...,0)„ ^ y (2,2,...,0)„ ^ y (4,0,...,0)„ 

The terms on the second line corresponds to non-spherical representations. 
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As mentioned before, the coefficient of the constant function in j^"^""!^ is then always equal to jjx—- O 



For the spaces SU(n)/SO(ri) and SU(2n)/USp(n), the remaining coefficient a can be found by evalu- 
ating the spherical functions at eg. Thus, 



SU(7i)/S0(n) : 
SU(2n)/USp(n) 



2,0,...,0)„_i 

;i,i,o,...,o)2„-i 



' _ 2 ^ ,(4,2,...,2)„_i. 



+ n 11^ + n 
2 1 2n2 - n - 1 



(2,2,1:...,1,0)2 



But in the other cases, the values of the spherical functions appearing in the right-hand side of the formulas 
of Lemma 2G are unfortunately not known a priori, which makes finding the coefficients a,b,c, . . . quite 
difficult. However, since one only needs to compute Et[(0'^""")^], and since c^-'''""' is explicit in terms of 
matrix coefficients, one can use the following Lemma (c/. [Levll, Proposition L4]). 

Lemma 27. Let k > 1 be any integer, and (gf)tgR^ be the Brownian motion on SO(n) or SU(n). The 
joint moments of order k of the matrix coefficients of gt are given by 

E[5n=exp(t^(J„r-+t ^ ri,,,{CM (4.8) 
\ i<i<]<k ) 

where is the coefficient introduced on page 10; Cg is the Casimir operator; and rji j is the linear map 
from M(n, k)*^^ to M(n, k)'^'^ defined on simple tensors X ®Y by 

In the complex case, one has also: 

\ l<i<j<k+l 

with 

7j,^j{X (g>Y)^ I -Vt,3{X (E) r*) ^/^ G [1, kj and j e {k + 1, k + Ij ; 

[r,,,,(x*(^r*) tft,jelk + i,k + ij. 

Proof. In the complex case, recall the stochastic differential equation satisfied by gt, and therefore by g^: 

dgt = gt dBt + ^gtdt ; dg; = -gt dB^ + ^gidt. 
Ito's formula yields then the following stochastic differential equation for (gt)'^'' <E) (fft)*^': 

d{g®' ® 5^'), = igtr' ® (5^)^' I ^±±1^ dt+ ^^A^Bt ® dBt 

y l<i<j<k 

Ck k+l 
Y,{Inf'-^® dBt {In)'®'-'®dBl®(Uf 
4=1 i=k+l 



The quadratic variation of Bt is given by the Casimir operator: dBt ® dBt ~ Cg dt. Taking expectations 
in the formula above leads now to a differential equation for £[((74)®'^ ® (fft)*^'], whose solution is the 
exponential of matrices in the statement of this lemma. The real case is the specialization Z = of the 
previous discussion, though with a different Casimir operator. In the quaternionic case, one has to be 
much more careful. In particular, since the quaternionic conjugate of pq is q*p* instead of p*q* , in the 
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previous argument the SDE for does not take exactly the same form. A way to overcome this problem 
is to use the doubling map (1.2). Thus, we write an equation for gt instead of gt'. 

\ l<i<j<k J 

where the Casimir is now considered as an element of (End(C^"))'*^. As we shall see later, joint moments 
of the entries of g and g* are combinations of the joint moments of the entries of Ij, so the previous 
formula will prove sufficient to solve our problem in the quaternionic case. □ 



It turns out that in each case important for our computations, the matrix X)i<i<j<fc '^i.jC'-^g) '^^'^ t)e 
explicitly diagonalized, with a basis of eigenvectors that is quite tractable (to be honest, with the help 
of a computer). In the following, we describe the eigenvalues and eigenvectors of these matrices, and 
leave the reader check that they are indeed eigenvalues and eigenvectors: this is each time an immediate 
computation with elementary matrices, though quite tedious if A: = 4 or A; + Z = 4. In what follows, we 
write e[ii, 12, . . . , v] = e^^ (Ki e^^ (g) • • • (g) e^,,. 



4.1.2. Quotients of orthogonal groups. For special orthogonal groups, set ^^M^^k = Si<i<j<fe Vi,j{C^so{n))i 
viewed as an element of End((M")'^'^). If fc = 2, then the eigenvalues and eigenvectors of A/„ 2 — 
Ei<i<,<niEij - E,ir^ are: 



eigenvalue 


multiplicity 


eigenvectors 


n — 1 


1 




1 


n(n— 1) 
2 


- * < j 


-1 


(n+2)(n-l) 
2 


+ e[i,«], i < j 
e[i,i\ — e[i + l,i + 1], i<n—l 



This allows to compute exp(— ^" „"'"''* ) exp(|-M„_2), which is the right-hand side of Formula (4.8) in the 
case of SO(n,E) and for k = 2. One obtains: 

E[{guf]^-+(l--)e-' ; E[(5,,)2] ^ A (l - e"*) ; 
n \ n J n ^ ' 

= \ (f' + c-"^*) ; E\g^,g,^\ = ^ (e"* - e"^*) ; 

and all the other mixed moments vanish (e.i?., W\giigij\ or W\gijg-ki\). Now, if fe = 4, then the eigenvalues 
and eigenvectors of M„^4 are: 



60 



PIERRE-LOIC MELIOT 



eigenvalue 


multiplicity 


eigenvectors (not exhaustive, some repetitions) 


2n-2 


3 




n 


3n(n — 1) 


T.2=i4hj,k,k]- e[j,i,k,k], i < j, * 


n-2 


3(n + 2)(n- 1) 


Efe=i «, fc, fc] - e[i + 1, i + 1, A;, fc], i < n - 1, * 


6 


n(n-l)(n-2)(7i-3) 
24 


Eaeei^Wef^^i-fc^^]"- i<j<k<l 


2 


3n(n+2)(n-l)(n-3) 

8 


Dl{i,j,k,l), D^^{i,j,k,l), D^,{i,j,k,l), i^j^k^l 





n(n+l)(n+2)(n-3) 
6 


Si{i,j,k,l), 82(1, j,k, I), i^j^ky^l 
Ki{i,j,k,l), K2{i,j,k,l), i j ^ k ^ I 


-2 




/ ^'1^2 1 „r- t,l®2 1 ^[1, -7(2)2 \ 

/ e\i,j\^ -\-eyj ,K,\^ -\-e[K \ • / • / 7 
\-elj,i]»^-elk,j]«'''-eli,k]'^y ' ^ 3 T * 

Dl{i,j,k,l), Dl{i,j,k,l), Dliz,j,k,l), i^j^k^l 


8 


-6 


n(n-l)(n+l)(n+6) 
24 


Eaee^4hj,k,l]'', i<j<k<l 
e[i,i,i,i]+e[j,j,j,j]-Y^'^^Q^e[i,i,j,j]'', i <j 



The star ★ means that the eigenvectors of a basis are listed up to action of S4; and the symbols ^^664 
mean that we take the sum of all distinct permutations of the tensors. For the eigenvectors associated to 
the value 2, denote ©4,(1) = ((1, 3, 2, 4), (1, 2)), ©4,(2) = ((1, 2, 3, 4), (1, 3)) and 1)4,(3) = ((1,2,4,3), (1,4)) 
the three dihedral groups of order 4 (hence cardinality 8) that can be found inside 64. Each dihedral 
group of order 4 has for presentation 

S)4 = (r,s I = s'' = {rsf = 1), 

so the parity rj^cr) of the number of occurrences of s in a reduced writing of G 'D4 is well-defined, 

and provides a morphism -q : ^4,{v) — ^ {±1} for ^ = 1,2,3. Then, it can be checked that for every 
i j k ^ I and any v, 

D^,{i,j,k,l)= Yl v{cT)e[i,j,k,ir 

is in V2. The eigenvectors Df{i, j,kJ), D^ii, j^kd) and D^{i, j,k,l) associated to the eigenvalue —2 are 
defined exactly the same way, but with the morphism 9 : 2)4, (w) — ^ {il} associated to the parity of the 
number of occurrences of r in a reduced decomposition of a GD4 (again it is well defined) : 

Dl{i,j,k,l)= 0i<^)e[t,J,k,ir. 
fe2)4,(„) 

For the eigenvectors associated to the value 0, Si{i,j, k, I) is defined by 

e[i,i,k,k] + e[k,k,i,i] + e[j,j,l,l] + e[l,l,j,j] - e[i,i,l,l] - e[l,l,i,i] - e[j,j,k,k] - e[k,k,j,j] 
- e[i, k, k, i]-e[k,i, i, k] - e[j, I, I, j] - e[l, j, j, I] + e[i, I, I, i] + e[l, i, i, I] + e[j, k, k, j] + e[k, j, j, k] 

and 82(1, j,k, I) is obtained by replacing each term a06060a by a060a06 in the previous formula. On 
the other hand, if J?4 = {id, (1,2)(3,4), (1,3)(2,4), (1,4)(2,3)} denotes the Klein group, then Ki{i,j,k,l) 
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and K2{i,j, k, I) are defined as 

cre.ft4 CTe(i,2,3)j?4 

o-ej?4 cre(l,3,2)J?4 

That said, tlie deduction of the mixed moments of order 4 of the coefficients of g goes as follows. One 
notices that 



(n + 2)^ef = J2 e[k,k,l,l]+e[k,l,k,l]+e[k,l,l,k]] +^ 



,k,l=l 



(TeS4 



with the first sum in the eigenspace V2n-2 and the second sum in the eigenspace V-q. On the other hand, 
for any iy^j, 

(n + 4)(ef4-ef )=6e[*,z,z,*]+ ^ ^ e[i,i,k,kr - ^ ^ e[j,j, k,kr - 6e[j,j,j,j] 

+ X! I e[i,i, + e[fc,fc,fc,fc] - ^ e[i, i, fc, A:]'^ j - I e[j, j, j, j] + e[fc, /c, fc, fc] - ^ e[j,j,k,k]'^\, 

with the ffist hue in K-2 and the second in V-e- Since ef ^ = ^ ^ + ^ EJ=i(ef - ef one 

concludes that 



1 " 



fc,;=i 



— — — ^ ^ (e[i,i,fc,fc] -e[Z,Z,fc,/c])" 

^ ' ae&i k,l=l 



n+1 



(?i + 2)(n + 4) 
1 

(n + 2)(n + 4) 



5^ I e[i, i, i] + e[A:, /c, fc, A;] — e[z, z, /c, fc]'^ 



o-eS4 



I e[fc, fc, fc, fc] + e[l, /, Z, /] — 5^ e[k,k,l, 



each line being in a different eigenspace: V2n~2, Vn-2; V-q and t^_6- This leads to the following formulas^^': 

mm,)'] = 

n{9^3f{g^kf] = 



3 




6(n - 




(n 




1) 


n{n 4 


2) 


+ -) 

nyn H 




(n 


+ 2)(n-l 


4) 


3 




6 






3 






2) 






(n 


+ 2)(?i-l 


4) 


1 










(n + 1) 




n{n -f- 


2) 






[n 


+ 2)(n-l 


4) 


1 




2 






1 




n(n -f- 


2) 


n{n H 


-4)' 


[n 


+ 2)(n4 


4) 



^^The technique is the following; to compute ]E[5iji9ij29»i392j4]i counts the number of occurrences of e{ji, j2, jz, ji] 
in each term of the previous expansion. 
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and one sees also that the other expectations V,[gijgikgugim] vanish {e.g., E[(7ij5ift(5i;)^] with i ^ j ^ k ^ 



I). Similar manipulations yield the decomposition in eigenvectors of ef ^ (8) ef^: 



n + 1 " 1 " 

e[fc, k, I, I] - _ -.^ . , Yl ^[^' ^' ^' ^1 + ^[^' ^' ^] 



n(n - l)(n + 2) ^j^/ ' " ^ (n - l)n{n + 2) 

1 / " " \ 

+ / Y Xl^t^'^'^'^J ~ e[fc,A;,Z,Z] + ^e[Z,/,j,i] - e[l,l,k,k] 

^'^^ ^1 k^^,J \l=l 1=1 ) 

^ f / n n \ 

' k^i,j aeei\i=i 1=1 J 

+ ^7 ^ VI Si{i,k,j,l) + Si{i,l,j,k) + S2{i,k,j,l) + S2{i,l,j,k) 

bin — l)(n — 2) 

+ ^ Y + e[i, +e[A;,fc,i,i] - e\j,j,i,i] - e[k,k,j,j] - e[i,i,k,k] 



^ 6(n 



^ f E + er-t: e[i, i, k, kA + E (ef + ' t e\J,j, k, kA ] 

-5 k + t -M."l')- 3(n + 2K„ + 4) i: k + t *M,<A 

\ 0-664 / ^ ^ (fe<o \ <Tee4 / 

where the eigenspaces associated to each part are V2n-2, Vn-2, Vq, V-2 and V-q. As a consequence, 

.0, .0. n + 1 2(n + 3) _f n — 3 2^-2 ^ n — 2 

^[(^'^^ ] = (n-lMn + 2) + ^^kr^ ^ + ^ " + ^ ^ 

+ 4n + 6 2n+4 . 

-I e n 

6(n + 2)(n + 4) 

] = n(n-l)(n + 2) " ' + W^) ' ^ 

n^ + 4n + 6 _2n+4j 
^ 6(n + 2)(n + 4)^ " ' 
r, \2/' \2i _ n + 1 n^ - 8 _( n — 3 2^-2 ^ 

H[9ii) {9jk) J - n{n-l){n + 2) + n(n-2)(n + 4) " 3(n-l)(n-2) ^ " 
_Jl^^_2t_ n _2»+4, 



2n 6(n + 2)(n + 4) 

. 2/ \2i n + 1 2 _j n — 3 2ra-2 ^ 

mgij) im) J = + 2) " (n-2)(n + 4) " 3(n-l)(n-2) " 

1 Of n 2n+4^ 

H e~ e — ~ • 

2n 6(n + 2)(n + 4) 
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n 



Hgttdijdjjgjt] = - TTT — ^ — 7 — — TTC 



Ti 


1 1 

1 -L 




•,{n — 


l){n + 

1 

L 


2) 




f 2)(71 

1 

L 


+ 4)' 


-— — 

TLyTL - 


- l)(n 
+ 4r^^ 


-X- 0\ 

+ 

^6 


6(n 


f 2)(ri 
1 


+ 4) 


n[n - 


-l)(n 
1 


+ 2) 




f 2)(ri 


+ 4)' 



2(71 + 2) _^ 



e - 



2 _f 71—3 2n-2 ^ 

- * e " ' 



(n-1) 



E[ffifc5*i5jfe5ii] = ^ TTT + —( ^TT — 7T 6 - TT, rvy — C 



2n-2 



e " . 



Finally, the elementary tensor (E) Cj (S) e^- (8) e; with i ^ j ^ k ^ I can be expanded as a combination 
of eigenvectors in Vg, V2, Vq, VL2 and V-^- This expansion is related to the following remarkable identity 
in the group algebra C(34, which can be considered as a relation of orthogonality of characters, but that 
only involves one-dimensional representations. Denote 



o-e£>4,(i) 

and similarly for ZJj, -D3, D\, D\ and D\. We also introduce / ~ Tliae&t E ~ Ylaeei ^^^) 

^i-E - - E - ; ^^-E - - E - 

cre.ft4 cre(l,2,3).S4 o-g.t?4 cre(l,3,2).ft4 

As explained before, all these sums correspond to eigenvectors in Vg, ¥2-, Vq, V-2 and V-q. Then, 

id[Ml - ^ ^ + 8 (^1 + ^2 + ^3') + ^{K, + K2) + l {Dl + Dl + Dl) + ^^E. (4.9) 

Indeed, ^ / + ^ i? is one twelfth of the sum of all even permutations in the alternate group 2t4. By 
adding j2^^^ ^ -^2) to this quantity, one removes all the permutations that are in 2I4 and not in .S4, so 



— / + — (Xi + i^2) + — £^ = - y ^T. 

24 12^ ^ ^'^ 24 4 ^ 

cre.l?4 



On the other hand, 



\ {Dl + D<i)^\ idii,4i + \ (1,2) (3, 4) - \ (1,3) (2, 4) - \ (1,4) (2, 3) 

and similarly for the two other dihedral groups. Thus, the whole sum on the right-hand side of (4.9) is 
indeed id|i.4|. As a consequence, e[i,j^k,l] is equal to 



cre©4 

^\{Dl{i,j,k,l) + Dl{i,j,k,l) + Dl{i,3,k,l)) + ^ E 



24 
o-ee4 
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with each term respectively in Vq, V2, Vo, K_2 and V-q. This leads to 

r 1 1 2n-8 ^ 1 2.1-4 ^ 1 1 2.1 + 4 ^ 

E[5^J5Jfc5fc/3^^] = ^ + 8 ° 8 24 ^ ' 

1 2n-8 ^ 1 2.1-4 ^ 1 1 2.1 + 4 ^ 

^[9ti9339ki9ik] = e " - g e " + g e + — e - *; 

1 2.1-8 ^ 1 2n-4 ^ 1 2.1-2 J 1 1 _2ii+4, 

'^[9^39]^9kl9lk] = " -gC " +ge " -g^ +^e 

and we are done with the computations in the case of special orthogonal groups. 

Remark. It would be very nice to find the analogue of Equation (4.9) for symmetric groups of larger 
order, in connection with the diagonalization of M„.2fc- An interesting feature of this identity is that it 
involves only generating functions of one-dimensional characters of subgroups H C 64: 



crG//C64 

Indeed, this is obvious for most of the terms, and notice that Ki + K2 = S(^4,x) + '^(2t4,X^)i where 
X is the quotient map 2I4 — >■ 2I4/.S4 ~ Z/3Z = {1, e^''^^'^, e"^''^/'^}. This combinatorial property of the 
reduction of M„.2fe seems profound and quite mysterious for the moment. 

Proposition 28. For the real Grassmannian varieties Gr(rt, g,M) and the spaces S0(2n)/U(n), the 
coefficients of Lemma 26 are the following: 



^ , . 2 2nV 1 1 

Gr(n, M) 



v? + n-2 3 V ("--!)("■- 2) pq{n - 2) 



,(2,2.0,....0)l^j 



(n-2)(n + 4) 3 \{n + 2){n + A) pq{n + 4) 

S0(2n)/U(n) : + + \ ^i^^-^o...,). . 

2n^ — n 3n (3n)(2n— 1) 



Proof. One expands the square of the sum given by Proposition 25, and one gathers the joint moments 
of the coefficients according to the possible identities between the indices. For real Grassmannians, 
^^(2,0, ...,0)^™ J ^ has for expansion: 

{\ E ^9.3? + \ E C9».)^ - 1) - f ^ E (^.)^ + \ E - 2^(^'--°)^t^ - 1 

\ iJ<P i,3>P J \ 'h3<P i-3>V J 

= {-) T[{9^^f] +U--\ m.^g^^f] + f 2 - -) (4T[(.gn5i2)'] + T[{g^2f] + r[(ffi252i)']) 

/ An \ f An \ 

+ 4n - 16 + — r[(5n523)'] + 2n - 12 + — {T[{g^29i:i?] + ^[(312.923)']) 

„2 _ g„ + 24 - ^"j T[(.9i2534)^] - 2 0(2^O-'O)LtJ - 1. 



where by T[{gii)'^] we mean a linear combination of products (ga)'^, whose expectation is therefore 
E[(gii)*]; by T[{giig22)'^] we mean a linear combination of products {gugjj)^ whose expectation will be 
E[(ffiig22)'], etc. 
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, 2 



Thus, the expectation of ^,^^^^°' ' '°''l3j ^ jg 



1^ - 16 2nV 1 1 



n2 + n-2 (n-2)(?i + 4) 3 \{n~l){n-2) pqin - 2) 



e - + — I 7 7T7 ^ I C " * 



1 2 



(n + 2)(n + 4) pq{n + 4) 



e " \ 



n 




1 




n 




n 




1 




n 




n 




1 


n 



and by identifying the Casimir coefficients of the spherical functions, one deduces from this the expansion 
of the square of the discriminating zonal function in zonal functions. 

For the spaces SO(27i)/U(7i), one computes again the square of the homogeneous polynomial of degree 
2 given by Proposition 25: 

1 / " V 

I X] 92i,2j 92i-l,2j-l - g2i,2j-l g2i-l,2j j 

= - (T[(5ii.922)^] +^[(512321)^] - 2T[5ii5i252252i]) 
(2T[(5i2.g34)^] -2T[gi3gi4g23924]) 
(2T[gi252i.934543] - 2 r[gi2g23g34g4i]) 

(^[311322333544] +^[512321334343] - 2r[3ll322334343]) 

+ remainder, 

with the same notations as before, and where the remainder is a combination of products of coefficients 
whose expectation vanish under Brownian (and Haar) measures. More precisely, terms with a certain 
symmetry cancel with one another when taking the expectation: for instance, 

(321, 2j 32i-l,2j-l — 92i,2j-l 32i-l,2j) X {g2k,2l g2k-l,2l-l — .92fc,2;-l 32fe-l,2i) (4-10) 

with i ^ i ^ k ^ I IS equal toa + fe — c — d, where a, &, c, d are products of type gijgki9mngopi £^nd have 
therefore the same expectation. Consequently, every product of type (4.10) will not contribute to the 
expectation of (0''^'^'°'''°''Lt j rpj^g following sets of indices have the same property of "self-cancellation": 

; ; ; ; {i,i,j,k) ; {j,k,i,i) ; 

{hj,i,k) ; {j,i,k,i) ; {i,j,k,i) ; {j,i,i,k) ; {i,j,k,l) ; {i,j,k,l) ; 

so it suffices to consider products with sets of indices ihjjjji) or — these 

are the four lines of the previous expansion. Using the formulas given before for the joint moments of 

the entries (beware that one has to use them with the parameter 2n), we obtain: 

,n 1 n n1 n9i 1 n — 1 2n-4, 4(71^ — 1) 2,1-1 J 

n{2n-l) in (3n)(2n-l) 
and it suffices then to identify the coefficients of the negative exponentials. □ 



4.1.3. Quotients of unitary groups. For special unitary groups, set -i^Mn^k,i — J2i<i<j<k+i''li,j{C^su(n)), 
viewed as an element of End((C")®*'''^'). li k — I — 1, then the eigenvalues and eigenvectors of M„ i i = 

i/„ (E) iln + n I]"j=i ^ij ® a-re: 
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eigenvalue 


multiplicity 


eigenvectors 


n2-l 


1 




-1 




e J 1 < i 7^ J < « 
e[i,i] — e[i + l,i + 1], i<n—l 



Consequently, 



n\9u\'] = 



n 



l--le-* ; E[|ft,f ] = ^ (1 - e-*) ; E[fe^] = e"*. 



If fc = Z = 2, then the eigenvalues and eigenvectors of -M„_2,2 are: 



eigenvalue 


multiplicity 


eigenvectors (not exhaustive, some repetitions) 


2n2 -2 


2 




n2-2 


4(n + l)(n- 1) 


YX=i ^[*' ^] - e[i + 1, fc, i + 1, fc], i<n-l 
YJk^i e[k,i,k,i] - e[k,i + l,k,i + l\, i <n-l 
J2k=i ^[*' k,k,i]- e[i + l,k,k,i + l\, i <n-l 
J2k=i e[k,i,i, k] - e[k,i + l,i + l,k\, i <n-l 


2n- 2 


n^(n+l)(n-3) 

4 


Ae[^j]-eb-.j])®2-(eb-,fc]-e[fc,i])®n ■ / • / / / 
\+(e[k,l]-e[l,k])«'^-(e[l,i]-e[i,l])»^J ' J T >■ 


-2 


(rt+2)(n+l)(n-l)(n-2) 
2 


I e[i,j,i,j]+e[j,fc,j,fe]+e[fc,i,fc,i] \ i ^ A ^ U 
\^-e[j,i,j,i]-e[fe,j,fe,j]-e[j,/s,i,fe]y ' '' ^ J 1^ 

( e[i,j,j,i]+e[j,k,k,j]+e[k,i,i,k] \ A ^ A ^ U 
\-e[j,i,i,j]-e[k,j,j,k]-e[i,k,k,i]) ' * ^ J ^ 


-2n- 2 


n^(n-l)(n+3) 
4 


e[i, i, i, i] + e[j,j,j,j] - {e[i,j] + e[j, i])®^, i < j 



For the eigenvectors associated to the value 2n — 2, we shall write 

S{z,j,k,l) = {e[i,j] - e\j,i]r^' - {e[j,k] - e[k,j]f^ + {e[k,l] - e[l,k]r^ - {e[l,z] - e[i,l]r\ 

Again, we can use the previous table to decompose some elementary 4-tensors in eigenvectors of M„^2,2- 
Thus, is equal to 

^ \^ p\k 1 k n-i-plk 1 1 k\ \ ( {e[i,l-:i,l\-e[k.l,k,l]) + (e[l,i,l.i]-e[l,k.lM)\ 

n(n + 1) ^ ^ ' 'J ^t'^' '^J n(n + 2)^^ \,+(e[»,z,(,i]-e[fe,;,z,fc])+(e[i,«,i,;]-e[i,fe,fe,i]) J 

^ ' fe,J=l ^ ' k^i J=l 



k<l 
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with the two first terms respectively in V2„2_2 and V^2_2, and the second Une in V-2n-2- It follows that 



n(n+ 1) 


+ 


n{n + 2) 


2 




4 


n(n + 1) 




n(n + 2) 


1 




n- 2 


n(n+ 1) 


+ 


n(n + 2) 


1 




2 



n(n+l) n(n + 2) c», _l i ^/■^ _l o\ 

1 2 
uytji J n(n + 2) 



(n 


+ l)(n 


+ 2) 




2 




(n 


+ l)(n 


+ 2) 




n 




(n 


+ l)(n 


+ 2) 




1 




(n 


+ l)(n 


+ 2) 



Similarly, Ci (g) (g) is equal to 



^ n n ^ n n 

^ ^ e[fc, Z, fc, I] - E E ^' ^' ^1 



1 ( e[i,i,i,i]-e[i,i,j,i] \ ] 1 (S^ { eli,l,i,l]-elk,l,k,l] \ 

-1-2) I Z-' v+^t''-''''-'l-«['''''''V i (n-2)(n + 2) \+e[i,3,i,3]-e[i,k,i,k]j 



n(n-2)(n + 2) ^^^^ 



I eli,l,l,i]+e\j,l,l,j]-2e[k,l,l,k] \ 

I Z-^ l^+e[Z,i,i,i]+e[(,jJ,i]-2e[i,fe,fe,i]y 



+ 4(n-l)(n-2) E 2^(*' J' ^' " ^('' ^' J' 

(fc</)7^i,j 

^ \ ^ / e[i,j,i,j\+e[j,k,j,k]+e[k,i,k,i] \ 
\fl 2-^ V-eb'>»>i>»]-«[fe>i>'=ij]-e[»,fc,i,fe]y 



k^id 



" + V^^* '^^i I 

- \ (ef + ef - (e[i, j] + e\j, i])^^) - ^^^^^^^^^^^ f E + " (^[^' ^1 + ^[^' 



^40 



with the parts of this expansion respectively in V2„2_2, V^2_2, V2„_2, 1^-2 and V-2n-2- This yields 
Ffin |2|„ |2i _ 1 I 2(n + l) t n-3 



n — 2 _2« . n + n + 2 



2ti + 2 , 



2n ^ 4(n+ l)(n + 2) ^ 

n — 2 71^+71 + 2 2,1+2 ^ 

~ 2n ^ 4(n + l)(n + 2) ^ " ' 
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i2i i2i 1 — 2n — 2 _j n — 3 2,1^2 ^ 

J " (n-l)(n + l) + n(n-2)(n + 2) ~ 4(n-l)(n-2) " 

1 „2f n — 1 2n + 2 ^ 

~ 2n ^ ^ 4(n+ l)(n + 2) ° " ' 

i2| |2i 1 2(n-l) t n-3 -2:i:^t 



E[l5'yl Iffjfcl ] - _ ^ ^) _ 2)(n + 2) 4(n - l)(n - 2) ^ 

1 n — 1 _2n±2^ 

"^271^ ~ 4(n+ l)(n + 2) ^ " ' 
^['^-I'l^-I'] = (n-lKn + 1) ' (n-2Kn + 2) °" + 2(. - iKn - 2) ^""^^ 

2^j 



2(n + l)(n + 2) 



e 



Proposition 29. For the symmetric spaces with isometry group SU(n) or SU(2n), the coefficients of 
Lem^ma 26 are the following: 

n 2 

Q 2/1 1 \ 

71^ — 1 — 4 2 \(n— — 2) pq(n — 2) 

n2 / 1 1 



2 +!)(" + 2) M(n + 2) 



.2,...,2)„_i. 



SU(n) /SO(n) • ^ + " ^ w,(4,2....,2)„-i . 



SU(2n)/USp(n) 



1 2n2 - n - 1 



(2,2,1, ...4,0)2 



2n2 — n 2n^ — n 



Proof. For SU(n)/SO(n) and SU(2n)/USp(n), the only missing coefficient has already been computed. 
For complex Grassmannians, ((^(2.1.- •■.i)n-i)2 j^g^g exactly the same expansion as in the real case, but 
with square modules. From the computation of the joint moments E[|(7ijgfc/p] performed previously, one 
deduces that the expectation of the square of the discriminating zonal function is 

n2-l n2-4 2 \{n-l){n-2) pq{n - 2) 

/ 1 1 

Y \ {n + l){n + 2) ^ pq{n + 2) 



2,1+2 f 

c 



whence the expansion in zonal spherical functions by identifying the coefficients. □ 



4.1.4. Quotients of symplectic groups. Finally, set ^Mn.k — X]i<i<j<fc+/ Vi,j{C^usp{n)), which is considered 
as an element of End((C^")®'^). Recall that the diagonalization of these matrices will yield the joint 
moments of the entries of g, the matrix obtained from g by the map (1.2). Again, as a warm-up, let us 
compute the joint moments of order 2. If fc = 2, the the eigenvectors and eigenvalues of Af„^2 are: 
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eigenvalue 


multiplicity 


eigenvectors 


2n+l 
2 


1 


ELie[2i- l,2i] -e[2z,2i- 1] 


1 
2 


{n- l)(2n + 1) 


(e[2i - 1, 2i] - e[2i, 2i - 1]) - {e[2i + 1, 2i + 2] - e[2i + 2,2i + 1]), i<n-l 
e[2i ~ 1, 2j - 1] - e[2j - 1, 2i - 1], e[2i, 2j] - e[2j, 2z], 1 < ^ < j < n 
e[2^-l,2J]-e[2J,2^-l], l<z^j <n 


1 
2 


n{2n+ 1) 


fife (8i 6/ + e; (g) e/c, 1 < fc < ^ < 2n 



It follows that 

1 rj — 1 1 

E[|feP]=- + e-* ; E[|5,,f] =-(l-e-') V*,je[l,nl; 

E[(5,,)2] = e-^* ; E[(5,,)2] ^ Vz,j e 11,2^1 . 

Now, for fc = 4. it is a little more tedious than before to find a complete list of "simple" eigenvectors of 
■Wn,4 (or at least a sufficient list to expand simple tensors). The list of possible eigenvalues of Af „ 4 is 

{2n + l,n + l,n,3,l,0,~l,-3}, 

and on the other hand, one can easily identify a basis of V2„+i: it consists in the three vectors 

n 

El e[2i-l,2i,2j"-l,2j]+e[2i,2i-l,2j\2j-l] \ 
\^-e[2i,2i-l,2j-l,2j]-e[2i-l,2i,2j,2j-l]J ' 

2,i = l 
n 

El e[2i-l,2j"-l,2i,2j]+e[2i,2j,2i-l,2j-l] \ 
\^-e[2i,2j-l,2j-l,2j"]-e[2i-l,2j,2i,2j-l]J ' 

n 

E/ e[2i-l,2j"-l,2j\2i]+e[2i,2j,2j-l,2i-l] \ 
\^-e[2i,2j-l,2j,2i-l]-e[2i-l,2j,2j-l,2i] j ■ 

But then, it becomes really difficult to describe the other eigenspaces. However, one can still find^' the 
eigenvector expansion of simple tensors such as ef"^, ef^e^^, or e[i, j, fc, l\\ hence, in the following, we just 
give these expansions (again it is easy to check that each part of an expansion is indeed an eigenvector) . 
The tensor e[«,i,i,i] is an eigenvector in T^_3, so one obtains readily 

E[(5n)1 = e-"^* ; ^{{mf \ = ^{^9.^~9^if \ = K(9^J9^kf] = 0. 
On the other hand, e[2i — 1, 2i — 1, 2j, 2i] decomposes into the following eigenvectors 

1 , , 

2n(2r.+ l) ("2n+l,2+«2„+l,3) 



V2n+1,2 — 



^^2n+l,3 — 



To guess what these expansions were, we simply computed (Af„ 4)'" v for a sufficiently big number of distinct powers 
r of the matrix. 
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n-2 



4n{n + 1) 



1 

+ 



/ e[2i-l,2j-l,2i,2j]+e[2i,2j,2i-l,2j-l] Y 
/ . / . \-e[2i-l,2j,2i,2j-l]-e[2i,2j-l,2i-l,2j]J 



1 \ " \ ^ / e[2j-l,2k,2j,2k-l]+e[2j,2k-l,2j-l,2k] ^ 

4n(n+l)^ ^ \-e[2j-l,2k-l,2j,2k]-el2j,2k,2j-l,2k-l]J 



, n—1 / 2e[2i-l,2i-l,2i,2i]+2e[2i,2i,2i-l,2i-l]-e[2i-l,2s,2i-l,2i]\ 
2n{n+ 1) \-Bl'^h^i-'^''^i>'^i-'^]-B[2i,2i-l,2i-l,2i]-el2i-l,2i,2i,2i-l]j 

i \^ \^ e[2i-l,2j-l,2i,2j]+e[2i,2j-l,2i-l,2j] V 

<l)Z^Z^ {-e[2i-l,2j,2i,2j-l]-el2i,2j,2i-l,2j-l]J 



2n- 1 

+ 



■Sr^ / e[2i-l,2j,2j-l,2i]+e[2i,2j-l,2j,2i-l] Y 

/ J / J \-e[2i-l,2j-l,2j,2i]-e[2i,2j,2j-l,2i-l]J 



+ 



2(2n+ l)(2n + 2) 

1 / e[2j-l,2fc-l,2j,2fc]+e[2j\2fc,2j"-l,2fc-l] 

2(2n + l)(2n + 2) ^ \-el2j-l,2k,2j,2k-l]-el2j,2k-l,2j-l,2k] 

{2n — l)(2n — 2) / 2e[2i-l,2j-l,2i,2i]+2e[2i,2i,2i-l,2i-l]-e[2i-l,2j,2i-l,2i]\ 
6(2n + l)(2n + 2) V-e[2i,2i-l,2i,2i-l]-e[2i,2i-l,2i-l,2i]-e[2i-l,2i,2i,2i-l] J 



n—1 

+ (e[2i - 1, 2i - 1, 2i, 2i] - e[2i, 2i, 2i -l,2i- 1]) 

, 1 \^ e[2i-l,2j,2j-l,2i]+e[2i,2j,2j-l,2i-l] Y 

A{n+l) ^ ^ V-«P«-l,2j-l,2j,2i]-e[2i,2j-l,2j,2i-l]j 



+ fi 51 e[2i-l,2i-l,2i,2i]^ 



6 

(T6S4 

with the parts of this expansion respectively in V2n+i, Ki+ij K)) l^-ij and ^-3. In these expansions, 
S = (Z/2Z)^ denotes the group of permutations {id, (1, 2), (3, 4), (1, 2)(3, 4)}. As a consequence, 

~ n2i 1 n — l_f 1 _i!±it (2n— l)(2n — 2) 2n+i ^ 

m.92i-i,2i-i 92i,2i) = —77. — — 7T H — 7 — —TV e H — - e " + --- — — — rr e 

n(2n + l) n(n + l) n+1 3(2n + l)(2n + 2) 

n—1 2n+2 J 1 2n+4f 

-I e — ~ -I — e — ~ • 

^2(n+l) ^6 

itnr/~ ~ \2i 1 »^-l _t 1 -^!±it (2n-l)(2n-2) 2r.+i ^ 

Jc- (52i-l,2i52i,2i-lj = — 77^ "I ? TT^ " + oTo . iN/o ^ 

n(2n + l) n(n + 1) n+1 3(2n + l)(2n + 2) 

n—1 2n+2 . 1 2n+4 . 

e ~* + - e 



2n+l. 



2(n+l) 6 

E[fe-i,2,-i ^2,2,)^] = E[(^2.-i,2, ?2,2,-i)^] = ^^^^ - e-* + e- 

and the other moments of type E[(g2i-i,a92i,b)'^] vanish. The other moments of type ^[{gab9cd)'^] with 
{a, b} 7^ {2i — 1, 2i} are much easier to compute: indeed, the expansion in eigenvectors of e[2i, 2i, 2j, 2j] 
is 

1 /2e[2i,2i,2i,2i]+2e[2j,2j,2i,2i]-e[2i,2i,2i,2i]\ 1 / e[2i,2i,2i,2i] ^ , 1 „rr,- r> • 9 • r, -la 

g \^-e[2j,2i,2i,2i]-e[2i,2i,2i,2i]-e[2i,2i,2i,2i]^ 2 V-e[2i,2i,2i,2i]^ g ^-[^^^4,^7,^./] 

0-664 
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with each part respectively in Vq, V^i and F_3; and similarly for the expansions of e[2i — 1, 2i — 1, 2j, 2j] 
or e[2i — l,2i — 1, 2j — 1, 2j — 1]. Thus, assuming that {a, b} is not a pair {2i — 1, 2i} in |1, 2n|, one has 

^r/~ ~ ^9^ 1 ^n+lf 1 2n+2^ 1 2n+4 ^ 

nigaa 9bbf] = 3 e-— * + - e-— * + - e-— *; 

„r,_ _ ,0, 1 + 1 _2ri + 2, 1 _ 2n + 4 ^ 

E[(5a6 56a) ] = 36 " + 6 ' 

and the other moments of type ]E[(5fo6 Scd)^] with {c, d} ^ {a, 6} vanish. The same expansions allows one 
to compute many moments of type ^[\gij 5fe;P], namely, all those that write as ^[\gij Qikl'^]- For instance, 
since 

\9ii\'^ = {g2i-l,2i-l g2i,2i — 52i-l,2i fl'2i,2i-l)^, 

its expectation is a linear combination of the expected values of (52i-i,2i-i 52i,2i)^) {g2i-i,2ig2i,2i-i)^ and 
52i-i,2i-i52i-i,2i52i,2i52i,2i-i- This last expectation is 

1 n-1 (2n- l)(2n- 2) ^_2^j _^ 1 ^_2^t 



2n(2n + l) 2n(n+l) 6(2n + l)(2n + 2) 6 

Thus, with a few more computations, one gets 

4 _ 3 3(n-l) (2n-l)(2n-2) 

J~ n(2n+l) ^n(n + l) ^ (2n + l)(2n + 2) 

■mn i4i ^ ^ — i ^ 2n+i 

^ " n(2n+l) " n(n+l) + (2n + + 1) " ' 
wri |2i _ 2 (n- 2) _^ 2(2n - 1) 2„+i^ 

E[|5»5i,M- ^(2n+l) +n(n + l)^ (2n + l)(2n + 2) ^ " ' 

|.. 2 2 _^ 2 _ 2n + l 

J = ^2n+1) "n(n+l)'' + (2n+l)(n+l)^ " ' 



Finally, the expansion in eigenvectors of e[2i — l,2i,2j — l,2j] enables one to compute the other 
moments of type ]E[|5ij which are useful in the case of quaternionic Grassmannians. As it is more 

than one page long, we skip the eigenvector expansion and just give the result: 

TT-,ri i2i 2n—l 2 _f n — 3 2,1-2 ^ 

Mgu gjj] = —I + e + e -> 

n{n—l){2n + l) n + 1 6(n — 1) 

n — 2 9t 2n^ — n + 3 2,.+i . 

-I e~ -I e ^ • 

2n 3(n+l)(2n + l) 

TffW i2i _ 2n — 1 2 _j n — 3 2n-2 j 

H\gij gji\ \ = n{n- l)(2n + 1) ~ n(n + 1) ^ +6(n-l)^ " 

n—2 2r7,^ — n + 3 2,.+i ^ 



2n 3(n + l)(2n + l) 



c 



i2i 2n — 1 — 3n + 1 _( n — 3 2n-2 ^ 

^ " n(n-l)(2n+l) + n(n + l)(n-2) " 6(n-l)(n-2) ® 
1 _2t 2n — 3 _2n+ij 

"271*^ ~ 3(n + l)(2n+l)® " ' 
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^[\9^3 9]k\ 



2n — I 2n — 3 _f n — 3 2n-2 ^ 

■ e " 



n(n-l)(2n + l) n(n+l)(n-2) 6(n-l)(n-2) 



1 _2t 2n — 3 2,1+1 ^ 

"^2^^^ ^ 3(n + l)(2n + 1) " ' 
21 2n — 1 2n — 2 _( 1 2^-2 ^ 

^ " n{n-l){2n + l) ^ n(n+l)(n-2) ^ 3(n- l)(n-2) 

4 2„+l. 

e " . 



c 



3(n+ l)(2n + 1) 



Proposition 30. For the quaternionic Grassmannian varieties Gv(n,q,M) and the spaces USp(n)/U(n), 
the coejficients of Lemma 26 are the following: 

1 ^ 1 1 ^ /,(i",o,...,o)„ 



^ 2n2-n-l 3 \{n - l){n - 2) pq{n - 2) 



JEl j,(1^0,....0)„-i , ^ / ^ I ^ ^ ,(2,2,0,. ..,0)„. 



USp(n)/U(n) . 



(n-2)(n + l) 3 V(« + l)(2^^+ 1) + 

i , 4(n-l)(n + l) (2,2,0,. ..,o)„ , " + 1 ^(4,o,...,o)„ 

2n2 + n 3n(2n + 1) 3n ^ 

Proof. The case of quaternionic Grassmannians is again done by using the expansion on top of page 
64, with square modules instead of squares. One obtains the following formula for the expectation of 

(^(1,1,0,.. .,0)„)2. 



f„-^ t f 1 1 
+ 7 7T e"* + — 7 7T7 ^ ^ I e 



2n2-7i-l (n-2)(n + l) 3 \(n - l){n - 2) pq{n - 2) 

.2/4 1 



271-2 



211 + 1 J 

c 



3 V('^+ l)(2n+ 1) pq{n+l), 
hence the expansion in zonal functions by identification of the coefficients. Finally, for the spaces 
USp(n)/U(n), (0(2,0,. ..,o)„)2 jg ^^^^^^ to 

^ +T[(5n522)'] +^[(312321)']) + — (^[(313324)'] +T[(gn533)'] +T[(5i353i)']) 

In n 

plus some remainder whose expectation under Brownian measures will be zero. Hence, 

IE[(0(2,o,...,o).i)2j ^ 1 + A{n- l){n + l) n + 1 ^ 

^ n(2n+l) 3n(2n + l) 3n 

and 2"^^ is the exponent corresponding to the spherical representation of label (2, 2,0,..., 0)„, whereas 
^"^^ is the exponent corresponding to the spherical representation of label (4, 0, ... , 0)„. □ 

4.2. Proof of the lower bound on the total variation distance. The proof of the lower bound is 
now a simple application of Bienaymc-Chebyshev inequality. First, under the Haar measure, we have: 

Proposition 31. If Ea is the event > a}, then the Haar measure of Ea satisfies the inequality 

Vx{Ea) < 4 

for every classical simple compact Lie group X = K and every classical simple compact symmetric space 
X = G/K. 

Proof. The previous computations ensure that ]Eoo[|i^P] = 1 in every case, so 



a2- 



□ 
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Next, let us estimate Et[^] and Vart[r2] for t = a (1 — £)logn. The exact values are listed in the 
following table: 



K or G/K 




VaiiP 


SO(n) 




1 + "("2"^) e- 


\ / 


o U yfi ) 


n e 2ti^ 


l + ( 


— 1) c 




USp(n) 




l + (2n + l)( 


u — 1) c^* + (2n + 1) ne~^* — 4n^ e~^^* 


(jr(n, g, M) 


/(ri+2)(n-l) „_t 

V 2 e 




2n^ 0' 

pq 


\ (n-l)(n+2) t „2 / (n+2)(«-l)\ „-2=^t 
/ (n-2)(ri+4) ^ ' 3 \^n-2 p(/(n-2) ^ 

2(n+2)(n-l)\ (n+2)(«-l) 2t 
pq(n+4) y L " 2 


Gr(n,o,C) 


\/n^ — 1 e~* 


' 2 


2n^ 

Vn+2 ^ 


\ n^-l „-t 1 f n+1 \ 2^^^ 
j n^-i'^ 1 2 l^n-2 pq(n-2) ) ^ 


Gr(n,g,H) 


^{2n + l){n-l)e-* 


' + ( 

2 


pq J 

/4(n-l) 


(«-l)(2n+l) t „2 /2„+i (2n+l)(«-l)\ 2^j 
(n-2)(n+l) ^ 3 n-2 pg(n-2) ^ ^ 


S0(2n)/U(n) 


^n{1n - l)e-V^* 




n-l)(2n- 
3 




SU(n)/SO(n) 




1 + i 


n+2)(n 1) 2„+2j „(„+!) (^-i)(2^+4)^ 
2 " 2 ^ " 


SU(2n)/USp(n) 


, ^ (,.-l)(2r,+ l) 

V2ri^ — ne 


1 + (2n2 - n 


2,1-1, , ^ . (n-l)(2„ + l) 

-l)c —*^-{2n^-n)c --^ * 


USp(n)/U(n) 


v/n(2n + l)e-'^* 


1 1 4(«-l)(n+ 
^ ^ 3 





In the following we assume e < \] indeed, Lemma 5 ensures that it is sufficient to control the total 
variation distance around the cut-off time. We shall use a lot the inequality of convexity 

e^ - 1 

expfa;) < 1 H x Va; G (0, y) . 

V 

Lemma 32. Under the usual assumptions on n, for groups and spaces of structures (but not for Grass- 
mannian varieties), \ait[^] is uniformly hounded for every t = a{l — e) logn with e G (0, 1/4). Possible 
upper bounds are listed below: 

SU(n), SU(n)/SO(n), SU(2n)/USp(n) : 1 ; 
S0(2n)/U(n), USp(n), USp(n)/U(n) : 3 ; 
SO(n) : 8. 

Proof. We proceed case by case, and denote A((A,/i) = e~^* — e~^*. Notice that At(A,/i) < if A > /i. 
On the other hand. At (A, /i) is always smaller than 1 for A, /i > 0. 
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• SO(n): 

rr^i ^ fr. n(n-l) ^ /n-4 n-l\ /n(n + 1) ,\ ^ n-1 

n(n — 1) . / n — 4 n — 1\ , n(n — 1) 



. /n — 4n— 1\ , n(n — 1) ^-i, / at \ 



< 1 + — n log n e - 
since < 1.382 when n > 10, and ^^j^^ < f • Then, 



I 3(n-l)logn (n-3) log n 14 . , , ■■ 

< e 2n =n e < — (niosn) 



for n > 10, so Vart[0] < 1 + 7 = 8. 

• SU(n): 

Var,[f]] = A, (o, ^) + (n^ - 1) A, (^1, ^) < 1. 

• USp(n): 

<l + (2n + l)(n-l)At(^l,^^^ < 1 + 271^ e'^^^^* (e^^ - l) 
< 1 + - n log n e 2,. 
since ^ < 0.367 when n > 3, and ^'^'^ < |. Then, 

2n + l f 3 1ogn _3 4 , , . _ -, 

e 2,1 < e 2 = n 2 < - log n) 

5 

for n > 3, so Va,VtM < 1 + 2 = 3. 

• S0(2n)/U(n): 

2n-2^ , (n-l)(2n-l) ^ /^2n-4 2n-2^ ^ 4(n2 - 1) ^ /2n-l 2n - 2 



Vart[O]=AJ0, +^ AA , ]+ ' 

\ n J 3 \ n n J S 

(n-l)(2n-l) , /2n-4 2n-2\ , 271^ 2„-2. / 2t \ 

20 , „ 2ra-2 J 

< 1 + — nlog2ne " ^ 
since ^^^^ < 0.922 when 2n > 10, and < I- Since 

n — — ' 0.922 — 3 

2n-2 ^ 3(ti-1) log 2,i 1 (n-3) log 2n 3, , „ ,_-| 

e^ se 2» =_^e 2n — ^(^^ ^og 2n) 

for 2n > 10, one concludes that Vart[0] < 1 + | < 3. 
• SU(n)/SO(n): 

Var,[l^] = a/o, '(--l){n + 2)\ ^ (n + 2)(n - 1) ^/2(n + l)^ 2(n-llJn + 2)^ ^ 



n 
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SU(2n)/USp(?i): 

, / (n-l)(2n + l)\ , ^ /2n-l (n-l)(2n + l)\ 



• USp(n)/U(n): 

Vart r! = At 0, + ^- ^ A* , + A*' 

\ n J 3 \ n n J A 

4(n2 - 1) /2n+l 2n + 2\ An^ 2„+2^ 

< 1 -I At ^ 1 I „ 



n 



\ 4n^ 221 + 2, / f \ 



O , 2T1 + 2 J 



< 1 + ^ n log n e " * 



by using the same estimate on ^^fiii as in the case of USp(7i). Since 

_22i + 2j _3_log^ _3 4 , 

6 " <e 2 =n2<_ (nlogn) 

5 

for n > 3; one obtains Vart[ri] < 1 + | < 3. 

It is not possible to prove such uniform bounds for Grassmannians, because of the term c^* that appears 
in the variance. We shall address this problem in Lemma 34. □ 

Proposition 33. Denote Kx the hound computed in the previous Lemma for the variance of the dis- 
criminating zonal function VL associated to a space X. Then, 

-iiKx + 1) 

dTv(Mt,Haar)>l-^^^^. 
Proof. Assuming a smaller than m = Et[ri], if \Q — m\ < a, then |f2| > m — a. Consequently, 

Ht[\n\ > TO - a] > 1 - V[\n - m| > a] > 1 - Vart^ 
Next, take a = ^. The combination of Lemma 31 and of the previous inequality yields 



dTv(Mt,Haar) > ^it{Ea) ~ VxiE,,) > 1 - ^^^^ = 1 ^^^"^ ^ 



Since behaves as n^^, this essentially ends the proof of the lower bounds in the case of compact Lie 
groups and compact spaces of structures. More precisely: 

• SO{n): m? > so the constant c in our main Theorem 8 is 4(8 + 1) = 36. 

• SU(n): again, to^ > n^^, so the constant is 4(1 + 1) = 8. 

• USp(n): here, > 4n^^ e~^^ > n'^^ for n > 3, so the constant is ^ 4(3 + 1) = 5. 

• S0(2n)/U(n): to^ > {2nf^ > ^ {2n)^^ for 2n > 10, whence a constant ^ 4(3+1) = f < 8. 



SU(n)/SO(n): m'^ > ^ e ^ > ^ for n > 2, so a possible constant is 3 x 4(1 + 1) = 24. 



SU(2n)/USp(n): > {2n)^^ > § {2n)^^ , and a possible constant is | 4(1 + 1) = f < 22. 
USp(n)/U(n): to^ > 2n2= e"^^ > if n^^ for n > 3, whence a constant i| 4(3 + 1) = 17. 



□ 



Unfortunately, for Grassmannian varieties, the variance of O at time t = {1 — e)logn can only be 
bounded by a constant times n^. However, since the mean of fl is also of order n^, this will still ensure 
that the discriminating zonal spherical function has not at all the same behavior under Haar measure 
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and under Brownian measures before cut-ofF time. The only downside is the loss of a factor rt^ in the 
estimate of the total variation distance. 



Lemma 34. Under the usual assumptions on n, for Grassmannian varieties, 

5 i/ k = C or H, 



VartP ^ J 3 if\ 



for every t — a{\ — e) logn with e S (0, 1/4). 



Proof. The quantity ^ is bounded by 



4 1 

< — < 



pq n — 1 

the extremal values corresponding to p — q — ^ and to p = n — 1 or g = n — 1. In particular, in the 
expansions hereafter, all the coefficients that precede differences of exponentials At(A,/i) are positive. 
Now, we proceed case by case: 

• Gr(n,g,M): 

V.r,|0| . A.(0, 2, + - S) ''-'!;-:;' A,,l, 2) ^ (!1±| - ™^) A, f^, 

\ pq J (n — 2j(ri + 4) 3 \n — 2 pq[n — 2) J \ n 

fn-1 ^ 2{n + 2){n- l)\ / 2n + 4 



6 \n + A pq{n + A) J \ n 

/ 2n — 2 

< l + 2nAt(l,2) + y — — ,2 

For the difference Af(l,2), one cannot obtain a better bound than e^* = rf^^. The second 
difference At( ^"~^ , 2) is bounded from above by 



9t / 2t \ 3 8 logn „ 9 

e"2* [e^ " ^) - - 



3n 

by using similar arguments as in the proof of Lemma 32, and the inequality n > 10. So, 

2 

Yart[n] <l + -+2n^ < Sn" . 

O 



• Gr(n,g,C): 

Var417 = At 0,2 + A* 1,2 +^ ^ aJ ,2 

/n- 1 - 1 \ , / 2n + 2 
' • ^AJ^— ,2 



2 V"- + 2 pq{n + 2)J \ 2 



n 



2 



< l + 2?iAt(l,2) + yAt(^?^^,2 

The second difference is controlled exactly as in the case of real Grassmannians, but under the 
constraint n> 2: 

^ , 2n — 2 \ 9f / 2t A 3 3 log n 9 9 

Ati ,2 I < e-2* (e^ - 1 1 < — < - n"^ 

n 4 



Hence, Varf [f^] < 1 + | + 2n'^ < 5n^. 
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. Gr(n,9,H): 

Var.lQ] = A.(0,2) + - - „(, 'l 2) ) ^'( — 

/ (n-2)(n+l) 3 pg(n + 1) / \ n 

2n^ /2n — 2 \ 3 

< l + 2nAt(l,2) + — At ,2 < 1 + 2n^ + - < hrf . 

6 \ n ) 2 

□ 



Now, Proposition 33 still holds, but with Kx varying with n and equal to irf or ^rf according to the 
field k = M, C or H. Thus: 

Proposition 35. For Grassmannian varieties Gi{n,q,k), if t — {1 — e)logn with e £ (0, 1/4), then 



c?TV ifJ-t , Haar) > 1 ^ with L 



Ln^ , ^ 16 z/Ik = M, 
— TT- with L = < 

m? 1 24 z/Ik = Cor] 



Finally, the deduction of the constants in Theorem 8 goes as follows: 

• Gr(n, (?, M): m > so the constant can be taken equal to 2 x 16 = 32. 

• Gr(n, C): m > "r?^ > | so a possible constant is again | 24 = 32. 

• Gr(n, q,m): m> ^""^'^^ n"^^ > | n"^^ for n > 3, whence a constant | 24 = 16. 
These computations end the proof of the cut-off phenomenon. 
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